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1  1*.  ABSTRACT  j 

This  report  documents  the  computer  programs  of  the  Bent  Ionospheric  Model  and 
briefly  describes  the  development  of  the  model.  The  FORTRAN  Program  is 
designed  for  general  use  and  can  generate  ionospheric  data  on  a  world-wide 
basis  for  any  past  or  future  date.  For  a  given  condition  consisting  of  station, 
satellite  and  time  information,  the  electron  density  versus  height  profile  is 
computed  from  which  range,  range  rate,  and  angular  rofraction  corrections  as 
well  as  vertical  and  angular  total  electron  content  are  obtained.  The  model  has 
the  additional  capability  of  improving  its  predictions  by  updating  with  actual 
ionospheric  observations.  Considerable  tests  in  the  past  have  proved  this 
empirical  model  highly  successful.  Also  included  in  the  documentation  is  an 
alternate  version  of  the  ionospheric  program  to  be  used  when  stringent  space 
and  time  requirements  are  Imposed  by  the  operating  system.  However,  several 
options  of  the  standard  program  are  not  incorporated  and  the  accuracy  of  the 
results  is  somewhat  reduced. 
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ABSTRACT 


This  report  documents  the  computer  programs  used  in  the  Bent 
Ionospheric  Model  and  briefly  describes  the  development  of  the  model. 

The  FORTRAN  Program  is  designed  for  general  use  and  can  generate 
ionospheric  data  on  a  world-wide  basis  for  any  past  or  future  date.  For  a 
given  condition  consisting  of  station,  satellite  and  time  information,  the 
electron  density  versus  height  profile  is  computed  from  which  range,  range 
rate,  and  angular  refraction  corrections  as  well  as  vertical  and  angular  total 
electron  content  are  obtained.  The  model  has  the  additional  capability  of 
improving  its  predictions  by  updating  with  actual  ionospheric  observations. 
Considerable  tests  in  the  past  have  proved  this  empirical  model  highly  success¬ 
ful.  Also  included  in  the  documentation  is  an  alternate  version  of  the  ionospheric 
program  to  be  used  when  stringent  space  and  time  requirements  are  imposed 
by  the  operating  system.  However,  several  options  of  the  standard  program 
are  not  incorporated  and  the  accuracy  of  the  results  is  somewhat  reduced. 
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This  specification  establishes  the  requirements  for  complete  identification 
of  Items  #0001  and  0002,  "Documentation  and  Description  of  the  Bent  Ionospheric 
Model,  "  to  be  formally  accepted  by  the  procuring  activity. 

The  Bent  Ionospheric  Model  is  an  empirical  world-wide  computerized 
alogrithm  capable  of  predicting  the  ionospheric  electron  density  profile  and 
the  associated  delay  and  directional  changes  of  a  wave  due  to  refraction.  The 
following  documentation  of  this  model  is  formatted  in  accordance  with  Para¬ 
graph  60.  5,  computer  program  product  specifications,  MIL-STD-483, 
"Configuration  Management  Practices  for  Systems,  Equipment,  Munitions, 
and  Computer  Programs," 

Sections  3.  1  and  3.4  outline  the  overall  program  structure, 

Section  3.2  gives  a  detailed  description  of  each  program  component, 

Sections  3.3  and  4.  1  incorporate  the  program  operation  description, 

Section  6.  1  and  6.  2  outline  the  ionospheric  model  development  and  present 
its  accuracy  and  limitations,  and  Appendix  I  contains  the  program  listings. 


2.0  Applicable  Documents 

The  documents  of  exact  issue  shown,  form  a  part  of  this  specification 
to  the  extent  specified  herein.  In  case  a  conflict  occurs  between  the  referenced 
reports  and  the  detailed  content  of  sections  3,4,5,  and  10,  the  detailed  content 
shall  be  considered  a  superseding  requirement  for  this  CPCI. 

Reference  > 

1.  1C.  V.  Appleton  &  W,  J,  G,  Eeynon,  Proc,  Phys.  Soc,  52,  Pt.  I, 

518(194.0);  Proc.  Phys.  Soc.  59,  Pt.  II,  58(1947) 

2.  R.  B.  Bent,  S.  K.  Llewellyn,  M.  K.  Walloch,  "Description 

and  Evaluation  of  the  Bent  Ionospheric  Model,  "  Vol.I, 

SAMSO  TR-72-239  (Oct.  1972) 

3.  S.  Chapman  &  J.  Bartels,  "Geomagnetism,"  Vol,  II,  Oxford  at 

the  Clarendon  Press  (1962) 

4.  D.  C.  Jensen  &  J.  C.  Cain,  "Iterim  Geomagnetic  Field,  " 

J.  Geogr,  Res.,  No.  9.  3 56 8- 3 56 9  (Aug.  1 96 2 ) 

5.  W.  B.  Jones,  R.  P.  Graham,  M.  Leftin,  "Advances  in  Ionospheric 

Mapping  by  Numerical  Methods,"  ESSA  Technical  Report 
ERL  107. ITS  75,  (May  1969) 

6.  W.  B.  Jones  &  D.  L.  Obitts,  "Global  Representation  of  Annual 

and  Solar  Cycle  Variation  of  f0F2  Monthly  Median  1954-  1958,  " 
OT/ITS  Research  Report  No,  3  (Oct.  1970) 

7.  A.N.  Kazantsev,  Tr.  IRE  AN  SSSR,  2,  36,  (1956  ) 

8.  R.  G.  Maliphant,  "The  Refractive  Deviation  of  Radiowaves  that 

Penetrate  the  Earth's  Ionosphere,"  DRTE  Report  No.  1000, 
(Sept.  1962) 

9.  F.  G.  Stewart,  M.  Leftin,  "Relationship  Between  10.7  cm  Uttowa 

Solar  Radio  Noise  Flux  and  Zurich  Sunspot  Numbers,  1 
ESSA  Technical  Report  (Oct.  1970) 
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In  the  area  of  satellite  communications  the  refraction  incurred  u,  u. 
wave  propagating  through  the  ionosphere  is  most  important.  The  Rent 
Ionospheric  Model  is  an  empirical  world. wide  algorithm  capable  of  accurate 
estimating  the  electron  density  profile  and  the  associated  delay  and  dire  ciioi.a’ 
changes  of  a  wave  due  to  refraction.  The  model  computes  the  electron  density 
versus  height  profile  from  which  the  range,  range  rate,  and  the  angular 
refraction  corrections  for  the  wave  are  obtained  as  well  as  the  vertical  and 
angular  total  electron  content.  Although  the  model  is  presented  for  ground 
to  satellite  communications,  it  is  readily  adaptable  for  ground  to  ground, 
or  satellite  to  satellite  communications. 

The  only  required  inputs  to  the  model  are  satellite  and  station  position 
and  time  information  and  a  limited  amount  of  Bolar  data.  For  the  model's 
additional  capability  of  improving  the  ionospheric  predictions  by  use  of  actual 
ionospheric  obse rvations ,  measured  values  of  electron  content  or  the  critical 
frequency  of  the  F2  layer,  foF2,  can  be  incorporated  along  with  the  observation 
Btation  and  time  information.  This  update  option  uses  a  weighted  mean 
technique  that  can  accept,  for  the  update,  several  measurements  from 
different  stations  separated  in  time  and  space  from  the  time  and  location  at 
which  the  ionosphere  is  to  be  evaluated. 

The  updating  process  is  generally  used  for  predicting  ionospheric 
conditions  or  refraction  corrections  after  the  fact,  when  observations  are 
generally  available.  However,  the  model's  prediction  accuracy  without  update 
accounts  for  approximately  75  to  80  percent  of  the  ionosphere  which  can  improve 
with  update  to  approximately  90  percent.  Thr.  model,  therefore,  may  be 
applied  for  future  predictions  or  after  the  fact  calculations.  Since  the  model  has 
been  developed  on  a  world-wide  basis,  predictions  are  not  limited  to  any 
particular  land  mass  or  segment  of  the  world.  The  updating  technique  does, 
however,  require  that  ionospheric  observations  be  from  stations  within 
2000  km  radius  of  the  evaluation  site.  The  model  is  applicable  for  determining 
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wave  refraction  and  ionospheric  characteristics  up  to  2000  km  in  height  and 
for  all  radio  wave  frequencies  as  long  as  the  vertical  component  is  slightly 

higher  than  critical  frequency. 

Built  into  the  model  are  the  combined  influence*  of  geographical  and 
geomagnetic  effects,  solar  activity,  local  time,  and  seasonal  variations. 

These  combined  effects  are  the  results  of  an  extensive  investigation  of  a 
vast  ionospheric  data  base  that  included  over  50,  000  topside  soundings,  6,  000 
satellite  measurements  of  electron  density  and  related  foF2,  and  over  400,  000 
bottomsidn  soundings.  The  data  base,  which  formed  the  basis  of  the  model, 
extended  over  the  period  of  1962  to  1969,  covering  the  minimum  to  maximum 
of  a  solar  cycle,  For  further  information  regarding  the  development  and 
evaluation  of  the  Bent  Ionospheric  Model,  see  Reference  2  and  Section  6.  0. 


3.  1  Functional  Allocation  Description 


The  ionospheric  PROGRAM  ION  is  written  in  FORTRAN  IV  code  and  lias 
a  simple  load  structure  with  no  overlay  requirements,  The  following 
program/subroutines  comprise  the  CPCI,  and  the  attached  diagram  identifies 
the  calling  routine*  and  the  subprograms  called  for  each  computer  program 
component; 

CPCs  :  PROGRAM  ION,  and  SUBROUTINES  REFRAC,  PLOTNH,  PROFL1, 
PROFL2,  BETA,  SICOJT,  DKSICO,  MAGFIN,  GK,  DKGK. 

The  following  library  subprograms  are  required  : 

ABS,  AMOD,  AT  AN,  COS,  EXP,  LOG1Q,  SIGN,  SIN,  SQRT. 

All  internal  data  transfer  between  tho  individual  CPC*  occurs  through  labeled 
common  blocks  and  through  the  calling  sequences,  which  are  both  described 
under  Section  3,  2.  1,  3  for  each  CPC.  The  external  data  transfer  consists  of 
input  coming  from  the  data  card  deck  into  PROGRAM  ION  and  from  the  iono¬ 
spheric  coefficient  data  tape  into  SUBROUTINE  REFRAC,  and  of  output  of  the 
reaulta  from  PROGRAM  ION  and  SUBROUTINE  PLOTNH  to  the  line  printer; 
these  files  are  described  in  detail  in  Section  3,  3,  1,  The  functions  performed 
by  the  program  are  described  in  Section  3.  4  and  referenced  to  the  CPCs  to 
which  they  are  assigned. 

An  alternate  vercion  of  the  ionospheric  program  is  included  in  this 
documentation,  consisting  of  a  preprocessor  TABGEN  and  a  reduction  program 
ION  1 .  Both  programs  are  written  in  FORTRAN  IV  code,  have  a  simple  load 
structure  with  no  overlay  requirements,  are  run  as  separate  entities,  and  are 
only  linked  by  the  data  file  (disc  or  tape  )  produced  by  the  preprocessor  and 
utilized  in  ION1.  PROGRAM  TABGEN  requires  the  following  SUBROUTINES 
SICOJT,  DKSICO,  MAGFIN,  GK,  DKGK,  and  the  library  functions  AMOD, 
ATAN,  COS,  SIN,  SQRT.  All  internal  data  transfer  occurs  through  the  calling 
sequences;  the  external  data  transfer  consists  of  input  coming  from  the  data 
card  deck  and  the  ionospheric  coefficient  tape  and  of  output  of  f0F2-hi  tables 
to  disc  or  tape,  all  in  PROGRAM  TABGEN.  PROGRAM  ION1  requires  the 
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following  SUBROUTINES  REFRC.l,  PROFL2,  BETA  and  the  library  functions 

ABS,  AMOD,  ATAN,  COS,  FLOAT,  SIN,  SQRT.  Ail  internal  data  transfer 
occurs  through  the  labeled  common  blocks  and  through  the  calling  sequences. 
The  external  data  transfer  consists  of  input  coming  from  the  data  card  deck 
into  ION1,  from  the  preprocessed  disc  or  tape  file  with  f0F2-b,  tables  into 
SUBROUTINE  REFRC1,  and  of  output  of  the  results  from  PROGRAM  ION1  to 
the  line  printer.  The  second  attached  diagram  shows  the  program  structures, 
the  data  files  are  described  in  Section  3,  3.  1,  and  the  functions  performed  by 
the  preprocessor  and  reduction  program  are  outlined  in  Section  3.  4. 

Whenever  ionospheric  predictions  are  desired,  PROGRAM  ION  should 
have  pref«  rence  over  the  program  set  TABGEN-IONI,  ION  will  yield  more 
accurate  results  than  ION  l  where  approximations  are  introduced  through  inter, 
polating  the  ^FZ-h,  tables  and  through  bypassing  the  iteration  on  the  height 
estimate  nf  the  ionosphere.  ION  also  has  the  additional  features  not  included 
in  ION  1  of  computing  range  rate  corrections  for  range  differencing,  of  plotting 
the  ionospheric  profile,  and  of  updating  the  predictions  with  actual  ionospheric 
observations.  For  many  applications  ION  will  be  suited  even  for  real-time 
processing,  The  program  set  TABGEN-IONI  should  only  be  used  when 
stringent  core  space  and/or  run  time  requirements  are  imposed  that  cannot 
be  met  by  PROGRAM  ION,  or  when  program  modifications  for  special 
applications  are  attempted.  Running  PROGRAM  TABGEN  in  a  preprocessing 
mode  results  in  the  significant  core  space  and  run  time  reduction  of  PROGRAM 


Thi*  paragraph  contain*  the  detailed  technical  description*  of  the 
computer  program  component*  identified  in  Paragraph  3,  1  of  thi*  specificati' 
The  instruction  Hating*  contained  in  Appendix  I  *pecify  the  exact  configuraH.  a 
of  the  Bent  Ionospheric  Program  ION  *  .id  the  alternate  version  TABGEN  -  ION  l; 


Following  are  specifically  the  descriptions  fort 


CPC  No.  1  -  PROGRAM  ION 

CPC  No.  2  -  SUBROUTINE  REFRAC 

CPC  No.  3  >  .  SUBROUTINE  PLOTNH 

CPC  No.  4  .  SUBROUTINE  PROFL1 

CPC  No.  5  .  SUBROUTINE  PROFL2 

CPC  No.  6  -  SUBROUTINE  BETA 

CPC  No.  7  -  SUBROUTINE  SICOJT 

CPC  No.  8  -  SUBROUTINE  DKSICO 

CPC  No.  9  -  SUBROUTINE  MAGFIN 

CPC  No.  10  -  SUBROUTINE  GK 

CPC  No.  11  -  SUBROUTINE  DKGK 


Particular  to  all  subroutines  Is  the  fact  that  none  of  the  input  variables 
transferred  through  common  or  the  calling  sequence*  are  modified  during 
excution  of  the  program  code.  The  units  internal  to  all  subroutines  are 
kept  in  meters  for  distances,  radians  for  angle*  and  times,  meters /second 
for  linear  velocities,  radians /second  for  angular  rates,  MHz  for  frequencies 
and  Gauss  for  magnetic  field  strength. 


Included  are  also  the  description*  of  the  routines  that  are  required 
in  addition  ot  the  ones  listed  above  for  the  alternate  version  of  the  ionospheric 
program,  consisting  of  separate  preprocessor  and  reduction  programs: 

CPC  No,  12  -  PROGRAM  TABGEN 

CPC  No.  13  -  PROGRAM  10N1 

CPC  No.  14  -  SUBROUTINE  REFRC1 


9 


3,  2.  1  Computer  Program  Component  1 


CPC  No<  1,  main  PROGRAM  ION,  is  written  in  FORTRAN  code.  It 
handlea  the  card  input  and  the  printing  of  the  reaulta  for  the  entire  program, 
except  for  the  Uat  and  plot  of  the  profile  which  is  done  by  SUBROUTINE  PLOTNH 
upon  call  from  ION.  ION  transfer*  the  input  conditions  through  commons  /'EVAL/ 
and/UPDT/  and  by  calling  SUBROUTINE  REFRAC  receives  the  computed  profile 
parameters  and  refraction  corrections  through  common /CORR/. 

‘  ’  j 

3,2.  1.1  CPC  No,  1  Description 

ION  reads  the  selections  for  the  output  and  update  options  from  cards,  i 

it  reads  the  station,  satellite  and  time  information  for  the  condition  to  be  j 

evaluated,  and  as  needed,  reads  the  solar  data  from  cards.  If  the  option 
for  updating  the  predictions  with  measured  ionospheric  data  was  chosen,  the 
number  of  observations  to  be  used  for  tho  update  and  the  corresponding  observation 
along  with  station  and  time  information  are  read  from  cards.  Up  to  eight 
measurements  can  be  used  simultaneously  for  updating  any  one  evaluation 
condition,  All  input  data  is  listed  for  reference  in  the  print  out. 

The  input  data  is  converted  to  the  internal  units  of  meters  for  distances 
and  radiar  s  for  angles  and  times.  The  variables  specifying  the  evaluation  j 

condition  are  transferred  through  common/EVAL/,  the  update  conditions  through  : 
l-.ommon/UPDT/to  SUBROUTINE  REFRAC.  Through  REFRAC  and  other  routines  t 

i 

called  by  REFRAC  ionospheric  profile  parameters,  vertical  and  angular  electron  ; 
content,  refraction  corrections  to  elevation  angle,  range,  and  instantaneous  i 

range  rate  are  computed  as  desired  and  returned  to  ION  through  common/CORR/,  i 
ION  print!  the  results  as  requested  and  calls  SUBROUTINE  PLOTNH  for  an 
electron  density  profile  plot  and  list,  whenthis  type  of  output  is  specified. 

i 

lithe  refraction  correction  to  range  rate,  obtained  by  range  differencing 
over  a  fintc  time  during  which  the  ionosphere  can  undergo  changes,  Is  requested, 
tho  input  lor  the  evaluation  condition  obove  relates  to  the  first  range  observation, 
and  additisnal  satellite  and  time  information  that  is  read  from  card  relates  to  .he 
last  ranga  observation  used  in  the  differencing  technique.  Upon  rcturr  the 

J 

/ 
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second  range  correction  from  F.EFRAC,  ION  computes  the  request!  .  r 

correction  by  differencing  the  two  range  correction*  and  dividing  by  th 
interval;  the  result  is  printed, 

Any  number  of  evaluation  condition*  can  be  proce**ed  by  supplying 
additional  input  data  and  repeating  the  program  *tep*  outlined  above,  For 
more  detail*  about  the  input  data  and  output  option*  refer  to  the  input  data 
description  under  3,  3,  1. 

3,2,  1,2  CPC  No,  1  Flowchart 

The  flowchart  is  shown  on  the  page  following  3.  2,  1.  5, 

3.  2,  1,  3  CPC  No,  1  Interface* 

a)  Library  *ubprogram<  required!  none 

b)  Other  *ubprogram*  called;  SUBROUTINES  REFRAC,  PLOTNH 

c)  Calling  program;  none 

d)  Calling  sequence*  PROGRAM  ION 

e)  Common  blocks*  EVAL,  UPDT,  CORR 
Variables  in  common; 

See  description  for  EVAL,  UPDT,  CORR  under  SUBROUTINE  REFRAC, 
CPC  No,  2. 

f)  File  requirements:  card  reader,  line  printer 

The  requirements  for  the  input  datacard  file  are  specified  under  3,  3.  1, 


3,  2.  1,  4  CPC  No,  1  Data  Organisation 
Variables  defined  in  data  statements; 


Name 

Dimension 

Description 

LYRMO 

1 

=  0,  initialization  constant  for  (year  >  1  QO+month) 

IDRD 

1 

=  0,  default  condition;  range  rate  correction  for 
observation  over  finite  time  is  not  desired 

I  OPT 

1 

=  1,  default  condition;  computation  of  critical 
frequency  and  corresponding  height  is  desired 
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Name 

MEAS 


Dimension  Description 


Array  containing  hollerlth  data  for  print  out 


Other  constant*  defined  in  data  statements 

QO=0,  01000=1000,  Q3600=36  00;  DR=1",  HR=1  houir,  PI2=360*  converted  to  radians, 
Important  variable*  are  described  under  3.  2,  1,  3  «)'Of  SUBROUTINE  REFRAC, 

CPC  No.  2. 


3,2,  1,  5  CPC  No,  1  Limitation* 

Up  to  e.ght  measurement  entrie*  can  be  u*ed  simultaneously  for  updating 
the  prediction*  for  any  one  evaluation  condition.  If  update  with  more  than  eight 
condition*  i*  requested,  the  program  u*e*  the  first  eight  entrie*,  ignoru*  any 
additional  input  data  and  print*  a  message  to  that  effect. 

Error  test*  on  the  sequence,  unit*  and  format*  of  the  input  data  are  not 
performed,  except  on  the  date*  of  the  solar  data  card*,  However,  mistake* 
in  the  set  up  oi  the  card  deck  are  revealed  in  the  printout  of  the  input  data 
that  is  listed  along  with  the  result*. 
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PROGRAM  ION  (continued) 


range  ameren ^ 
sCing  wanted,-^  V _ y 


VRe*d  &£  print 
\  additional 
\  condition 


/Call  REFRa"^ 
/for  AR  for  new  \ 
\  condition  / 


Compute  aA 
for  range 

diffe  rencing 


L 


Print  AR 


3,2,  1  Computer  Program  Component  2 


CPC  No.  2,  SUBROUTINE  REFRAC,  i*  written  in  FORTRAN  code  and 
is  called  from  the  main  PROGRAM  ION,  REFRAC  prepare!  the  coefficient 
and  aolar  input  data,  it  obtain*  the  ionospheric  profile  parameters  via  PF.OFL1 
and  PROFL2,  it  performs  an  optional  update  using  up  to  eight  observation 
entries,  it  computes  the  ionospheric  refraction  corrections  AR  for  range, 

A  A  for  instantaneous  range  rate  and  obtains  the  refraction  correction  AE 
for  the  elevation  angle  via  LZTA, 


3,  2,  1,  1  CPC  No,  2  Description 

REFRAC  prepares  the  coefficients  to  be  used  in  SUBROUTINE  DKSICO  for 
the  computation  of  the  time  dependent  coefficients  which  in  turn  are  required 
for  the  computation  of  critical  frequency  f0F2  and  M(3Q00)F2,  At  first  it  is 
checked  if  the  coefficients  are  already  available  for  the  desired  date,  and  if  not 
available,  the  proper  coefficients  are  read  from  tape.  These  general  coefficients 
are  valid  for  any  condition  and  do  not  have  to  be  updated  or  replaced,  but  can 
be  adjusted  for  any  time  in  the  past  or  future, 


The  general  f0F2  coefficients  were  derived  using  the  work  of  Jones  and 
Obitts  (Reference  6);  they  provide  annual  continuity  and  are  valid  for  approx¬ 
imate  10  day  periods,  for  the  spans  from  day  1  to  10,  day  11  to  20,  and  day  21 
to  30  (or  28,  29,  3  1  )  of  each  month,  There  are  coefficients  for  36  periods  to 
cover  the  whole  year.  The  general  f0F2  coefficients  W.,  ,  ,  „  represent  the 
coefficients  to  a  second  order  polynomial  in  the  12-month  running  average  of 
solar  fluxF_a  (observed  Ottowa  10.  7  cm  solar  flux).  They  are  evaluated  for 
the  specific  Fia  of  the  evaluation  date  to  yield  the  specific  f0F2  coefficient  set 
U;  ,  „  (stored  in  array  U)  used  in  SUBROUTINE  DKSICO; 


^  1 1  >t  =  Wx  .  i .  K  +  Wg  , ,  , ,,  «  F19  +  W 3  ,  i  , 

k=  0,  1, _ 75. 


,  for  i  =0, 1  ,  ...  12  and 
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The  general  M(3QQ0)F2  coefficiente  available  from  NOAA,  Boulder, 
are  valid  for  monthly  periods,  There  are  coefficents  for  12  periods  to  cover 
the  whole  year,  and  for  each  period  there  are  two  sets  V,  ,  k  (0)  and  V[  ,  „  (100), 
one  for  a  12-month  running  average  of  sunspot  number  Sia=  Q  and  the  other 
for  =  100.  The  coefficients  are  adjusted  by  interpolating  or  extrapolating 
the  two  sets  to  the  specific  Sj a  of  the  evaluation  date  yielding  the  specific 
M(3000)F2  coefficient  set  U,  ,k  (stored  in  array  UM)uaed  in  SUBROUTINE  DKSICO; 


U 


l  >  it 


V 


l  >  k 


;o)  + 


s. 


is 


.V,  ,  k  (100)  -  Vj  ,  „  (0)  j  *  —  ,  for  i  =0,  1,  ...  8  and 
k=  0,  1,  ...  48. 


The  10,  7  cm  Ottowa  solar  flux  data  is  prepared  for  use  in  SUBROUTINES 
PROFL1  and  PROFL2.  The  difference  AF  between  the  daily  value  F  and  the 
12-month  running  average  of  the  solar  flux  is  formed,  AF  =  F-F-.  s.  If  the  daily 
solar  flux  in  not  available,  F;  a  is  substituted.  If  the  daily  solar  flux  is  greater 
than  130,  liO  is  substituted  which  is  a  limit  imposed  by  the  data  base  on  which 
development  of  the  model  was  founded. 


The  first  parameters  for  the  ionospheric  profile,  the  critical  frequency 
f0F2  and  th<s  corresponding  height  h,  are  obtained  via  SUBROUTINE  PROFL1, 


On  option  REFRAC  updates  the  predicted  fnF2  with  observations  of  fQF2 
or  with  ver.ical  or  angular  electron  content  reduced  from  Faraday  rotation 
measurements  from  other  stations.  Up  to  eight  update  observations  of  cither 
type  separated  by  different  amounts  in  time  ana  space  from  the  evaluation  time 
and  station  can  be  accepted.  To  obtain  the  best  possible  update,  the  observation 
times  and  stations  should  be  the  closest  to  the  evaluation  condition  available, 

In  any  case,  the  update  station  should  be  within  2000  km  of  the  evaluation 
site. 


If  the  c  bse  r  vat  ion  is  angula  r  ele  ctr  on  c  ontent  Nr  K ,  it  is  reduced  to  total 
vertical  electron  content  NT  by, 


Nr  =  NT» 


R, coiE ^  a 
R,  +  h,' 
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E  being  the  eievation  angle  of  the  observation,  and  R,  the  mean  earth  r.  diuL, 

For  each  update  observation  the  predicted  fQF2  is  obtained  by  calling  SUBROU1  (NL 

PROFJ-1,  and  the  update  ratio  r  is  formed  for  fQF2  observations , 

foF2  oba. 

r  =  — - 

f0F2  pred. 

If  the  obaervation  is  electron  content,  the  additional  profile  parameter  NT  /N, 

1b  obtained  via  SUBROUTINE  PROFL2,  and  the  following  ratio  1b  formed, 


1.24  *  1 01 0  f0F2ipred. 


where  f0F2  ia  in  MHz 


and  since  the  maximum  electron  density  ia  N,=  1.24*  1  013  f0F2a  ,  and  NT  ia 
approximately  proportional  to  f0F22,  the  electron  content  information  ia 
reduced  to  a  f0F2  ratio, 

/  N,  obi.  foF2  obB. 


/  Nr  obB.  foF2  obB. 

Nt  pred.  f0F2  pred. 

If  there  iB  only  one  update  condition,  the  ratio  r  ia  used  for  the  final 
ratio  R  to  update  f0F2.  If  several  n  conditions  are  uaed  for  the  update, 
a  weighted  mean  technique  combines  all  n  ratios  rt  to  the  final  ratio  R 
having  aa  weights  w,  the  time  differences  At,  between  observation  and 
evaluation  timea  and/or  the  earth  central  angles  a,  between  the  ionospheric 
points  at  which  the  rays  from  the  observation  and  evaluation  stations  pass 
through  the  ionosphere, 


R  =  1  =  1 


l  i- 

U  w, 
i  =  i _ 

i  -l 

Ld  w, 


w,  =  At.  , if  observations  are  from  one  station  at  different  times, 
w,  =  a,  ,  if  obse  rvations  are  from  several  stations  at  the  same  time, 
w,  =  At.  a,,  if  observations  are  from  several  stations  at  different  times. 

At  =  It -to  I  and 

cos  a  =  sinb  sinb0  +  cos<£  cos  b0  cos  ( X- X0 ), 
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where  t,  <p,  \  and  to  ,  <^0  ,  \0  are  the  time,  latitude  and  longitude  of  the  ior.o- 

•pheric  point*  for  evaluation  and  ob»ervation  condition  respectively,  The 
final  ratio  R  up-data*  the  critical  frequency  by  the  tame  overall  percentage 
by  which  the  predictions  deviate  from  the  ionospheric  observations, 

£qF2  upd,  =  fQF2pred.«  R. 

By  calling  SUBROUTINE  PROFL2  the  remaining  profile  parameters 
are  obtained:  y,  the  half  thickness  nf  the  bottomside  bi-parabolic  layer, 
yt  the  half  thickness  of  the  topside  parabolic  layer,  lq,  Iq,  kj  the  decay 
constants  for  the  lower,  middle,  and  upper  section  of  the  topside  exponential 
layer,  Nr/N,  the  ratio  of  the  total  integrated  electron  content  to  the  maxi¬ 
mum  electron  density,  rr.  the  multiplier  of  the  ,  rate  of  change  in  height, 
term  in  the  range  rate  equation. 

The  one-way  ionospheric  refraction  correction  AE  to  the  elevation 
angle  E  is  calculated  via  SUBROUTINE  BETA.  The  total  integrated  electron 
content  N7  a.  ons  a  vertical  path  through  the  ionosphere  and  the  angular 
content  alonp  the  line  of  sight  NTAare  computed  as; 

_ n. _ 

NT  =1.1)4  *  10lofQF2a  (Hil— ),  N(,  =  J  1  .  (  R«  cobE  f 

N«  '  R,+h,  ' 


The  one-way  ionospheric  refraction  correction  to  range  A&  is  given  by 
the  equation 


n  40.  3  X  !CTiaN,  4Qt  3  >  1,  24  «  10~a 

4  '  e  n  (R-c°,E)a  =  /  (R. 

AR,+h,  '  V  " '  R,  +h,  ' 


where  f  is  the  transmission  frequency, 
are  uplink  and  downlink  frequencies. 


(-pr  +  jr)  fu  *nd  £d 

fU  id 


The  one-way  ionospheric  refraction  correction  to  range  rate  AR  consists 
of  two  termi,  one  multiplied  by  the  altitude  rate  h,  the  other  by  the  elevation 
rate  t,; 


IB 


This  range  rate  correction  foxmulation  applies  only  to  instantaneous 
range  rate  measurement*,  since  it  assumes  that  the  only  variation  in  the 
total  electron  content  over  the  time  of  the  observation  is  due  to  the 
positional  change  of  the  satellite  and  that  the  ionosphere  between  station 
and  satellite  remains  constant  for  the  duration  of  the  measurement. 
Corrections  to  range  differencing  are  discussed  under  3,  2,  1,  5, 

The  signs  of  the  refraction  corrections  are  set  for  the  corrections 
to  be  subtracted  from  their  respective  observations,  The  units  in  all 
equations  above  are  kept  in  meters,  met  ' /second,  radians,  radians/ 
second  and  MHz, 

3.2,  1,2  CPC  No.  2  Flowchart 

The  flowchart  is  shown  on  the  page  following  3.2.  1.  5, 

3.2,  1,3  CPC  No,  2  Interfaces 

a)  Library  subprograms  required:  ABS,  ATAN,  COS,  SIN,  SQRT 

b)  Other  subprograms  tailed*  SUBROUTINES  PROFL1,  PROFL2,  BETA 

c)  Calling  programs*  PROGRAM  ION 

d)  Calling  sequence;  SUBROUTINE  REFRAC 

e)  Common  blocks;  EVAL,  UPDT,  CORR 
Variables  in  common; 


Common 

Name 

Variable 

Name 

Dimension 

I/O  Description 

EVAL 

FS 

1 

I  Transmission  frequency  (MHz  ) 

EVAL 

FLAT 

1 

I  Latitude  of  station  (radians) 

EVAL 

FLON 

1 

I  Longitude  of  station  (radians) 

EVAL 

ELEV 

1 

I  Elevation  to  satellite  (radians) 
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Common  Variable  Dimension  I/O  Description 

Name  Name 


EVAL 

AZ 

1 

I 

Azimuth  to  satellite  (radians) 

EVAL 

HS 

1 

I 

Height  of  satellite  (m) 

EVAL 

EDOT 

1 

I 

Ele  vation  rate  (radians /sec  ) 

EVAL 

HDOT 

1 

I 

Altitude  rate  (m/sec) 

EVAL 

TIME 

1 

1 

Universal  time  (radians ) 

EVAL 

FLXD 

1 

1 

Daily  solar  flux 

EVAL 

SIS 

1 

I 

12-month  running  average  of  sun., 
spot  number 

EVAL 

SIF 

1 

I 

12-month  running  average  of  solar  flux 

EVAL 

IYR 

1 

I 

Year  (last  2  digits  ) 

EVAL 

MON 

1 

I 

Month  (  =  1  through  12) 

EVAL 

IDAY 

1 

I 

Day  ( =  1  through  3 1  ) 

EVAL 

IOPT 

1 

I 

Control  constant  for  optional  computa- 

tions;  =1  to  compute  f0F2  and  h,  , 

-2  to  also  compute  remaining  profile 
paiAme'ters  and  electron  content, 

=  3  to  compute  AH  in  addition,  -4  to 
also  compute  AR 


EVAL 

IDEL 

1 

I  Control  constant  to  compute  AE  be¬ 
sides  profile  parameters  and  electon 
content,  -0  compute,  =1  not 

EVAL 

IDRD 

l 

I  Flag  to  eliminate  unnecessary  comput 

tions  during  calculation  of  the  second 
range  correction  used  in  the  different 
ing  for  the  range  rate  correction, 

=  0  for  first,  =1  for  second  calculation 

EVAL 

IUPDT 

1 

I  Update  flag,  -0  no  update,  -1  update 

EVAL 

IT  P 

1 

I  Unit  assignment  of  general  iono¬ 

spheric  coefficient  tape 

UPDT 

ULAT 

8 

I  Latitudes  of  update  stations  (radians) 

UPDT 

ULON 

8 

I  Longitudes  of  update  stations  (radians 

ULEV 


Elevation  angles  of  observations  (radians) 


UPDT 


8 
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Common  Variable  Dimension  I/O  Description 


Name 

Name 

UPDT 

UZIM 

8 

I 

Azimuth  angles  of  observations  (radianB) 

UPDT 

UT 

13 

I 

Universal  time  of  observations  (radians) 

UPDT 

CBS 

8 

1 

Observation  of  f0F2,  vertical  or  angular 
electron  content  (MHz  or  electrons/ma  ) 

UPDT 

ITYP 

8 

1 

Observation  type,  si  for  f0F2,  =2  for 
vertical,  =3  for  angular  electron  content 

UPDT 

NUPDT 

1 

I 

Number  of  update  conditions 

CORR 

DRANG 

1 

0 

Range  correction  (m) 

CORR 

DRATE 

1 

O 

Range  rate  corroction  (m/sec) 

CORR 

DELEV 

1 

0 

Elevation  angle  correction  (radians) 

CORR 

F0F2 

1 

0 

Critical  frequency  (MHz) 

CORR 

HM 

1 

o 

Height  at  maximum  electron  density 
(mete  rs ) 

CORR 

YM 

1 

o 

Half  thickness  of  the  bottomsida  hi- 
parabolic  layer  (meters) 

CORR 

YT 

1 

0 

Half  thickness  of  the  topside  parabolic 
layer  (meters  ) 

CORR 

XK 

3 

o 

Decay  constants  of  lower,  middle  and 
upper  section  of  the  exponential  top¬ 
side  layer  (1 /meter) 

CORR 

TOTN 

1 

o 

Total  vertical  electron  content 
(e /m*  column ) 

CORR 

TOTNA 

1 

o 

Total  angular  electron  content 
(e/ma  column) 

f)  File  requirements*  general  coefficient  input  tape,  line  printer 

The  format  of  the  general  coefficient  tape  is  described  under  3.  3,  1, 
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3,2.  1.4  CPC  No.  2  Data  Organization 

Variables  defined  in  data  statements; 


Name 

Dimension 

Description 

R 

1 

Mean  earth  radius  (meters) 

STRS 

1 

Approximate  height  of*tationary  satellite  used 
when  updating  with  observed  electron  content 
(meters ) 

TOL 

1 

Tolerance  for  difference*  in  positions  or  obser. 
vation  times  of  multiple  update  stations  below 
which  they  are  assumed  to  be  identical  (radians 

MONDY 

1 

Initialization  constants  for  last  and  first 

MOND 

1 

(month x  100  +  day)  for  which  coefficients  are 
in  core 

LYRMO 

1 

Initialization  constant  for  (year  *  100  +  month) 

Other  constant  a  defined  in  data  statements; 

QO=0,  Ql  =  l,  Q 1 00=  1  00,  Q130=130,  QP1  =  .  1,  QNM=1,  24  *  1  O10,  RN3  =  . 49972; 
PI=180',  P[2=360*  converted  to  radians. 

Other  important  variables  are  described  under  3.2,  1.3  e). 

3,  2.  1.  5  CPC  No,  2  Limitations 

The  daily  value  of  solar  flux  transferred  to  SUBROUTINE  PROFL2 
for  the  computation  of  the  decay  constants  for  the  topside  exponential 
profile  is  t.-uncated  at  a  maximum  value  of  130.  This  is  the  boundary 
that  was  imposed  by  the  data  base  during  the  model  development  and  is 
thus  a  limit  to  the  model  since  the  extension  of  solar  flux  beyond  130  could 
result  in  in/alid  profiles. 

The  dimensions  of  several  arrays  restrict  the  update  procedure  to  be 
applied  to  tie  predictions  of  any  one  evaluation  condition,  to  not  include 
more  than  light  observation  entries. 

The  ra  ige  rate  correction  formula  in  this  routine  applies  only  to 
instantaneo  is  range  rate  measurements,  since  it  is  assumed  that  the  only 
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variation  in  electron  content  over  the  time  of  the  obaervation  ia  due  to  the 
positional  change  of  the  satellite  and  that  the  ionosphere  between  station 
and  satellite  remains  constant.  If  the  range  rate  corrections  are  desired 
for  observations  obtained  by  range  differencing  over  a  finite  time  interval 
during  which  the  ionosphere  can  undergo  significant  changes,  a  range 
correction  differencing  technique  should  be  used  over  the  same  time 
interval.  This  type  of  correction  can  optionally  be  requested,  it  requires 
additional  satellite  and  time  information  and  is  handled  directly  in  PROGRAM 
ION,  CPC  No.  1. 

If  the  ionospheric  coefficients  are  not  found  on  the  tape  for  the  evaluation 
date,  an  error  condition  has  occurred,  -  ...e • faage  is  printed  out,  and 
control  is  transferred  to  PROGRAM  ION  to  proceed  with  the  next  data  case, 
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CPC  No.  2  Flowchart,  SUBROUTINE  REFRAC 
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i,'  Jr>.OUlTNE  REFRAC  (contiruridj 


3.  2.  1  Computer  Program  Component  3 


CPC  No. 3,  SUBROUTINE  PLOTNH,  ia  written  in  FORTRAN  code.  It 
ia  called  from  the  main  PROGRAM  ION  and  liata  and  plot*  the  electro- 
denaity  veritua  height  profile, 

3,  2.  1.  1  C IKC  No.  3  Deacription 


PLOTM-I  plota  a  graph  of  electron  denaity  N  veraua  height  h  at  2  5  km 
height  inert menta  from  25  km  to  1000  km  and  it  printa  a  list  of  electron 
denaitiea  for  corroaponding  height  valuea  from  25  km  to  2000  km  at  25  km 
incrementa 

The  elt  ctron  denaity  ia  modeled  differently  in  five  height  layera  (aee 
Figure  2  in  section  6.1).  k^  ,  kg  ,  kg  denote  the  decay  conatanta  for  the  lower, 
middle  and  upper  section  of  the  exponential  topside  profile,  and  yt ,  y,  are 
the  valuea  of  half  thickness  for  the  topside  parabolic  layer  and  for  the 
bottomaide  bi. parabolic  layer  reapectively.  The  height  limits  for  each 
layer  are  firat  determined  and  the  value  of  electron  density  at  the  start 
point  of  the  various  layers  N, ,  N0,  N:  ,  Na  .  The  height  increments  measured 
from  the  atirt  point  of  the  various  layers  are  denoted  an  variables  b^by, 
ai  i  Rg  ,  a3  .  The  electron  density  equations  are; 


N  =  N,  1 


h.  -y. 


N  =.N0  e 


N  =  Nx  e 


N  =  No  e 


'  **9  *3 


h,  s.  h  ^  hoaht+d 


^  *  h  <  ht  =ho+(  1  01  2km-h0  )/3 


hj  jc  h  <  h^^hj^+f  1 01 2  km-ho  )/3 


hg  *  h  *  2000  km 
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where  h,  is  the  heighi  at  the  maximum  electron  density,  d  is  the  distance 

above  h#  at  which  the  lower  exponential  layer  start*,  and  the  electron  densitie 
at  the  start  points  of  the  various  layers, 

N,  =  1.  24  *  1010  i0TZa 

No  =  N.  (l  .  il ) 

'  y*- 

Nx  =  N0  e’ki(hl  -h°  ) 

Na  =  N\  e_ka(ha  } 

3,  l,  1.2  CPC  No,  3  Flowchart 

The  flowchart  as  shown  on  the  page  following  3,  Z,  1,  5, 

3.  2.  1.  3  CPC  No.  3  Interfaces 

a)  Library  subprograms  required:  EXP,  LOGIO.SQRT 

b)  Other  subprograms  called:  none 

c)  Calling  program;  SUBROUTINE  REFRAC 

d)  Calling  sequence:  CALL  PLOTNH  (F0F2,  HM,  YM,  YT,  XK) 

Variables  in  calling  sequence: 

Name  Dimension  I/O  Description 

F0F2  1  1  Critical  frequency  (MH* ) 

HM  1  1  Height  at  the  critical  frequency  (meters) 

YM  1  I  Half  thickness  of  the  bottom  bi-parabolic  layer 

(meters ) 

YT  1  I  Half  thickness  of  the  topside  parabolic  layer 

(mete  rs  ) 

XK  3  I  Decay  constants  for  lower,  middle,  and  upper 

section  of  the  topside  exponentia  1  laye r  ( 1  / mete r ) 

e)  Common  blocks:  none 

f)  File  requirements;  line  printer 


3.2.  1  4  CPC  No.  3  Data  Organization 


Variables  defined  in  data  statements 
Name  Dimension  Deacription 

IBLANK  1  Hollorith  "blank"  symbol  uaed  for  plotting 

MARK  1  Hollorith  "  *  "  lymb  )1  used  for  plotting 

Other  constants  listed  in  data  statement; 

Q0=0,  Qlsl,  Q3  =  3,  Q124E  =  1.  24x1c10,  Q 1 01 2E  =  1 012000,  Q1025E  =  1025000, 
Q25E=25000,  Q10*10,  Q27  =  27,  Q2025E=2025000, 

Other  important  variables  ere  described  under  3,  2,  1,  3  d), 

3,  2.  1.  5  CPC  No,  3  Limitations 

If  electron  density  values  are  computed  smaller  than  1010  or  larger 
than  5  x  li)13  (electrons /meter3  ),  they  exceed  the  limits  of  the  graph  and 
automaticnlly  are  not  plotted.  Since  these  cases  do  not  normally  involve 
error  conditions,  a  message  is  not  required  and  the  values  are  printed 
as  computed  in  the  electron  density  versus  height  list, 
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CPC  No.  j  Flowchart,  SUBROUTINE  PLOTNH 
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3.2.  1  Computer  Program  Component  4 


CPC  No.  4,  SUBROUTINE  PROFL1,  ia  written  in  FORTRAN  code.  It 
la  called  from  SUBROUTINE  REFRAC  and  computes  the  ionospheric  profile 
parameters  critical  frequency  f0F2  and  the  corresponding  height  h,  at  the 
location  where  the  wave  passes  through  the  ionosphere. 


3.  2.  1.1  CFC  No.  4  Description 

PROFL1  computes  the  ionospheric  characteristics  fQF2  and  M(3000)F2 
following  the  analysis  of  Jones,  Graham  and  Leftin  (Reference  5).  First 
the  trigonometric  functions  of  the  multiples  of  the  Greenwich  hour  angle  t, 
-180,fit*  180°,  t  =  0  at  Greenwich  noon,  are  computed  via  SUBROUTINE 
SICOJT  for  use  in  DKS1CO.  The  time  dependent  coefficients  are  computed 
via  SUBROUTINE  DKSICO  based  on  the  coefficient  sets  U^  *  prepared  in 
REFRAC.  Utilizing  the  f0F2  and  M(3000)F2  coefficient  sets  (in  arrays  U 
and  UM )  the  time  dependent  coefficients  respectively  for  the  fQF2  and 
M(300C)F2  evaluation  are  prepared, 


Defined  by  the  latitude  d  and  longitude  \  at  which  the  ray  from  station  to 
satellite  passes  through  the  ionosphere  is  the  ionospheric  point.  It  iB 
calculated  as  a  function  of  the  station  latitude  b>s,  longitude  and  the 
elevation  an<le  E  and  azimuth  angle  A  to  the  satellite; 


</> 

\ 


arc  sin  ^sind>5 
\s  +  arc  sin  ( 


cob  a  j-  cos  bj 

sin  A  sin  a  '' 
cosii 


sin  a  cos  A 


\ 

J 


whe  re 
point, 


is  tie  earth  central  angle  between  the  station  and  the  ionospheric. 


E  -  arc  sin 


R, cosE 
R,  +  h. 


)■ 


R,  is  the  mi  an  earth  radius,  and  h,  is  the  height  of  the  ionosphere  at  the 
maximum  electron  density  above  the  surface  of  the  earth.  Since  h,  i.  to 
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be  determined  later  or  in  this  subroutine,  a  first  estimate  of  ht  is  required 
and  assumed  as  h„=300  km.  After  computing  the  actual  hB  prediction,  the 
new  value  is  compared  with  the  estimate  and  if  it  deviates  by  more  than  1  km, 
all  computations  starting  with  the  determination  of  the  ionospheric  point 
are  repeated  using  the  new  hB  . 

The  position  dependent  functions  requi:  ed  for  the  fQF2  and  M(300Q)F2 
computations  are  all  evaluated  at  the  ionospheric  point  which  can  differ  by 
up  to  21*  from  the  station  position.  First  the  earth's  magnetic  field 
components  X-north,  Y-eaet  and  Z-vertical  up  are  computed  at  the  ionospheric 
point  via  SUBROUTINE  MAGFIN,  and  they  form  in  turn  the  modified  magnetic 
dipx  as  a  function  of  the  magnetic  dip  I; 

I  -Z 

x  =  arc  sin  ,  ,  I  =  arc  tan  — 

V  Ia+cosd  Vxa+Ya 

Eased  on  the  following  coordinates,  ionospheric  latitude,  longitude  and 
modified  magnetic  dip,  SUBROUTINE  GK  evaluates  the  geographic  coordinate 
functions  for  the  fnF2  computation.  Extracted  from  these  functions  is  the 
subnet  which  forms  the  geographic  coordinate  functions  needed  for  the 
M(3000)F2  computation. 

SUBROUTINE  DKGK  multiplies  and  sums  the  proper  sets  of  time 
dependent  coefficients  and  position  dependent  functions  and  forms  M(3000)F2. 
With  the  Appleton-Beyr.on  equations  (Reference  1  ),  a  second  order  poly¬ 
nomial  in  M(3000)F2,  the  height  of  the  maximum  electron  density  is  obtained 
in  mete  r  a ; 

h,  =  1  1  346 .  92  -  526.  40  *  M(3000}F2  +  59.  825  [M(3000)F2|2  }  «  1 03 

h„  is  compared  with  its  estimate  and  if  the  difference  is  greater  than  1  km, 
the  computations  above  starting  with  the  ionospheric  point  determination 
are  iterated  on  using  the  new  value  for  hB  , 
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Using  the  proper  time  dependent  coefficient*  and  position  dependent 
functions,  SUBROUTINE  DKGK  computes  the  10  day  mean  of  the  critical 
frequency  which  then  is  adjusted  for  day  to  day  changes  in  the  ionosphere 
and  for  additional  magnetic  latitude  variations,  following  the  model 
description  in  Section  6,  1.  The  magnetic  latitude  of  the  ionospheric  point 
is  determined  as, 

dv  =  arc  sin  [sind  sin^  +  cosd  cos  dp  cos(A-\p)], 

where  dF ,  \p  ire  the  latitude  and  longitude  of  the  magnetic  north  pole  and 
Interpolating  the  model  conitants  (array  CENT)  tod,  results  in  ca, 

The  daily  vaiiation  from  the  mean  value  is  dependent  on  t\ F,  the  difference 
between  the  daily  value  and  the  12 -month  running  average  of  the  solar  flux 
and  on  the  model  constant  cx  (variable  PER).  The  f0F2  computed  by  DKGK 
is  multiplied  bv  the  adjustment  factor  (cjAF+Cj  )  to  yield  the  final  predicted 

f0F2. 

The  uniti  in  the  above  equations  are  kept  in  meters,  radians  and  MHz. 

3,  2,  1,  2  CPC  No,  4  Flowchart 

The  flowchart  i*  shown  on  the  page  following  3.  2.  1,  5, 

3.  2.  1,  3  CPC  No.  4  Interfaces 

a)  Library  subpr  grams  required:  ABS,  ATAN,  COS,  SIN.SQRT 

b)  Oth.»r  Bubprograms  '.ailed:  SUBROUTINES  SICOJT,  DKSICO,  GK, 
MAGFIN,  DKGK 

c)  Calling  program-  SUBROUTINE  REFRAC 

d)  Calling  uequence:  CALL  PR OFLl(FLAT,  FLGN,  ELEV,  AZ,  TIME, 

DF LUX,  U,UM,  OLAT,  OLON,  F0F2,  HM,  HLAT  ) 
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Variables  in  calling  sequence: 


Name 

Dimension  I/O 

Description 

FLAT 

1 

I 

Station  latitude  (radians) 

FLON 

1 

I 

Station  longitude  (radians) 

ELEV 

1 

I 

Elevation  angle  to  satellite  (radians) 

AZ 

1 

I 

Azimuth  angle  to  satellite  (radians) 

TIME 

1 

I 

Universal  time  (radians) 

DFLUX 

1 

I 

Difference  between  the  daily  value  and  the 
12-month  running  average  of  the  solar  flux 

U 

13  >76 

I 

Art  -y  containing  coefficients  used  for  the  foF2 
computation 

UM 

9  *49 

I 

Array  containing  coefficients  used  for  the 
M(3000)F2  computation 

OLAT 

1 

o 

Latitude  of  the  ionospheric  point  (radians) 

OLON 

1 

o 

Longitude  of  the  ionospheric  point  (radians) 

F0F2 

1 

0 

Critical  frequency  f0F2  (MHz) 

HM 

1 

o 

Height  at  the  maximum  electron  density  h, 

(meters ) 

HLAT 

1 

0 

Magnetic  latitude  of  the  ionospheric  point 
(radians ) 

e)  Common  blocks: 

:  none 

f)  File 

requirements:  none 

3.  2.  1.  4 

CPC  No.  4 

Data  Organization 

Variables  defined  in  data  statements; 

Name 

Dimension  Des 

cription 

K 

KN 

KM  10 

NFF 

NMF 

10 

10 

1 

1 

1 

Interger  indices  and  index  arrays  used  for  the 
’  computation  of  fQFZ  and  M(3000)F?.  in  SUBROUTINES 
DKSICO,  GK  and  DKGK 

R 

1 

Mean  earth  radius  (meters) 

SPLAT 

1 

Sine  function  of  the  geographic  latitude  of  the  magnetic 
north  pole 
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Name  Dimension  De»criptlon 

CPLAT  1  Cosine  ox  the  geographic  latitude  of  the  magnetic 

north  pole 

PLOIS  1  Geographic  longitude  of  the  magnetic  north  polo 

(radians) 

HI  1  Coefficient*  used  in  the  formula  expressing  h, 

H2  1  as  a  aecond  order  polynomial  of  M(3000)F2 

H3  1 

PER  1  Model  conatanta  uaed  for  adjusting  f0F2  for  daily 

CENT  3  variation,  dependent  on  the  dadby  value  and  the 

12-month  running  average  of  solar  flux  and  magnetic 
latitude 


Other  constants  listed  in  data  statementsi 

QU1,  01000  =  1000,  Q1P999  =  1.  999999,  Q3T6  =  3>  1 0®;  D180=180‘, 
DG(l)=594  ,  IiG(2)=28*,  DG(3)«-33*  converted  to  radians. 

Other  imporlant  variables  are  described  under  3,  2.  1.3  d). 


3.2.  1,5  CPC  No,  4  Limitations 

There  a  *e  no  program  restrictions  connected  with  this  subroutine, 
and  the  limitations  to  the  accuracy  of  thy  result*  obtained  from  the 
formulas  art  discussed  in  Section  6.  2. 
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CPC  No.  -  Flowchart,  SUBROUTINE  PROFL1 
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3.2.  1  Computer  Program  Component  5 

CPC  No.  5,  SUBROUTINE  PROFL2,  if  written  in  FORTRAN  code.  It 
is  called  from  SUBROUTINE  REFRAC  and  computei  the  following  iono- 
•  pheric  profile  parameter*!  the  value*  of  half  thickness  y,,  yt  for  the 
bottomside  bi-parabola  and  the  topsido  parabola  respectively,  the  decay 
values  k^k,  ,  k,  for  the  topside  exponential  layers,  the  ratio  NT/N,  of  the 
total  content  ;o  the  maximum  electron  density,  and  the  multiplier  m  for  use 
in  the  range  .rate  computation, 

3.2.  1.1  CPC  No.  5  Description 

PROFL2  evaluates  the  ionospheric  profile  based  on  the  model  constants 
presented  in  graphic  form  in  Section  6,  1.  The  local  time  is  computed  from 
the  universal  time  t  and  the  longitude  \  of  the  ionospheric  point, 

t  1  OS  =  t  +  \  i 

The  half  thiccness  y,  of  the  bottomside  bi-parabola  varies  with  critical 
frequency  f0;r2  and  local  time.  Values  of  the  half  thicknes*  are  tabulated 
in  array  YMUAB  at  1  MHz  increments  for  foF2«»2,  3, ...  10  MHz  and  at 
2  hour  interval*  for  tloe  =0,  2, ...  22  hours.  To  obtain  y,  for  the  given 
condition*,  toe  table*  are  interpolated  in  two  dimensions  between  the 
fixed  values;  local  time  interpolation  is  carried  continuously  aero**  the 
0/24  hour  murk,  and  the  boundary  values  are  anumed  whenever  f0F2 
is  outside  tho  limits  2  and  10  MHz. 

For  »ea  tonal  adjustments  computation  of  the  parameter  Ay  (variable 
DSZA)  is  recuired.  Ay  is  the  deviation  of  the  daily  value  X  from  the 
yearly  average Y  of  the  noontime  solar  zenith  angle.  First  the  solar 
declination  6  is  evaluated  for  the  given  day, 

6  --  6„,  sin  [  il  (JDAY-80)  ], 

36  5 

6 =  23,  4444 ’  l*  the  maximum  solar  declination,  JDAY  is  the  day  of  the 
year,  For  stations  in  the  northern  hemisphere  and  outside  the  tropic*, 
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with  latitudes  >  23.4444*,  A X‘b\  for  stations  in  the  southern  hemisphere 
and  outside  the  tropics,  A\  =  -6.  In  the  tropics  the  yearly  average  of 
the  noontime  solar  zenith  angle  is  computed  as, 

X  =  —  -da  +<6  arc  sin  £ —  )  , 

r  *  »»* 

dbeing  the  latitude  of  the  ionospheric  point.  The  daily  noontime  solar 
zenith  angle  is  x  a|d-6l  »  and  the  difference  A6=7-X* 

The  half  thickness  of  the  bottomside  parabola  y,  is  multiplied  by 
a  seasonal  adjustment  factor  that  varies  with  Ay>  Ioc*l  time  and  magnetic 
latitude  dfc,  Adjustment  factors  are  tabulated  in  array  YRAT  at  8*  incre- 
ments  for  Ax-24,  16,  8,  0,  .8, « 16,  .24  de^’’*™*,  at  6  hour  intervals  for 
tl64  =5,  5,  1  1,  5,  17.  5,  24.  5  hours  where  the  absolute  value  of  the  magnetic 
latitude  is  greater  or  equal  15',  and  at  12  hour  intervals  tlc,  “3,  15  hours 
where  |d,|<  5*.  The  seasonal  adjustment  factor  for  the  given  conditions 
is  obtained  by  three  dimensional  interpolation;  the  local  time  interpolation 
is  carried  continuously  across  the  0/24  hour  mark  and  the  magnetic  latitude 
interpolation  is  only  performed  between  5  and  15  degrees. 

The  decay  constants  k^k^lo,  for  the  lower,  middle  and  upper  layer 
of  the  exponential  topside  are  related  to  the  daily  solar  flux  F  through 
the  first  order  polynomial, 

k(  a  S(  «  F  +  C(  ,  i  =  l,  2,  3. 

The  slopes  S,  ,  stored  in  array  SLOP,  and  the  intercepts  C,  in  array  CEPT 
uf  this  straight  line  relationship  vary  with  magnetic  latitude  d,  and  with 
f0F2.  For  each  of  the  three  topside  layers,  S,  and  C,  are  tabulated  at 
30°  Intervals  for  |da|  =  15,  45,  75  degrees,  and  at  3  MHz  increments  for 
f0F2  =  2,  5,  8,  11  MHz,  To  obtain  the  decay  constants  for  the  given  conditions, 
the  tables  for  S,  and  C$  are  interpolated  in  two  dimension  between  the 
fixed  values, and  whenever  f0F2  is  outside  the  limits  2  and  11  MHz  or  id.,! 
is  outside  15  and  75  degrees,  the  boundary  values  are  used, 
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Seasonal  effect*  are  imposed  on  the  topside  by  multiplying  the  decay 
conetante  by  leaaoa  adjustment  factors  that  vary  with  the  deviation  Ax  in 
the  eolar  aenith  angle  and  with  local  time,  The  adjustment  factor!  are 
tabulated  in  array  RAT  It  for  each  of  the  thr**  topside  layers  at  8*  increments 
for  Ax  =24,  16,  8,  0,  -8,  -16,  -24  degrees,  and  at  6  hour  Intervals  for 
tt  0  0  =2,  8,  14,  20  hours.  They  are  interpolated  f a*  each  kt,  1=1,2,  3  in  two 
dimensions  to  the  given  conditions;  the  local  Urn*  interpolation  is  carried 
continuously  across  the  0/24  hour  mark. 

The  half  thickness  of  the  topside  parabola,  extending  from  the  point 
of  maximum  electron  density  to  the  lower  exponential  layer,  is  dependent 
on  y,  and  f0I'2  through  the  relationship, 


Yx  * 


y.  ,  for  f0F2  <  10.  5 

y,  [1  +  C.  133333  (foF2-10.5)]  ,  for  f0F2  >  10.5 


The  distance  d  above  the  height  at  maximum  electron  density  h,  where  the 
slopes  of  the  parabola  and  the  lower  exponential  layer  are  the  same  is, 

d  .  -L  U  i+k?y=  -  1  )  . 


The  total  vertical  election  content  Nt  is  obtained  by  integrating  the 
electron  demity  profile  from  zero  to  the  height  of  the  satellite  hs .  The 
program  computes  the  ratio  of  total  electron  content  to  the  maximum  electron 
density  NT  /M,  (variable  XNTNM)  by  one  of  the  following  six  equations  de¬ 
pending  on  tlie  upper  integration  limit.  At  the  samo  time,  the  multipler 
m  (variable  R.RM )  required  for  the  instantaneous  range  rate  computation 
is  evaluated  and  its  formulation  also  varies  depending  on  the  height  of  the 
satellite,  1  ir  a  satellite  below  the  bi-perabolic  layer  of  the  icnospherei 

N, /N,  =  0  . 
m  *  0  . 
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For  a  satellite  in  the  bottomside  bi-parabolic  layer  with  half  thickness  ya  t 


3  N»  /  y*  *  (h, 


,  -h, )  + 1  i?  .  1  } . 


For  a  satellite  in  the  topside  parabolic  layer  with  half  thickness  yt: 


M.  a  M  f  ®  ..  M.  l_  v  .  1  (k»  )  1 

INt  i§  y*  -  (h,  -hj)  +  y  —^jj -  ; 


For  a  satellite  in  the  lower  exponential  i-yer  of  the  topside  with  decay 
constant  kj  t 

N,  •  N,(l  ■  -£-)  l-ji-  (l  -  .-k'(h‘-ho,)}  +  N.  , 
and  the  height  of  the  bottom  of  the  lower  exponential  layer  is  ho  =  h,  +d,  and 


N'  =  N,  {  TT  Vs  -  >’*1.  -h,,  ) 


i  ^.-^1 

3  y2  * 

>  t 


(  ,  #1*  \ 

m  =  (  1  -  -h—  je 

y? 


da  '  -K  (hj  -*k>  ) 


For  a  satellite  in  the  middle  exponential  layer  of  the  topside  with  decay 
constant  kg  : 


Nt  =N.  (l  -  ~  ){_!+  e'kl(hi‘ho)  [.  _L  +  _J_  (l  .  e-Vh»-bi  )yjJ+N,  , 

y%  hi  hj  kj 
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and  the  height  of  the  bottom  of  the  middle  exponential  layer  is; 
hj  »  he  4-  -i-  (1.012  *  10® -h©  ). 

m-  (l  -  . 

y? 

For  a  aatellite  in  the  upper  exponential  layer  of  the  topside  with  decay 
constant  kg : 

N,  .N,(l  .  *)(J_  +  .■*.(*.*.)[.  1  +  1 

and  the  heigh:  of  the  bottom  of  the  upper  exponential  layer  is, 

hg  =  ho  +  L  (1.  012  «  1 0®  -  h<j  ). 

m=  (l  .  )e-M^-^)e-Mh.-h1)e.ka(h,.h,) 

y* 

3.  2,  1.  2  CPC  No.  5  Flowchart 

The  flowchart  is  shown  on  the  page  following  3,  2.  I,  5. 

3,2,  1.3  CPC  No.  5  Interfaces 

4 

a)  Library  subprogram  required)  ABS,  AMOD,  ATAN,  EXP,  SIN,  SQRT 

b)  Other  subprograms  called)  none 

c)  Calling  program)  SUBROUTINE  REFRAC 

d)  Calling  nequencei  CALL  PROFL2  (OLAT,  OLON,  HS,  TIME,  IDAY,  MON, 
FLUX,  F0F2,  HM,  HLAT,  YM,  YT,  XK,  RRM,  XNTNM) 


Variables  in  calling  sequences 


Name 

Dimension 

I/O 

Description- 

OLAT 

1 

I 

Latitude  of  ionospheric  point  (radians) 

OLON 

1 

I 

Longitude  of  ionospheric  point  (radians) 

HS 

1 

1 

Height  of  satellite  above  earth's  surface  (meters) 

TIME 

1 

I 

Universal  time  (radians) 

IDAY 

1 

I 

Day  (»1  through  31 ) 

MON 

1 

I 

Month  ( a  1  through  12  ) 

FLUX 

1 

I 

Daily  solar  flux  value 

F0F2 

1 

I 

Critical  frequency  (MHe) 

HM 

1 

I 

Height  at  maximum  electron  density  (meters) 

HLAT 

1 

I 

Magnetic  latitude  of  ionospheric  point  (radians) 

YM 

1 

0 

Half  thickness  of  the  bottom  bi-parabolic  layer 
(meters ) 

YT 

1 

0 

Half  thickness  of  the  topside  parabolic  layer 
(meters ) 

XK 

3 

0 

Decay  constants  for  lower,  middle  and  upper 
section  of  the  topside  exponential  layer  (1 /motor) 

RRM 

1 

0 

Multiplier  of  the  li  term  i  n  the  range  rate 
formula  (dimensionless) 

XNTNM 

1 

0 

Ratio  of  total  vertical  electron  content  to  the 
electron  density  (meters) 

e)  Common  blockn  none 
I)  File  requirements1  none 

3,  <1,  1,4  CPC  No.  5  Data  Organisation 

Variables  defined  in  data  statement*! 

Name  Dimension  De scriptlon 

501  1  Maximum  solar  declination  (radians) 

502  1  Multiplier  to  convert  365  days  to  2tt  radians 

(radians/day) 
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Name  Dimension  Description 

RN4  1  Average  frequency  to  which  topaide  sounders  me&eured 

the  ionospheric  profiles  is  RN4*fQF2 

H1012  i  Average  height  to  which  topside  exponential  layer 

was  modeled  (meters) 

CEPT  4*3*3  )  Model  constants  used  for  computing  the  decay  constants 

SLOP  4*3x3  for  the  lower,  middle  and  upper  section  of  the  topaide 

1  exponential  layer,  dependent  os  daily  solar  flux, 
critical  frequency  and  magnetic  latitude 

RATK  4*4*3  Model  constants  used  for  adjusting  the  computed  decay 

constants  for  the  lower,  middle  and  upper  exponential 
topsides  for  seasonal  effects,  dependent  on  the  differ, 
once  between  the  yearly  average  and  the  daily  value  of 
the  noontime  solar  zenith  angle  and  on  local  time 

YMTAB  12*  9  Model  constants  used  for  computing  the  half  thickness 

of  the  bottomside  bi. parabola,  dependent  on  local  time 
and  critical  frequency 

YRAT  7*6  Model  constants  used  for  adjusting  the  computed  half 

thickness  of  the  bottomside  bi. parabola  for  seasonal 
effacts,  dependent  on  the  difference  between  the  year  iy 
average  and  the  daily  value  of  the  noontime  solar 
zenith  angle,  on  local  time  and  magnetic  latitude 

Other  constants  listed  in  data  statements  for  convenience! 

QO**0,  Ql  =  l,  Q2=2,  Q3*3,  Q4«4,  Q5  =  5,  Q6=6,  Q8  =  8,  Q24=24,  Q37  =  37, 
Q1000=1000,  QP05  =  ,  05,  Q PI 333  =  .  1  3 33 33,  OP95  =  ,  95,  Q2P5=2.  5, 

Q 1 0P5  =  1 0.  5 ,  Q8015  =  .  533333333; 

D5  =  5  * ,  D7P5-7,  5* ,  D8  =  8\  D10=10C,  Dl6  =  l6*.  D30=30’,  D135=135°, 
D180=180‘,  PIH  =  90*,  PI2  =  36  0* ,  DEG(1)=75‘,  DEG(2)=45',  DEG(3)=15° 
converted  to  radians, 

Other  impor:ant  variables  are  described  under  3,  2,  1,3  d). 

3.  2,  1,5  CPC  No.  5  Limitations 

There  are  no  programming  restrictions  connected  with  this  subroutine, 
and  the  limitations  to  the  accuracy  of  the  results  obtained  from  the  formvlas 
are  discussed  in  Section  6,  2, 


CPC  No.  5  Flowchart,  SUBROUTINE  PROFL2 
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43 


3,  2,  1  Computer  Program  Component  6 


CPC  No.  6,  SUBROUTINE  BETA(  is  written  in  FORTRAN  code.  It 
is  called  from  SUBROUTINE  REFRA.C  and  computes  the  ionospheric 
refraction  correction  for  the  elevation  angle. 


3,  1,  1  CPC  No,  6  Description 

BETA  computes  the  angular  refraction  correction  to  the  elevation 
angle.  Usinj  the  results  of  Mai iphant1  a  work  (Reference  8), the  deviation 
angle  a  is  expressed  as  the  angle  between  the  true  ray  path  above  the 
ionosphere  and  the  apparent  ray  path. 


a  = 


l_ 

2 


tan  $0  secaegc 


Nr. 

N, 


I 


where  f  is  thu  transmission  frequency,  f0F2  the  critical  frequency,  Nr 
the  total  electron  content,  N,  the  maximum  electron  density; 

r0  =  R,  +  h,  t  0.  5333  yt  , 

and  R  is  the  earth  radius,  h,  the  height  of  the  maximum  electron  density, 
and  y,  the  half  thickness  of  the  bottom  layer  of  the  ionosphere; 


arc  sin 


cosE 


,  E  being  the  elevation  angle, 


and  5  is  a  function  of  the  squared  deviation  factor  (iecc,  *i0F2/lf  and 
is  interpolated  from  tabulated  values  f"1  ;  cp.  =  arc  sinf*  — cosE  )  . 

After  determinating  athe  following  two  auxiliary  equations  are  evaluated, 


Xv  =  I  (R,  +h5  f  .  R3  cos3eJ^  +  R,  cosE  tan  — 

2 

a 

Xa  =  R<  iinE  -  Rt  co#E  tan  —  . 
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The  elevation  angle  correction  AE  is  then  given  by, 


AE 


arc  coe 


X:  cose.  -Xg _ 

lx?  +x=  ^x,  to«a)i 


3.  2.  1,  2  CPC  No,  6  Flowchart 

The  flowchart  is  shown  on  the  page  following  3,  2,  1.  5. 


3,  2,  1,  3  CPC  No,  6  Interfaces 

a)  Library  subprograms  required*  ABS,  ATAN,  COS,  SIN,  SQRT 

b)  C>ther  subprograms  call.-d;  none 

c)  Calling  program:  SUBROUTINE  REFRAC 

d)  Calling  sequence:  CALL  BETA  (FRAT,  XNTNM,  HS,  HM,  VM,  SE,  CE, 
DELEV) 


Variables  in 

caJling 

sequence: 

Name 

Dimen 

Ci.cn  I/O  Inscription 

FRAT 

1 

I 

Square  ratio  of  critical  frequency  to  the  trans¬ 
mission  frequency 

XNTNM 

1 

I 

Ratio  of  total  electron  content  to  the  electron 
density  (meter ) 

HS 

1 

I 

Height  of  the  satellite  above  the  earth's  surface 
(mete  rs  ) 

HM 

i 

I 

Height  of  the  maximum  electron  density  (meters) 

YM 

1 

I 

Half  thickness  of  the  bottom  layer  of  the  ionosphere 
(meters ) 

SE 

1 

I 

Sine  function  of  the  elevation  angle 

CE 

1 

I 

Cosine  of  the  elevation  angle 

DELEV 

1 

O 

Ionospheric  refraction  correction  to  the  elevation 
angle  ( radians  ) 

e )  Common  bio 

cks;  rone 

f)  File 

requirements; 

line  printer 

3.  2.  I,  4  CPC  No.  fe  Data  Organization 

Variables  defined  in  data  statements; 


Name 

Dimensions 

Description 

XAX 

5 

Values  of  the  squared  deviation  factor 

(sec  Cpc  *  f0F2/f  f  for  which  file  function  ?_1  is 

tubulated 

YAX 

5 

Tabulation  of  the  function  ?  1  as  given  in  Refe  ence  8 

R 

l 

Mean  earth  radius  (meters) 

Other  constants  listed  in  data  statements; 

Q0=0,  QU1,  Q2=2,  0.5333*.  5333 

Other  important  variables  are  described  under  3,2.  1.  3  d), 

3,2,  1,5  CPC  No.  6  Limitations 

The  equations  for  the  deviation  angle  a  which  are  coded  into  SUB- 
ROUTINE  BETA  are  accurate  everywhere  except  right  about  reflection 
conditions,  Whenever  the  deviation  factor  (sec  cp,  *  f0F2/f)  is  less  than 
0.  9,  all  equations  are  valid;  this  means  the  results  are  correct  whenever 
the  '•ompone  it  of  the  wave  frequency  vertical  to  the  ionosphere  is  slightly 
la.  ger  than  the  critical  frequency  ( 1 ,  i  »  fQF2 ),  But  the  more  the  deviation 
factor  excee  is  0.  9,  the  larger  the  errors  might  be  in  the  computation  for 
a  and  therefore  AE.  An  error  check,  programmed  into  the  routine,  tests 
if  the  deviation  facto  is  greater  than  0.  9  in  which  event  a  zero  elevation 
angle  correction  is  returned  and  an  error  mcsBage  is  printed. 


CPC  No. 


Flowchart,  SUBROUTINE  BETA 
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3.  2.  1  Computer  Program  Component  7 

CPC  No.  7,  SUBROUTINE  SICOJT,  is  coded  in  FORTRAN  code.  It  is 
called  from  SUBROUTINES  PROFL1  and  MAC  FIN  and  performa  auxiliary 
computation!  by  expressing  the  multiple  angle  trigonometric  functions. 

3.2,  1.1  CPC  No.  7  Description 

SICOJT  computes  the  trigonometric  functions  for  multiples  of  the 
angle.  It  forms  ein(jT ),  coa(JT  )  for  J«l. ,,,  L  by  computing  sinT,  co«T 
for  the  single  angle  T,  and  by  using  for  multiple  angles  the  recursive 
equations: 

sin  ( j+1 )  T  "]  =  sinT  cos(jT)+cosT  sin(jT) 
cos  pj+ 1 )  T  J  =  cosT  cos(jT  )  -  sinT  sin(jT). 

3.2.1.Z  CPC  No.  7  Flowchart 

The  flowchart  iJ  shown  on  the  page  following  3,  2,  1,  5. 

3.2.  1,3  CPC.  No.  7  Interfaces 

a)  Library  subprograms  required:  COS,  SIN 

b)  Other  sucprogratns  called:  none 

c)  Calling  programs:  SUBROUTINES  PROFL1  *nd  MAGFIN 

d)  Calling  sequence:  CALL  SICOJT  (L,  C,  S,  T  ) 

Variables  in  calling  sequence; 


Name 

Dimsnsion 

I/O 

Description 

L 

i 

I 

The  largest  integer  by  which  T  is  to  be  multiplied 

C 

L 

O 

Array  containing  values  cos (jT  ),  j  =  l . L 

s 

L 

o 

Array  containing  values  sin(jT  ),  j  =  l,  ,  .  .  ,  L 

T 

1 

I 

The  angle  (radians) 
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e)  Common  blocks;  none 

f)  File  requirements)  none 

3,2,  1,4  CPC  No,  7  Data  Organisation 

Important  variables  are  described  under  3.  2,  1,  3  d). 

3,  2,  1.  5  CPC  No,  7  Limitation* 


None. 


CPC  No.  7  Flowchart,  SUBROUTINE  SICOJT 
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3,2,  1  Computer  Program  Component  8 

CPC  No.  8,  SUBROUTINE  DKSICO.  la  written  in  FORTRAN  code. 

It  ia  called  from  SUBROUTINE  PROFL1  and  calculates  the  time  dependent 
coefficient*  which  are  required  for  the  computation  of  critical  frequency 
and  aaaociated  height, 

3,  2,  1,  1  CPC  No,  8  Description 

DKSICO  forma  the  orthonnrmal  coefficients  D*  for  a  fixed  time  T 
represented  by  the  Fourier  aeries  representation, 

M 

D,(T):Usl+  I  [u8j,k  coa(jT)  *  Ua  .  k  ain(jT)  ]  ,k=l,...,K. 

j  =  i 

These  coefficients  are  to  be  used  for  the  computation  of  the  ionospheric 
characteristics  in  DKGK.  The  number  of  harmonics  retained  in  the  series 
is  H,  higher  harmonics  are  not  considered  since  they  are  produced  more 
by  noise  than  by  real  physical  variation.  For  the  f0F2  computation  H=6 
and  for  the  M(3000)F2  computation  H*4  are  sufficient.  The  coefficients 
Ui  „  are  either  a  monthly  predicted  coefficient  set  for  M(3000)F2  or  a 
ten  day  predicted  coefficient  set  for  f0F2,  which  are  both  specific  subsets 
derived  from  the  generalized  f0F2  and  M(3QOQ)F2  coefficients  in  SUB- 
ROUTINE  REFRAC.  The  Dk  coefficients  are  computed  for  each  term  in 
a  series  with  cutoff  point  K,  K  =  75  for  the  series  expressing  f0F2  and 
K=48  for  the  series  representing  M(3QQO)F2, 

3.  2.  1,  2  CPC  No,  8  Flowchart 

The  flowchart  is  shown  on  the  page  following  3.  2,  1.  5. 

3,  2,  1,  3  CPC  No,  8  Interface* 
a)  Library  subprograms  required!  none 
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b)  Other  subprograms  called:  none 


c)  Calling  program:  SUBROUTINE  PP.OFL1 

d)  Calling  sequence,  CALL  DKSICO  (MX,  LH,  D,  SITIME,  COTIME,  DK) 

Variables  in  calling  sequence. 

Name  Dimension  I/O  Description 

MX  1  I  Cutoff  index  =  cutoff  poiat  of  series  4-1 

LH  1  1  Number  cf  harmonics  retained  in  Fourier  series 

representation  of  Dk 

D  (LH»2+1)«MX  I  Predicted  coefficient  array  U,  k  for  f0F2  or 

for  M(3000)F2  ' 

SITIME  LH  I  Array  of  values  sin(jT ) 

COTIME  LH  I  Array  of  values  cos(jT) 

DK  MX  0  Array  of  coefficients  Dk  at  fixed  time  T,  k=0,  .  .  .  ,  K 

e)  Common  blocks,  none 

l)  File  requirements,  none 


3,  2,  1. 4  CPC  No,  8  Data  Organization 


Important  variables  are  described  under  3,  2,  1,  3  d), 


3,  2,  1,  5  CPC  No.  8  Limitations 
None. 
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CPC  No,  8  Flowchart,  SUBROUTINE  DKSICO 


\ 
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3,  2,  1  Computer  Program  Component  9 


CPC  No.  9,  SUBROUTINE  MAGFIN,  It  written  in  FORTRAN  code. 

It  is  called  from  SUBROUTINE  PROFIT  and  evaluate*  the  magnetic  field 
component*  at  the  point  whore  the  wave  penetrate*  the  ionoaphere.  The 
field  component*  are  required  for  the  computation  of  the  critical  frequency 
and  the  aaaoc  i&ted  height. 

3.  2,  1.  1  CPC  No,  9  De*crlption 

MAGFIN  compute*  the  earth1*  magnetic  field  component*  at  a  desired 
location  following  the  spherical  harmonic  analysis  of  the  magnetic  field 
by  Chapman  ind  Bartel*  (Reference  3)  and  uaiug  the  coefficient*  g* ,  h* 
given  by  Jem  en  and  Cain  (Reference  4)  for  Epoch  I960.  The  X-north  , 
Y-eaat,  and  Z-vertical  (up)  component*  of  the  magnetic  field  are  u*ed 
for  the  computation  of  the  modified  magnetic  dip  in  SUBROUTINE  PROFL1, 

U*ing  the  specified  point  (0,X,1\,')<  the  colatitude  is  introduced 
55=90* -0 ,  and  the  ratio  R=R,  /(R,+h^  ),  where  R,  i*  the  radiu*  of  the  earth 
and  h^  =  300  km  i*  the  F2  layer  height  on  which  the  coefficient  analyni*  was 
baaed,  The  irigonometric  functions  *in(mX  ),  co*(mX  )  for  the  multiple 
longitude  angle  \  aru  computed  via  SUBROUTINE  SJCOJT,  The  magnetic 
field  component*  are  defined  in  the  following  equation*  and  are  obtained 
by  first  expressing  the  multiple  of  the  aiaociated  Legendre  function  and  it* 
derivative,  then  accumulating  the  term*  of  the  inner  sum*  and  finally 

forming  the  outer  sum*. 

8  r 

X  =  -  [  R''+*  Pn,  .  (co*ra)  rg^  coe(mX)  +  h*  •  in(mX)"|] 

n  si  t=0 

6  . 

Y  -  —•  |  R*+a  X  m  P  (co*cp)  rg!!  iin(m\)  -  hj  co»(mX  )  "I) 

aino  v  _  1  J 

n=  )  1  “-0 

6  ^  ^ 

Z  =  -  (n+1 )  R,,+s  \  P  (co*C5)  ^gj  co*(m\  )  +  h;  *in(mX  )  1  ^ 

k-«  *  <— i  J  J 

n  =1  i=0 
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The  multiple  of  the  associated  Legendre  function  is  given  by, 


3n..(co«cp)  =  «in*©  [co»n'*o  - 


+  (n-niXn-m.l  Hn-rn-Z  Xn-rn.3) 

(2)  (4)  (2n- 1 ) (2u.3)  » 


3.  2.  1.2  CPC  No.  9  Flowchart 

mm  — i  n— <n  mt  mi 

The  flowchart  is  shown  on  the  page  following  3,  2.  1.  5. 


3,2,  1.3  CPC  No.  9  Interfaces 

a)  Library  subprograms  required)  ABS,  SIGN,  SIN,  SQRT 

b)  Other  subprograms  called;  SUBROUTINE  SICOJT 

c)  Calling  programt  SUBROUTINE  PROFL1 

d)  Calling  sequence;  CALL  MAG  FIN  (POS,  UNE) 

Variables  in  calling  sequar.ee: 

Name  Dimension  I/O  Description 

POS  3  I  Array  containing  latitude,  longitude  and  height 

(radians  and  meters) 

UNE  3  O  Array  with  Z  (vertical  up),  X(north),  and 

y(oa8t)  components  of  magnetic  field  (gauss) 
at  the  location  specified  by  POS 

e)  Common  blocks;  none 

f)  File  requirements:  none 


3.  2,  1, 4  CPC  No,  9  Data  Organization 
Variables  defined  in  data  statements; 


Name 

Dimension 

Description 

CT 

7  *  7 

Array  containing  coefficients  for  the  computation  of 
the  associated  Legendre  function 

G 

7*  7 

Array  of  g*  coefficients  given  in  Reference 
the  earth  magnetic  field  for  Epoch  I960 

4  for 

H 

7  » 7 

Array  of  h*  coefficients  given  in  Reference 
earth  magnetic  field  for  Epoch  I96  0 

4  for  the 

I 


Name  Dimension  Pea cription 

RE  l  Mean  earth  radius  (meter*) 

Other  constants  listed  in  data  statement*) 

P(  1 ,  1  )=  1 ,  DP(l,l)*=0,  SP(1)*0,  CP(1)»1,  Q0*0;  RB99*B9.  9® 
converted  to  radians 

Other  important  variables  are  described  under  3,2,  1,  3  d). 

3,  2,  1,  5  CFC  No.  9  Limitations 
None. 


CPC  No.  9  Flowchart,  SUBROUTINE  MA.GFIN 
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3.  2.  1  Computer  Program  Component  10 


i 

i 

4 

1 


1 


CPC  No.  10,  SUBROUTINE  GK,  la  written  in  FORTRAN  code.  It  ia 
called  from  SUBROUTINE  PROFL1  and  calculate*  the  geographic  coordinate 
function*  which  are  required  for  the  computation  o(  critical  frequency  and 
aaaocinted  height, 


3,  2,  1,  1  CPC  No,  10  Deacrlption 

GK  compute*  the  geographic  coordinate  function*  G„  a*  a  function  of 
latitude  <P,  longitude  X,  and  modified  magnetic  dip  xsx(0,  X  ),  which  itaelf 
i»  dependent  on  the  geographic  po*ition,  These  coordinate  function*  are 
to  be  u*ed  for  the  computation  of  the  ionoapheric  characteriatic  f0F2  in 
aubroutine  DKGK.  The  function*  Gk  repreaent  the  main  latitudinal  variation 
and  the  first  order  through  8th  order  longitudinal  variation  terms.  The 
main  latitudinal  variation  ia  expressed  a*, 

G„  =  sin‘x  for  k=0,  l,  .  ,  .  ,  1 1, 
and  the  jth  order  longitude  terms  are  computed  as, 


(s:nx/l|,,i^  cos-'^  c  ■  h ( _, X) 
(s  nx/*  1 J  ‘l^8  cos- sin( jX  ) 


,  for  k  even 
,  for  k  odd 


k=m. ,  (mj  +  1  ) . (m 


1), 


1 


i  t 
t 


The  longitude1  orde  rs  are  j  - 1 ,  2,  ,  .  .  ,  8  while  k  =  1  2,  13,  ,  ,  ,  75 ,  and  the  index  - 

ing  is  defined  by:  m.  =12,  0^=36,  my -54,  m4 -64,  ms=68,  me  =  70,  m,  -72, 

m,  =  74.  i 

3.  2.  1.2  CP3  No,  10  Flowchart 

The  flowchart  is  shown  on  the  page  following  3.  2,  1.  5.  I 

i 
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3 .  L .  1 .  3  C  PC  Xa,  10  I nte  r  fa c e  s 


a)  Library  subprograms  required;  COS,  SL\ 

b)  Other  subprograms  called;  none 

c)  Calling  program;  SUBROUTINE  PROP  LI 

d)  Calling  sequence;  CA  LL  GK  (K,C,G) 

Variables  in  calling  sequence; 

Xanm  Dimension  I/O  Description 

K  10  I  Integer  index  array  containing  (mj  - 1 ) 

C  3  I  Array  containing  modified  magnetic  dip,  geographic 

latitude  and  longitude  (radians) 

G  76  O  Array  w;th  geographic  functions  Gx  ,  k  =  Q,  ,  .  .  75 

e)  Common  blocks ;  none 

f)  File  requirements;  none 

3.2,  1,4  CPC  No,  10  Data  Organization 

Constants  defined  in  data  statements; 

Q  1  =  1,  N  =  8  =  Highest  order  of  longitude  included  in  Gk  computation, 

important  variables  arc  described  unde r  3.2.  1.3  d), 


3.  2.  1.3  CPC  No.  1  0  Limitations 


CPC  No.  10  Flowchart,  Sl'HROL'TINE  GK 


_ }l _ 

[compute  terms  1 
•  of  due  to 

main  latitudinal  ; 
variation  j 

_ i. _ 

Prepare  loop 
over  8  orders 
of  longitudinal 

variation 

- 1 - 
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3,2.  1  Computer  Program  Component  11 


CPC  No.  11,  SUBROUTINE  DKGK,  is  written  in  FORTRAN  code,  it 
is  called  from  SUBROUTINE  PROFL1  and  compute!  the  critical  frequency  or 
the  associated  height  depending  on  the  input. 

3.  2.  1.1  CPC  No.  11  Description 

DKGK  computes  the  ionospheric  characteristic  fl,  by  forming  a  series 
of  products  of  time  dependent  coefficients  Dk  and  position  dependent  geo- 
graphic  functions  Gk  , 

-  (*,X,T)  =)  Dk  (T)  G„  )  . 
x  -  o 

The  coefficients  D„  are  precomputed  for  a  fixed  time  T,  and  the  geographic 
functions  GK  are  for  a  fixed  latitude  and  longitude  \ ,  K  is  the  cutoff  point 
for  the  approximate  series  representation  of  Q.  For  the  determination  of 
the  ionospheric  characteristic  Ci  =  fQF2  the  cutoff  point  K-75  is  used  and 
for  the  calculation  of  n  =  M(3000)F2  the  cutoff  point  is  K=48,  The  inputs 
Dk  and  Gk  are  specifically  set  for  either  the  f0F2  or  the  M(3000)F2 
computation. 

3 .  2 .  1 .  2  C PC  No.  11  Flowcha rt 

The  flowchart  is  shown  on  the  page  following  3.  2.  1.  5. 

3.2.  1.3  CPC  No.  11  Interfaces 

a)  Library  uuoprograms  required;  none 

b)  Other  subprograms  called;  none 

c)  Calling  program;  SUBROUTINE  PROFL1 

d)  Calling  sequence;  CALL  DKGK  (MX,  G,  DKSTAR,  OMEGA) 


6  1 


Variables  in  calling  sequence: 


Name 

Dimens  ion 

I/O 

De  s  c  ription 

MX 

1 

I 

Cutoff  index=cutoff  point  K  of  Beries  +1 

G 

MX 

1 

Array  of  geographic  functions  Gk  ,  k  =  0,  .  .  .  ,  K 

DKSTAR 

MX 

I 

Array  of  coefficients  Dk ,  k  =  0,  ,  .  .  ,  K 

OMEGA 

1 

0 

Ionospheric  characteristic  f0F2(MHz)  or  M(3000)F2 
(dimensionle  s  s ) 

e)  Common  blocks*  none 

f)  File  requirements:  none 

3,  2,  1,  4  CP3  No.  1 1  Data  Organization 

Important  variables  are  described  under  3,  2,  1.  3  d). 

3,2,  1,5  CFC  No,  1 1  Limitations 
N  one . 
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CPC  No.  11  Flowchart,  SUBROUTINE  DKGK 
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3.2.1  Computer  Program  Component  12 

CPC  No.  12,  main  PROGRAM  TABGEN,  is  written  in  FORTRAN  code. 
For  any  specified  date  and  station  preprocessor  TABGEN  computes  values 
of  critical  frequency  and  corresponding  height  for  14  time  intervals  at  each 
of  25  locations  around  the  station  covering  the  visible  ionosphere.  The 
resulting  f0F2-h,  tables  are  written  onto  file  for  use  in  the  ionospheric 
reduction  program  ION1. 


3.  2.  1.1  CPC  No.  12  Description 

TABGEN  reads  the  date,  station,  and  solar  flux  information  from  card 
for  which  f(.Ir2-hi  tables  are  to  be  generated,  It  lists  the  input  data  for 
reference  in  the  print  out  and  converts  the  units  of  the  angles  to  radians. 

The  genera,  coefficients  are  read  from  tape  if  not  already  available  and 
the  specific  coefficient  sets  required  for  the  fQP2  and  M(3000)F2  computa¬ 
tion  are  prepared  as  well  as  the  solar  data.  The  applicable  procedures  are 
already  described  in  the  first  four  paragraphs  of  Section  3.  2,  1.  I,  CPC  No.  2, 

A  pattern  of  25  points  is  generated  around  the  station  as  shown  in 
Figure  1;  tie  point  distribution  covers  the  visible  ionosphere  in  fairly 
even  density.  The  earth  central  angle  a  between  station  and  ionospheric  point 
varies  in  7*  increments,  while  the  azimuth  A  is  0°  for  a=0‘,  and  rotates  in 
90°  steps  fc  r  a  =  7  ° ,  in  45”  steps  for  a  -  14 0 ,  and  in  30°  steps  for  a  -  2  1  5  out 
of  the  northern  position.  For  each  ionospheric  point  the  geographic  latitude 
1 1  and  longitude  \  and  the  magnetic  latitude  0t  are  reduced  from  the  station 
position  ,  \j ,  the  position  of  the  magnetic  north  pole  ,  K .  and  a  and  A; 
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The  position  dependent  functions  required  for  the  f0F2  and  M(3000)F2 

computation  are  evaluated  using  SUBROUTINE  MAGFIN  and  GK  at  described 
in  the  third  paragraph  of  Section  3.  2,  1.  1,  CPC  No,  4, 

The  diurnal  variation  at  each  of  the  25  points  is  produced  by  evaluating 
the  critical  frequency  and  corresponding  height  at  14  different  time  intervals 
at  0,  2,  4,  5,  6,  7,  8,  10,  12,  14,  16,  18,  20,  22  hours  of  local  time.  The  time 
pattern  is  densified  around  sunrise  to  properly  represent  the  rapid  change 
in  the  ionosphere  during  that  time.  f0F2  and  M(3000)F2  are  computed 
by  preparing  the  time  dependent  coefficients  via  SUBROUTINES  SICOJT 
and  DKSICO  and  combining  the  time  dependent  coefficients  and  position 
dependent  functions  by  calling  SUBROUTINE  DKGK,  The  height  of  the 
maximum  electron  density  h,  is  computed  with  the  Appleton  -  Beynon  equation 
(Reference  1)  in  units  of  km; 

h,  =  1346.  92  .  526. 40  «  M(3000)F2  +  59.  825  [m(3000)F2  }  , 

The  critical  frequency  is  adjusted  for  day  to  day  variations  as  a  function 
of  AF,  the  d.fferonce  between  the  daily  value  and  the  12-month  running 
average  of  the  solar  flux.  Using  the  model  constants  c;  (variable  PER) 
and  cs  obtained  by  interpolating  the  constant  table  (array  CENT  )to  the 
magnetic  la'  itude  of  the  ionospheric  point,  f0F2  is  multiplied  by  the  ad¬ 
justment  factor  (c;/*F*ca), 

Fur  eac  h  point  and  time  f0F2  and  h,  are  coded  into  one  8  digit  integer, 
the  first  four  digits  defining  h,  in  units  of  km  the  last  4  digitB  specifying 
f0F2  in  unit  )  JoQ  MHz.  The  f0F2-h,  table  ib  accurm.  ated  for  all  14  time 
intervals  and  all  25  points,  and  is  written  to  tape  or  uiic  file  along  with  the 
date,  station,  and  solar  flux  information.  The  process  can  be  repeated  for 
any  number  of  date  and  station  conditions  desired.by  specifying  additional 
input  data  and  repeating  the  steps  outlined  in  this  section. 
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3.2.  1.2  CPC  No.  12  Flowchart 


The  flowchart  la  shown  on  the  page  following  3,2,  1,  5, 

3,  2,  1.3  CPC  No.  12  Interfacea 

a)  Library  subprograms  required!  AMOD,  AT  AN, COS,  SJN,  SQRT 

b)  Other  subprograms  called!  SUBROUTINES  MAGFIN,  GK,  SICOJT, 
DKSICO,  DKGK 

c)  Calling  programs!  none 

d)  Calling  sequence;  PROGRAM  TABGEN 

e  )  Common  blocks;  none 

f)  File  requirements:  general  coefficient  input  tape,  output  disc  or  tape 
file  with  f0F2-h,  tables,  card  reader,  line  printer 
The  formats  of  the  general  coefficient  input  tape  of  the  fQF2-h,  table 
output  file  and  the  requirements  for  the  input  data  card  file  are  speci¬ 
fied  under  3.  3.  1, 

3,  2,  1,  4  CPC  No,  12  Data  Organization 
Variables  defined  in  data  statements! 


Name 

Dimension 

Description 

JAZ 

4 

Index  array  specifying  number  of  azimuth  angle  divisions 
for  each  earth  central  angle  used  in  25  point  pattern 

IT  P 

1 

Unit  assignment  of  general  ionospheric  coefficient  tape 

JTP 

1 

Unit  assignment  of  file  with  f0F2-h,  tables 

MONDY 

1  I 

Initialization  constants  for  last  and  first  (month*  100+day) 

MOND 

1  | 

for  which  coefficients  are  in  core 

LYRMO 

1 

Initialization  constant  for  (years  100+month) 

K 

10 

KN 

10 

Integer  indices  and  index  arrays  used  for  the  computa¬ 

KM  10 

1 

tion  of  f0F2  and  M(3000)F2  in  SUBROUTINES  DKSICO, 

NFF 

1 

GK,  and  DKGK 

NMF 

1 
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Name 

Dimension 

Description 

PER 

1 

Model  constants  used  foe  adjusting  f0F2  for  daily 

CENT 

3 

variation,  dependent  on  the  daily  value  of  the  12- 
month  running  average  of  solar  flux  and  magnetic 
'  latitude 

SPLAT 

1 

Sine  function  of  the  geographic  latitude  of  the  magnetic 
north  pole 

CPLAT 

1 

Cosine  of  the  geographic  latitude  of  the  magnetic  north 
pole 

PL  ON 

1 

Geographic  longitude  of  the  magnetic  north  pole  (radians 

HI 

H2 

H3 

1  ] 
1 

1 

Coefficients  used  in  the  formula  expressing  h„  as  a 
second  order  polynomial  of  M(30Q0)F2 

Other  constants  listed  in  data  statements: 

Q 1  =  1 ,  Q 1 0-1(1,  Q100=100,  Q130=130,  Q3T5  =  3<10®,  QP1=.  1,  QP5  =  ,  5; 

DR=  1  * ,  PI  2  =  3  6  0° ,  D7  =  7\  DHR1=1\  DHR2  =  2*,  0180  =  180*, 

DG(  1  }= 5  9  * ,  DG(  2  )=28  ° ,  DG(3)=.33°  converted  to  radians. 

3.2.  1.5  CPC  No,  12  Limitations 

The  daily  value  of  t  olar  flux  transferred  through  the  data  file  to  the 
ionospheric  reduction  program  for  computation  of  the  decay  constants  for 
the  topside  exponential  profile  is  truncated  a  i  maximum  value  of  130.  This 
is  the  boundary  that  was  imposed  by  the  data  base  during  model  development 
and  extension  of  solar  flux  beyond  130  could  reBult  in  invalid  profiles. 

Approximations  are  introduced  through  bypassing  the  iteration  on  the 
height  estimate  of  the  ionosphere.  In  this  case  the  latitude  and  longitude 
of  the  ionospheric  pointB  are  not  effected,  only  the  height  itself  at  which 
the  magnetic  field  components  are  evaluated,  Error  estimates  for  these 
approximations  are  not  yet  available. 

If  the  ionospheric  coefficients  are  not  found  on  the  tape  for  the  specified 
date,  an  error  condition  has  occurred,  a  message  is  printed  out,  and  the 
program  is  tirminated. 
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The  solar  input  data  cards  are  checked  for  consistency  of  the  date 
and  if  disagreement  is  found,  a  message  is  printed  and  the  program  is 
terminated. 


CPC  No.  12  Flowchart,  PROGRAM  TAEGEN 
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PROGRAM  TABGEN  (continued) 
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3.  2.  1  Computer  Program  Component  13 


CPC  No.  13,  main  PROGRAM  I0N1,  is  written  in  FORTRAN  code. 

It  handles  the  card  input  and  the  printing  of  the  result*  tor  the  entire 
program.  ION1  transfers  the  input  conditions  through  common  /EVA  L  l /, 
and  by  calling  SUBROUTINE  REFRC1  receives  the  computed  profile  para- 
meters  and  refraction  corrections  through  common/COARl /. 

3.  2,  1,  1  CPC  No.  13  Description 

ION1  reeds  the  station,  satellite,  and  time  information  for  the  condition 
to  be  evaluated  from  cards.  The  input  data  is  converted  to  the  internal 
units  of  meters  for  distances  and  radians  tor  angles  and  times.  The  variables 
specifying  th  ;  evaluation  condition  are  transferred  through  common/E VAI.  1  / 
to  SUBROUTINE  REFRC1,  Through  REFRC1  and  other  routines  called  by 
REFRC1,  ionospheric  profile  parameters,  vertical  and  angular  electron 
content,  refraction  corrections  to  elevation  angle,  range,  and  instantaneous 
range  rate  a:'e  computed,  They  are  roturned  to  ION1  through  comroon/CORRl  / 
and  are  prinled. 

Any  number  of  evaluation  conditions  can  be  processed  by  supplying 
additional  inout  data  and  repeating  the  program  steps  outlined  above.  Fn” 
more  details  about  the  input  and  output  data  refer  to  the  file  descriptions 
unde  r  3.  3,  1 , 

3.  2.  1.  2  CPU  No,  13  Flowchart 

The  flowchart  is  shown  on  the  page  following  3.2.  1,  3, 

3.  2.  1.3  CP3  No,  13  Interface 

a)  Library  subprograms  required;  none 

b)  Other  subprograms  called;  SUBROUTINE  REFRC1 

c )  Calling  programs;  none 
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d)  Calling  seiiu'jnu;:  PROGRAM  ION  1 

e)  Common  clocks;  EVAL1,  CORR1 

Variables  in  Common; 

See  description  for  EVAL1,  CORR1  unde  r  SUBROUTINE  REFRC1, 

CP'-  No.  14 

f)  I'ile  requirements;  card  reader,  line  printer 

The  requirements  for  the  input  data  card  file  are  specified  under  3.  3.  1, 

3.2.1.-  CT  No,  1 3  Data  Organization 

Con  tants  defined  in  data  statement; 

QQ  0,  Qi  D00--1C00,  Q3600-3600;  DR  =  1  ° .  MR-1"  converted  to  radians. 
Important  variables  are  described  under  3.2.  1.  3  e)  of  SUBROUTINE 
RE  PRC  1,  CPC  No.  14. 

3,  L,  ' M  CPC  No.  13  I. imitations 

trr  ir  tests  on  the  sequence,  units,  and  formats  of  the  input  data  are 
not  performed.  However,  mistakes  in  the  set  up  of  the  card  deck  arc 
r  mealed  in  the  pri  .out  of  the  input  data  that  is  listed  along  with  the  results. 

ION  1  is  a  program  for  special  applications  and  limited  use  compared 
to  the  general  purpose  PROGRAM  ION.  Not  inc  luded  in  ION  1  arc  the  addi¬ 
tional  matures  <>!  ION  of  plotting  the  ionospheric  profile,  of  updating  the 
predictions  with  actual  ionospheric  observations,  and  of  computing  range  rate 
corrections  for  range’  differencing.  For  the  purpose  of  saving  space  only 
four  digits  are  carried  lor  fQF2  and  h,  in  the  h^FZ-h,  tables  which  eliminates 
th*>  option  of  differencing  range*  cur  lections  where  the-  Sth  and  6th  digit  of 
fnF2  aie  significant  to  the  result.  Because  of  approximations  in  TABGEN 
and  KKFRfl,  I  ON  1  also  yields  less  accurate-  results  than  ION . 
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3.2.  1  Computer  Program  Component  14 


CPC  No.  14.  SUBROUTINE  REFRC1,  is  written  in  FORTRAN  code  and 
is  called  from  main  PROGRAM  ION1,  REFRC1  extracts  the  foFZ-h,  tables 
from  tape  or  disc  file  and  interpolates  the  values  in  the  tables  to  the  specified 
position  and  1'me,  The  remaining  profile  parameters  are  obtained  via 
SUBROUTINE  PROFL2,  the  ionospheric  refraction  corrections  to  range  AR, 
to  instantaneous  range  rate  \A  are  computed  and  SUBROUTINE  BETA  provides 
the  elevation  angle  correction  AE. 

3,2.  1,  1  CPC  No.  14  Description 

REFRC  1  retrieves  the  f0F2-hs  tables  from  the  tape  or  disc  file  that  was 
prepared  by  the  preprocessor  TABGEN,  if  the  tables  for  the  given  evaluation 
condition  are  not  already  available.  Data  for  up  to  four  station  and  date 
combinations  can  be  kept  in  core  simultaneously  which  greatly  reduces  the 
10  requirements  for  data  reductions  where  a  few  stations  are  observing  inter, 
mittently.  In  addition,  if  new  data  is  requested,  it  automatically  replaces 
of  the  four  tables  the  one  having  been  in  core  for  the  longest  time. 

The  earth  central  angle  between  station  and  ionospheric  point,  the 
geographic  latitude  and  longitude,  and  the  magnetic  latitude  of  the  ionospheric 
point  are  computed  using  the  equations  shown  in  Section  3.  2.  1.  1,  CPC  No,  4. 
Luc  time  t,  is  computed  from  the  universal  time  t  and  the  longitude 

t;  -  .  =  t  +  X 

Critical  frequency  and  corresponding  height  are  extracted  from  the 
f0F2-h,  table  containing  data  for  14  time  intervals  during  the  specified  day 
at  each  of  2^  locations  covering  the  visible  ionosphere  around  the  given 
station.  A  linear  interpolation  process  is  used  in  three  dimensions,  in 
azimuth,  earth  central  angle,  and  local  time.  First  it  is  arranged  for  index¬ 
ing  purposes  that  azimuth  lies  between  0  and  3bO  degrees,  central  angle  be¬ 
tween  0  and  dO  degrees,  and  local  time  between  0  and  24  hours.  The  indices 
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and  increments  fur  the  interpolation  are  computed  for  all  three  variables. 

Continuous  interpolation  is  insured  between  22  and  0  ho  rs  of  local  time 
and  between  the  highest  value  and  0  degrees  of  azimuth  for  each  central 
angle.  The  limiting  values  at  21  degrees  are  used  if  due  to  some  rare 
occasion  or  an  error  condition,  the  earth  central  angle  should  exceed  21  degrees 
this  value  wan  arrived  at  for  the  extreme  condition  of  an  observer  iooking 
horizontally  it  a  453  km  high  ionosphere. 

By  calling  SUBROUTINES  PR0FL2  and  BETA  the  remaining  profile 
parameters  and  the  refraction  correction  to  the  elevation  angle  are  evaluated 
respectively.  Vertical  and  angular  total  electron  content  as  well  as  the 
refraction  corrections  to  range  and  instantaneous  range  rate,  are  computed 
following  the  description  m  the  last  five  paragraphs  of  Section  3.2,  1.  1, 

CPC  No.  2. 


3.2.  1.2  CPC  No.  14  Flowchart 

The  flowchart  is  shown  on  the  page  following  3.2.  1.  5. 

3.  2,  1,3  CPC  No,  14  Interfaces 

a)  Library  subprograms  required;  ABS,  AMOD,  AT  AN ,  C  OS ,  F  LOAT ,  SIN  , 
SORT 

b)  Other  subprograms  called;  SUBROUTINES  PROF  i.2,  BETA 

c)  Calling  program;  PROGRAM  ION  1 

d)  Calling  lequence;  SUBROUTINE  REFRC1 
e  )  C  ommon  blocks:  EVAL1 ,  CORR1 

Variables  in  common; 

Common  Variable  Dimension  I/O  Description 
Name  Name 

EVAL1  F!i  1  I  Transmission  frequency  (MHz) 

EVA1.1  F  ^AT  1  I  Latitude  of  station  (radians) 

EVAL1  F  _ON  1  I  Longitude  of  station  (radiBr.s) 
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C  ommon 

Name 

Variable  Dimension 
Name 

I/O 

Desc  ription 

EVAL1 

ELEV 

1 

I 

Elevation  to  satellite  (radians  ) 

EVAL1 

AZ 

1 

I 

Azimuth  to  satellite  (radians) 

EVAL1 

HS 

1 

I 

Height  of  satellite  (m) 

EVAL1 

ED  or 

1 

I 

Elevation  rate  ( radians /sec  ) 

EVAL1 

HDOT 

1 

I 

Altitude  rate  (m/sec) 

EVAL1 

TIME 

1 

I 

Universal  time  (radians) 

EVAL1 

IYR 

1 

I 

Year  (last  2  digits  ) 

EVAL1 

MON 

1 

1 

Month  (  =  1  through  12) 

EVAL1 

IDAY 

1 

I 

Day  (  =  1  through  31 ) 

EVAL1 

JTP 

1 

I 

Unit  assignment  of  ionospheric 
file  with  fQF2-hB  tables 

CORR1 

DRANG 

1 

0 

Range  correction  (m) 

CORR1 

DRATE 

1 

o 

Range  rate  correction  (m/sec) 

CORR1 

DELEV 

1 

o 

Elevation  angle  correction  (radians) 

CORR1 

F0F2 

1 

0 

Critical  frequency  (MHz) 

CORR1 

HM 

1 

o 

Height  at  maximum  electron  density 
(mete  rs  ) 

CORR1 

YM 

1 

o 

Half  thickness  of  the  bottomside  bi¬ 
parabolic  layer  (meters) 

C  ORR 1 

YT 

1 

0 

Half  thickness  of  the  topside  parabolic 
laye  r  (mete  rs  ) 

CORR1 

XK 

3 

o 

Decay  constants  of  lower,  middle, 
and  upper  Bection  of  the  exponential 
topside  layer  (1/meter) 

C  ORR  1 

TOTN 

1 

0 

Total  vertical  electron  content 
(e/ms  column) 

CORR1 

TOTNA 

1 

o 

Total  angular  electron  content 
(e/m*  column) 

f)  File  requirements:  ionospheric  input  tape  or  disc  file  with  f0F2-h, 
tablci.  The  format  of  the  file  containing  the  f0F2-h,  tables  is 
described  under  3.  3,  1, 
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3,  2,  1.  4  CPC  No.  14  Data  Organization 


Variables  defined  in  data  statements: 


Name 

Dimension 

Description 

JAZ 

4 

Index  arrays  used  in  defining  the  25  point  pattern 

KAZ 

4  1 

around  the  station 

NO 

l  ! 

Initialization  constants  for  storage  condition  of 

NR 

1  | 

f0F2-ht  tables 

R 

I 

l 

Mean  earth  radius  (meters) 

SPLAT 

l 

Sine  function  of  the  geographic  latitude  of  the  magnetic 
north  pole 

C  PLAT 

l 

Cosine  of  the  geographic  latitude  of  the  magnetic 
north  pole 

PL  ON 

l 

Geographic  longitude  of  the  magnetic  north  pole 
(radians ) 

RM 

l 

Estimate  for  radial  distance  of  ionosphere  from 
earth  center  (meters) 

TOL 

i 

Tolerance  allowed  in  identifying  station  latitude 
and  longitude  (radians) 

Other  constants  listed  in  data  statements* 

00=0,  Q 1 

=  1,  Q2  =  2,  C 

-7,  0100=100,  Q3P5  =  3,  5,  Q4P5=4,  5,  QNM=  1 .  24« lO10 , 

RN3  =  ,  49972;  PI2-360”,  DR=1°,  HR  =  lr  convc rted  to  radians . 

Othe r  impo riant  variables  are  described  under  3,2,  1,3  e). 

3 .  2 .  1 ,  3  CP' 3  No,  14  Limitations 

Approxi  nations  are  introduced  into  the  computation  of  the  critic  al 
frequency  and  the  height  of  the  maximum  electron  density  by  two  facts; 
through  the  linear  interpolation  in  space  and  time  of  the  precomputed 
f0F2-h,  tabljs,  and  through  bypassing  the  iteration  on  the  height  estimate 
of  the  ionosphere.  Thus  caution  should  be  uged  ana  further  tests  of  accuracy 
requirement!  might  be  desired  when  using  this  program  version,  An 
estimate  of  «he  expected  errors  is  given  In  Section  6,  2. 
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Thu  range  rate  correction  formula  in  this  routine  applies  only  to  i-  ■■  r  '• 

eous  range  rate  measurements  Bince  it  is  assumed  that  the  only  variation  in 
electron  content  over  the  time  of  observation  io  due  to  the  positional  change 
of  the  satellite  and  that  the  ionosphere  between  station  and  satellite  remains 
constant.  Range  rate  corrections  to  observations  obtained  by  range  differ 
over  a  finite  time  interval  during  which  the  ionosphere  can  undergo  distinct 
changes,  cannot  be  computed  by  this  routine  because  the  f0F2-h,  tables  do 
not  carry  enough  significant  digits.  For  this  purpose  PROGRAM  ION  should 
be  used. 

If  the  f0F2-hs  table  for  the  specified  date  and  station  is  not  found  in  the 
data  file,  an  error  message  is  printed  out  and  control  is  transferred  to 
PROGRAM  ION1  to  proceed  with  the  next  data  case. 
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CPC  No.  14  Flowchart,  SUBROUTINE  REFRC  1 


Ente  r 


Initialize 

constants 


>les  in  core  fm 


Read  file  with  j 
,  foF2-hm  / 

\  tables  j 


^  proper  s 
f0F  2-hm  table 


found 


Print  error  j 
message  / 


Yes 


Fo i  m  azimuth 
&  earth  central 
angle  from  ata, 
to  on,  point _ 


Compute  geo. 
lat  tude  &  long. 
L  rnagn.  latitude 
for  ion.  point 


Compute 
local  time 


Compute  indices 
&  increments  for 
interpolation  in 
local  time, 
central  angle 
it  azimuth 


Pe rform  inter¬ 
polation  in  local 
time,  azimuth 
i  central  angle 


/Compute  2nd  ' 
part  of  profile 

\  ( PROFL2 )  / 


/  Compute  \ 
ilevation  angle  ' 
correction  AE/ 
\  (BETA)  / 


Compute  vertical 
K  angular 
electron  content 


Compute  cor¬ 
rections  to  range 
AR  and  to  inst. 
range  rati  \R 


Return 


i 


3.3  Storage  Allocation 


The  s  ize  requirements  and  storage  allocations  of  the  total  program 

and  the  individual  components  were  extracted  from  computer  runs  of  the 
programs  on  the  CDC  6600  computer  system,  In  the  load  maps  that  are 
shown  on  the  following  pages  the  starting  addresses  of  the  program  and 
system  functions  in  the  detailed  breakdown  are  listed  in  octal  words. 


The  total  core  space  requirements  are: 

37604  octal  =  1626  0  decimal  words  for  the  Bent  Ionospheric  PROGRAM  ION; 

24232  octal  =  10394  decimal  words  for  the  preprocessor  PROGRAM  TABGEN, 

63, :4  octal  =  3412  decimal  words  for  the  reduction  PROGRAM  ION  1 

of  the  alternate  version  of  the  ionospheric  program. 


Following  are  the  size  requirement':  for  the  individual  components- 


C  omponent 


COMMON 

/EVAL/ 

COMMON 

/U  PDT  / 

COMMON- 

/CORR/ 

PROGRAM 

ION 

SUBROUTINE 

REFRAC 

SUBROUTINE 

PLOTNH 

SUBROUTINE 

PROFL1 

SU  B ROUTINE 

PROFL2 

SU  B  ROUTINE 

BETA 

SUBROUTINE 

SICOJT 

SUBROUTINE 

DKSICO 

SUBROUTINE 

MAG  FIN 

SUBROUTINE 

GK 

SU  b  ROUTINE 

DKGK 

C  OMM  ON 

/EVAL  1  / 

COMMON 

/CORR  1  / 

PROGRAM 

TABGEN 

PROGRAM 

ION  1 

SUBROUTINE 

REFRAC 1 

Size  in  decimal  words 
20 
57 
12 
3841 
5426 
366 
624 
1  085 
1  80 
6  4 
72 
520 
147 
39 
13 
12 
10394 
341'’ 

1988 
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Load  Map  for  PROGRAM  ION; 


-PR0GRA4-- 

ION 


REFRAC 


PLOTNH 
PROFL  1 
PROF  L  2 
BETA 
SICOJT 
OKS IC  C 
MAGFIN 
G< 

DKG* 

AUGOFR* 

ABSS 

AMOOS 

SXRNJ 

AwMCGE 

AL  CGI  PS 

AT  ANS 

E>  PJ 

E  X  c  E 

SINCCSE 

FQPTS 

SORTE 

GET9A 

SICi 

CCS$ 

SINS 

AT  ANE 

SYSTEM* 

I F  E  NC  r i 

INPUTFJ 

INPUT  C  i 

KOOEPi 

K  R  A  K  E  -  i 

OUTPTCi 

REWINDS 


--ACORESS- 

--LABELED - COHPCN-- 

000231 

eval 

C001C0 

UPCT 

000124 

rORR 

000215 

007622 

eval 

0001CC 

UPDT 

QCC124 

CORR 

000215 

0  22  314 

9  2 1 0  7  ? 

02425? 

026347 

02j6?3 

a  26  7  3  3 

027042 

031053 

0  2  0  276 

030345 

P 30360 

0332b! 

u  ?  0  3  7  0 

030374 

3  ’  3  4  3  3 

u 30465 

0  3  j  5  0  5 

u  2  054  7 

I30bl5 

07067? 

u  3  J  7 1 6 

0  3  3  7  <♦  0 

030757 

0  3  2  37  2 

032424 

02245b 

022637 

033547 

0  3  j  b  2  b 

034107 

0  1 4  2  7  6 

075652 

037436 

037532 
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Load  map  for  PROGRAM  TABGEN: 


-PuiOC-RAH - ACQRESS-  --LA9ELED- 

JAS6EA - -O-udlOQ -  - -  .  .  _ 


sirojT 

-OKSIG-C . -- 

MAGFIS 

<JK 

0243J2 

-0644  £2 - 

024542 
n  9 St tt  9 

DKGK 

025775 

A&S*  - 

024044  - - - 

AM001  ' 

u  26  047 

■S1W4 . 

026-0*4-  - - 

AT  ANS 

0  26060 

-HMCCSS- 

-  - 

SCRT  * 

026155 

026201 

GET3A 

026223 

-si-o-i 

026242  . 

coss 

027655 

s I nt  ..  _ . 

«««■ - - 

AT  ANE 

f  27741 

S»STikJ 

530022 

tNCFILi 

031032 

I-F  EN-EF  i 

531103 

INPUTDJ 

0  31162 

WWC4 

•63144?  ■  - - 

K0CEF<» 

0  3 1  p  7  2 

033206 

OUTOTBS 

03*772 

outctci 

035251 

R  E  H I  K  *  J 

035  3*E 

-COHMON-- 
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Load  map  for  PROGRAM  I0N1: 


‘-“I 


-PROGRAM--- 

-ACDRESS- 

--LAPELED - CCHMCh-* 

- rani - 

0001*4  . . . . 

£  V  AL1 

000100 

C0RR1 

0  0  Q  1 1 5 

- R-EFR&l  - 

0-06655  -  - 

--EV4L1 

000100 

C0RR1 

000115 

_ _ 

-0-12*41  . -  ■  - 

9ETA 

014656 

45S-5 

015142 

A  MOO  f 

01p145 

- -F4-&44-*  - 

015152  . - 

AT  AN! 

015155 

- *X-P«— . 

0-15175  - 

EXPE 

015237 

...  -SU4CCSC 

013305 

SORTS 

015262 

SORTS- 

019406 

GETBA 

015430 

- - 

•  04*  4-47  ■  -  . . 

COSS 

017062 

W44- 

017114 

ATANE 

017146 

-S¥STf4* 

&17227 

IEENCFS 

0  2  C  2  37 

.  ..  T  LClLt  C-tf 

1  rrr’trr  “  w 

63(316 

INPUTCS 

02(577 

KOOERS 

02(726 

KRAKE9S 

0  2  2  742 

OOTOTC4 

02*126 

REWINDS 

02*  222 

«4 
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3.  3,  1  Data  Bane  Characteristic  -  File  Description 

External  data  transfer  in  and  out  of  the  ionospheric  model  ION  is 
handled  through  three  files;  the  input  data  card  deck  read  in  PROGRAM 
ION  ,  the  input  ionospheric  coefficient  tape  read  in  SUBROUTINE  REFRAC, 
and  the  output  to  the  line  printer  is  written  in  PROGRAM  ION  and  SUB¬ 
ROUTINE  PLOTNH. 


Prog  ram 

File  Type 

Mode 

I/O 

Fortran 
Unit  No, 

Description 

Details 

Under 

ION 

T  ape 

B  IN- 

I 

1 

Ionospheric  coeff.tape 

3.  3.  1.  I 

ION 

Line  printer 

BCD 

0 

6 

Output  listing  from  ION 

3.  3.  1 .  2 

ION 

Care;  react: r 

BCD 

1 

=> 

Input  data  deck  to  ION 

3.  3.  1.  3 

The  alternate  version  of  the  ionospheric  program  consists  of  two 
separate  entities,  the  preprocessor  TARGETS'  and  the  reduction  program 
ION1,  External  data  transfer  in  and  out  of  TABGEN  is  handled  through 
four  files:  an  Input  data  card  deck,  data  output  to  line  printer,  the  inovit 
ionospheric  coefficient  tape,  and  the  output  disc  or  tapo  file  with  f0F2-hr 
tables.  External  data  transfer  in  and  out  of  ION1  is  handled  through  the 
following  3  data  files;  nr,  input  uuta  card  deck,  the  input  data  file  with 
i0I  2-h.  tables,  and  output  to  the  line  printer. 


P  roc  ra  m 

File  Type 

Mode 

I/O 

Fortran 

Unit  No, 

Dear  riptior, 

Dcta  il  s 

Unde  r 

TABGEN 

Tape 

R’N 

1 

1 

Lm,  toeff,  tape 

3.  3.  1 ,  1 

TABGEN 

Disc  or  tape 

B  IN 

O 

2 

File  with  fnF2 -httable s 

3.  3.  1,  4 

TABGEN 

Line  printo r 

BCD 

o 

b 

Output  listing  from 
TABGEN 

3.  3.  1.  h 

TABGEN 

C a rd  roade r 

BCD 

I 

6 

Input  data  deck  to 
TABGEN 

3.  >,  1,6 

ION  l 

Di  s  t  o  r  t ape 

BIN 

I 

2 

File  with  f0F2-hrtables 

3.  3,  1.  4 

ION  1 

Card  reader 

BCD 

I 

S 

Input  data  deck  to  ION1 

oc 

I  ON  1 

Line  printer 

BCD 

o 

6 

Output  listing  from 

ION  1 

3.  3.  1.  7 

3.  3,  1,  1  Ionospheric  Coefficient  Tape 


There  are  36  fixed  length  records  on  the  cape  followed  by  a  double 
end-of-file.  Each  record  contains,  in  3848  words,  the  generalized  10-day 
f0F2  and  30-day  M(3000)F2  coefficients  to  be  used  for  one  third  of  one  month. 
The  36  records  are  in  time  sequence  and  valid  for  the  periods  January  1-  10, 
January  11-20,  January  21-31,  February  1-10,  February  11-20,  February  21- 


28  or 

29,  . 

.  .  .  ,  Dec* 

embe  r  21-31. 

Word 

Mode 

Fortran  Name 

Description 

1 

Integer 

LOND 

=  (month*  1  00+day  ),  fir  st  date  for  which 
coefficients  are  valid 

2 

I.itege  r 

LONDY 

=  (month<  1  00+day  ),  last  date  for  which 
coefficients  are  valid 

3-2966 

Real 

WCOEF 

Array  of  dimension  3*1  3*76  of  general¬ 
ized  f0F2  coefficients  valid  for  the  time 
interval  specified  by  words  1  and  2 

2967. 

3407 

Real 

UM 

Array  of  dimension  9*49  of  M(3  000)F2 
coefficients  valid  for  a  12 -month  running 
average  of  the  sunspot  number  =  0,  and 
to  be  used  for  the  time  interval  specified 
by  words  1  and  2 

3408- 

3848 

Ileal 

UM  1 

Array  of  dimension  9‘49  of  M(3000)F2 
coefficients  valid  for  a  12-month  running 
average  of  the  sunspot  number  100,  an: 
to  be  used  for  the  time  interval  specified 
by  word'1  1  and  2 

The  formation  of  the  specific  coefficient  setB  for  f,,F2  and  M(30Q0)F2  from 
the  gene  ra  1  i  oefficients  is  discussed  unue  r  3,1,  1,  1,  C  PC  N u,  2, 


3.  3.  1.  2  Line  Printer  Output  Listing  from  ION 


The  typical  output  format  of  the  results  from  ION  is  shown  for  some 
test  cases  unde  r  4.  1,  In  addition,  the  following  otror  messages  may  oi  cur; 

Printed  in  PROGRAM  ION,  "Error  in  solar  input  data  lor  year  ,  .  and 
month  --  where  upon  the  computer  run  is  terminated. 
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Printed  in  PROGRAM  ION,  "Remaining  update  data  not  used";  if  more 
than  eight  update  conditions  are  supplied,  the  first  eight  are  used,  the  re¬ 
maining  card*  are  skipped  over, 

Printed  in  SUBROUTINE  P.EFRAC,  "Coefficients  not  found  on  tape  for 

year,  month,  day  = . ",  where  upon  control  is  transferred  to 

PROGRAM  ION  to  proceed  with  the  next  data  case. 

Printed  in  SUBROUTINE  BETA,  "Ray  is  reflected  at  ionosphere  or  near 
reflection  condition,  elevation  angle  correction  is  not  computed,"  where  upon 
control  is  transferred  to  SUBROUTINE  REFRAC  to  proceed  normal  wl'h  the 
remaining  computations. 


3,  3,  1.  3  Input  Data  Card  Deck  to  ION 

The  input  card  deck  to  ION  specifies  the  output  and  update  options  and 
it  defines  the  evaluation  and  update  conditions  and  tho  required  solar  data. 
The  set  up  procedure  for  the  card  deck  is  described  below  followed  by  a 
description  of  the  solar  data  and  by  the  detailed  card  type  and  format  in¬ 
formation. 


a)  Procedure  to  Set  Up  Card  Deck  for  ION 


Specify  options  ** 

Card  type  1  ;  ISEL(l)  -  ISEL(5),  output  options  for  ionospheric  profile 

and  refraction  corrections.  =0  wanted,  =1  not  wanted. 


Card  type  2  •  IUPDT,  IDRDAV,  update  option  and  output  option  for  correc¬ 

tion  to  range  differencing.  =0  not  wanted,  =1  wanted.  If 
=  1,  additional  input  data  is  required,  cards  9  and  10  and/or 
card  11, 


+#  Specify  evaluation  condition 

Card  type  3  :  FS,  FLAT,  FLON,  station  information;  wave  frequency, 
latitude  and  longitude.  If  refraction  corrections  are  not 
desired,  FS  is  not  used  and  should  be  left  blank  or  set 
=  0  or  positive. 


Card  type  4  ;  ELEV,  AZ,  HS,  EDOT,  HDOT,  satellite  information;  elevation 

angle,  azimuth,  height,  elevation  rate,  altitude  rate.  If 
the  instantaneous  range  rate  correction  is  not  desired,  EDOT 
and  HDOT  are  not  used  and  should  bo  left  blank  or  set  to  any 
value. 


87 


Card  type  5  ;  IYR,  MON,  IDAY ,  TIME,  time  information:  year,  month, 

day,  time. 

*+  Specify  solar  data  ** 

#*  If  the  year  and  month  of  thia  condition  are  the  same  as  the  year  and  month 
of  the  previous  condition,  skip  cards  6,  7,  8. 


Card  type  6 


IYR,  MON,  FLX(1)-FLX(16),  date  and  daily  values  of  observed, 
solar  flux  for  the  first  16  days  of  the  month.  If  future  pre¬ 
dictions  are  to  be  evaluated,  leave  array  F LX  blank. 


Card  type  7 


IYR,  MON,  FLX(17)-FLX(31),  date  and  daily  values  of  ob¬ 
served  solar  flux  for  the  latter  part  of  the  month.  If  there 
are  less  than  31  days  to  the  month,  the  additional  spaces 
are  normally  left  blank. 

Card  type  8  •  IYR,  MON,  SIS,  SIF,  date  and  12-month  running  average  of 
sunfipot  number  and  solar  flux, 

Preparation  of  solar  data  is  discussed  under  b). 

**  Specify  update  data  ** 

++  If  update  is  not  desirod,  IUPDT  =  0  on  card  2,  skip  cards  9  and  10, 

Card  type  9  :  NUPDT,  number  of  observation  conditions  to  be  used  for 

updating  the  predictions  for  the  evaluation  condition, 
Maximum  =  8, 

**  If  update  is  not  desired  for  this  particular  evaluation  condition,  NUPDT'O, 
skip  card  1  0. 

Card  type  10:  ULAT,  ULON,  ULEV,  UZIM,  UT,  OBS,  ITYPE,  update  data: 

latitude,  longitude  of  observation  station,  elevation  and  azi¬ 
muth  of  observation,  observation  time,  value  of  measure¬ 
ment  and  type,  When  the  observation  is  critical  frequency 
set  elevation  to  90°  and  azimuth  to  0°.  For  vertical  and 
angular  content  use  the  appropriate  angles, 

**  Rcpc&t  card  10  until  all  NUPDT  conditions  are  defined. 

**  Specify  additional  data  for  range  differencing 

+*  If  corrections  to  range  rate  by  differencing  technique  are  not  desired, 
IDRDAV  *  0,  skip  card  11. 

Card  type  11  :  E  LEV,  AZ  ,  HS,  TIME,  satellite  information,  elevation, 

azimuth,  and  height  and  time  information  for  the  tecond 
observation  used  for  the  range  differencing. 

+  *  Repeat  cards  3  through  11  for  any  number  of  conditions  desired. 

+*  Terminate  vith  card  3  containing  a  negative  value  for  the  wave  frequency 
FS. 
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b)  Preparation  of  Solar  Data 

The  solar  data  can  be  extracted  from  the  "Solar-Geophysical  Data" 
monthly  publications,  issued  by  NOAA,  Boulder,  Colorado. 

The  daily  values  of  solar  flux  are  to  ba  copied  from  the  table  "Daily 
Solar  Indices"  (normally  page  7)  under  the  column  "Observed  Flux  Ottowa 
2800"  MHz  (corresponds  to  10.7  cm  wavelength ),  If  future  predictions  are 
to  be  evaluated  and  therefore  no  measurements  available,  the  daily  flux 
values  required  on  card  with  the  appropriate  year  and  month  are  to  be  loft 
blank,  The  program  automatically  checks  for  this  condition  and  inserts  the 
beat  eotimate  for  the  daily  flux  values  which  is  the  12-month  running  average 
of  the  solar  flux. 


The  12-month  running  average  ll8,  ,  for  month  j  of  a  solar  index  1  with 
a  mean  value  Tk  for  month  k  is  defined  as, 


I. 


(l 

12 


+5 


.  -eJ-tie*  ♦  > 
V  2 


lJ+! 


1--  6 


The  monthly  means  of  t.he  index  for  the  month  under  consideration,  for  1  through 
5  month  past  and  prior  and  half  the  value  of  the  monthly  mean  for  6  months 
past  and  prior  are  added  and  divided  by  12,  yielding  an  average  over  12  months 
centered  around  the  specified  month.  The  12-month  "unning  average 
( =amoothed )  of  the  sunspot  numbe  rs  Sx  2  ,  .  for  month  j  are  listed  in  the  "Solar 
Geophysical  Data"  publication  (normally  page  9)  in  table  "Smoothed  Observed 
and  Predicted  Sunspot  Numbers"  and  are  to  be  used  for  paBt  as  well  as  future 
evaluations , 

The  12-month  running  average  of  the  solar  flux  is  computed  from  the 
accumulated  monthly  means  using  the  formula  above,  The  monthly  means  arc 
listed  along  with  the  daily  values  of  solar  flux.  If  not  enough  advance  data  is 
available  to  form  the  12-month  running  average,  that  value  can  be  approxi¬ 
mated  with  a  11.5,  10,  5  or  9.  5-month  running  average; 
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(—  6-* 

1=  -5 

If  not  even  onough  data  ia  available  to  form  a  9.  5-month  running  average, 
an  eatimate  of  the  12-month  running  average  of  the  a.olar  flux  can  be  derived 
from  the  12-month  running  average  of  the  aunapot  number  for  which  tabulated 
predictiona  are  available.  The  relationahip  between  aolar  flux  and  aunapot 
number  was  arrived  at  by  Stewart  and  Leftin,  Reference  9. 

approx.  Fla  ,  3  =  63.  75  +  0.  728  *S;a,  3  +  0.  000&9  , } 

The  attached  tablea  contain  the  final  12-month  running  averages  for 
aunapot  number  and  aolar  flux  from  I960  on  and  the  monthly  me  ana  for  solar 
flux  from  1970  on. 


J+ii 


,  k  =  1 , 2,  or  3 
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Table  1.  12-Month  Running  Average  of  the  Zurich  Relative  Sunipot  Number 


I960 

1961 

1962 

1963 

1964 

1965 

1966 

1967 

1968 

1969 

Jan. 

128.  9 

80.  2 

45.  2 

29.4 

19.  5 

11.7 

27.  7 

75.  0 

102.  6 

no.  o 

Feb. 

125.  0 

74.  8 

41.  8 

29.  8 

17.  8 

12.  0 

31.3 

78.  8 

102.  9 

109.  6 

Mar. 

121.  6 

68.  8 

39.  8 

29.  7 

15.  4 

12.  5 

34.  5 

82.  2 

104.  7 

108.  0 

Apr. 

119.  6 

64.  3 

39.  4 

29.  0 

12.7 

13.  6 

37.4 

84.  6 

107.  2 

106.  4 

May 

117.  0 

60.  1 

39.  2 

28.  7 

10.  8 

14.  6 

40.  7 

87.  5 

107.  6 

106.  2 

June 

113.  9 

55.  8 

38.  3 

28.  2 

10.  2 

15.  0 

44.  7 

91.  3 

106.6 

106.  1 

July 

108.  6 

53.  1 

36.  8 

27.  7 

10.  3 

15.  6 

50.  3 

94.  1 

105.  2 

105.  8 

Aug. 

102.  4 

52.  5 

34.  9 

27.  2 

10.  2 

16.  4 

56.6 

95.  3 

104.  8 

106.4 

Sept. 

97.  9 

52.  3 

32.  7 

26.9 

9.  9 

17.  4 

63.  1 

95.  3 

107.  0 

105.  4 

Oct. 

93.  3 

51.4 

30.  8 

26.  0 

9.6 

19.  7 

67,6 

95.  0 

109.9 

104.  1 

Nov. 

87.  9 

50.  5 

30.  0 

23.  8 

10.  1 

22.  3 

70.  2 

97.  1 

no.  6 

104.  6 

Dec. 

83.  7 

48.  7 

29.  8 

21.3 

11.  0 

24.  5 

72.7 

100.  6 

110.  1 

104.  9 

1970 

1971 

1972 

Jan. 

105. 6 

80.  4 

70.  8 

Fob. 

106.  0 

77.  8 

71.  2 

Mar. 

106.  2 

74.  4 

72.  4 

Apr. 

106.  1 

7  0.  9 

73.  4 

May 

105.  8 

68,  1 

72.  9 

J  une 

105.  3 

66.  7 

70.  5 

July 

103.  8 

65.  5 

68.  1 

Aug, 

101.0 

65.  0 

65.4 

Sept. 

97.  2 

66.4 

62.  0 

Oct. 

93.  9 

o7 .  1 

60.  3 

Nov. 

89.  4 

67.6 

3  S.  5 

Dec. 

84.  1 

69.  9 

54.  8 

'll 


Table  2.  12-Month  Running  Average  of  the  Solar  Flux  at  10.  7  cm  Wavelength  (Ottawa) 


1960 

1961 

1962 

1963 

1964 

1965 

1966 

1967 

1968 

1969 

Jan, 

178.  7 

128.  9 

97.9 

81.8 

76.4 

73.  8 

85.  7 

128.  1 

150.  0 

150.  2 

Feb. 

174.  6 

124.  1 

95.2 

81.8 

75.  5 

74.  4 

88.  4 

131.  5 

149. 4 

150.  2 

War. 

170.  8 

119.  1 

93.  1 

81.  8 

74.  4 

74.  9 

91.2 

134.  3 

149.  3 

15  0.  1 

Apr. 

168.6 

115.  1 

91.7 

81.  5 

73.  3 

75.4 

93.  8 

136.  3 

150.  4 

150.  0 

May 

166.  2 

no.  8 

91.  1 

81.  2 

72.  5 

75.  8 

96.  5 

138.  8 

150.  8 

150.  8 

June 

162.  9 

106.  6 

90.4 

81.  0 

72.  2 

76.  0 

100.  1 

141,  7 

149.  9 

151.4 

July 

157.  8 

103.  7 

89.  2 

80.  6 

72.  3 

76.  4 

104.  6 

145.  0 

147.  8 

151.4 

Aug. 

151.  8 

102.  4 

87.7 

80.  3 

72.4 

77.  2 

109.  7 

147.  8 

145,  5 

152.  6 

Sept. 

147.  4 

102.  0 

85.  8 

80.  1 

72.  2 

78.  3 

115.  3 

148.  2 

146.  0 

152.  8 

Oct. 

143.  1 

101.  5 

84.  2 

79.  8 

72.  1 

80.  0 

119.6 

147.  4 

148.  3 

152,5 

Nov. 

137.  9 

101.  1 

83.  1 

78.  7 

72.  5 

81.  9 

122.  8 

147.  9 

149.  0 

153.  7 

Dec, 

133.  1 

100.  2 

82.  3 

77.  3 

73.  2 

83.  6 

125.  7 

149.  3 

149. 4 

154.  4 

1970 

1971 

1972 

J  an. 

154.  7 

135.  0 

120.  5 

Feb. 

155.  1 

132.  5 

121.2 

Mar. 

155.  2 

129.  9 

122.  1 

Apr . 

155.  2 

126.  6 

123.  1 

May 

155.  2 

122.  8 

123,  2 

June 

155.  8 

1  19.  7 

121.  7 

July 

136.  3 

1  16.  5 

120.  3 

Aug. 

155.  0 

1  14.  7 

118.  0 

Sept. 

151.4 

1  15.  5 

115.0 

Oct. 

147.  6 

1  16.  1 

113.  5 

Nov . 

143.  3 

1  16.  7 

111.8 

Dec. 

138,  6 

1  18.  9 

108.  6 
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Table  3. 

Monthly  Mean  of  the  Solar 

Flux 

1970 

19  n 

1972 

1973 

January 

158.  3 

162.6 

114.  8 

102.  2 

February 

175.  4 

13;.  8 

141.  8 

98.  7 

March 

158.  4 

111.9 

128.  5 

100.  4 

April 

162.  0 

1 16,7 

112.  9 

105.  0 

May 

168.  4 

109.  9 

129.6 

97.  0 

June 

154.  9 

101.  7 

135.  4 

91.2 

July 

152.  0 

117.  4 

122.  0 

August 

138.  2 

114.  1 

125.  7 

Septembe  r 

143.  2 

104.  0 

113.  7 

Octobe  r 

148.  3 

107.  2 

121.  1 

Novembe  r 

162.  0 

114.  0 

101.6 

Decembe  r 

152.  8 

124.  5 

102.  9 
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Card  Type  4 

Satellite  data  for  evaluation  condition 


10 


l  I  I  I  I  I  I  I  I  |  I  I  l  I  I 
— E  LEV—'1  l  AZ 


20  30  40  30  60 

!  i  ;  I  i  •  i  1 !  ■  •  .  1 1  ! 

•  -  HS  — '  -  EDOT  _J  -  HDOT  


80 


Format  (2  F10.  6,  F10.0,  2D1S.8) 


W  ord 
No. 

Prog  ram 
Variable 

Units 

Format 

Column 

Description 

1 

ELEV 

degreei 

F1C.  6 

1-10 

Elevation  angle  to  satellite 

2 

AZ 

degreei 

F  1  0,  6 

1  1-20 

Azimuth  angle 

3 

HS 

km 

F10.  0 

21-30 

Height  of  satellite  above  surface  of  earth 

4 

EDOT 

rad / se  < 

D  1  5.  S 

3  1-45 

Elevation  rate 

5 

HDOT 

m/lec 

D 1  5.  8 

46-60 

Altitude  rate 
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Cafd  Type  5 

Time  data  for  evaluation  condition 


Format  (315,  110.7) 


Word 

No. 

Program 

Variable 

JnitB 

F  ormat 

Column 

Desc  ription 

1 

IYR 

1  5 

1-5 

Year  (last  2  digits ) 

2 

MON 

I  5 

6-  10 

Month  (-1  through  12) 

3 

IDAY 

I  5 

11-15 

Day  (  =  1  th rough  3  1) 

4 

TIME 

\o\x  r  8 

F10.  7 

16-25 

Universal  time 

<18 


IYR 

MON 


Year  (last  2  digits  ) 

Mo  nth  (-1  through  12) 

Daily  solar  flux  xlC  for  aay  1  of  month 
Daily  solar  flux  xlG  for  day  2  of  month 


Daily  solar  flux  xlO  for  dav  lb  of  the  month 


Year  (last  Z  digits  ) 

Month  (=  1  through  12) 

Daily  solar  flux  x  1C  for  day  17  of  jnth 
Daily  solar  flux  x  1C  for  uay  18  of  month 

Daily  solar  flux  xlC  for  day  31  of  month; 
if  the  month  na  b  less  than  3L  Jays,  tne 
•  pare  locations  are  left  blank 


Month  (-  i  through  12) 

12-month  running  e  veragi  of  sunspot 
numbs r  x  1 0 

12-month  running  ei/erage  of  solar  flux  x  10 


Card  Typa  9 


Update  control  conitant 


10 


to 


>0 


1 1 1 1 1  1 1 1 1 1 1 1 1 1 1 


a 

a, 

a 

2 


Format  (I  5 ) 


1  i  i  l  ;  l  ■  :  i  I  i  ' 


40 


BO 

f 

t  ;  '  I 


60 


70 


to 


W  ord 
No. 


1 


Program 

Variable 


NUPDT 


Unit* 


Format! 


I  5 


Column  ! 


Description 


1-5  I  Number  of  update  conditions,  maximum 


i 


Update  data  condition 


/' 


10 


20 


so 


40 


60 


/ 

1  i  1  1 1  1  1  i  1 

1  1  I  I  1  1  :  '  |  i  .  I  ,  :  !  j  '  1  .  : 

-  1 

! 

1 — ULAT-- 

— ULON-J  —  ULEV-  L-UZIM  -  1 

i 

L-  U’T  •-> 

eo 


!  .  ,  I 


OBS 


70 


80 


a 

> 

h 


rormat  (2F10.  5,  2F10.6,  F10.7,  D1S.B,  I  5) 


W  o  r  d 

No. 

Prog  ram 
Variable 

U  nita 

Format 

Column 

De»  criptlon 

1 

U'LAT 

degree  I 

F  1  C ,  5 

1-  10 

Latitude  of  update  itatlon 

2 

ULON 

degree! 

F  1  0 ,  5 

11-20 

Longitude  of  update  nation  Ipontiva  east, 

0- 36C  degree*  ) 

3 

ULEV 

degree i 

FiO.  6 

21-30 

Elevation  angle  of  observation  (=90  for  f0F2 
data,  *  elevation  to  satellite  for  vertical 
and  angular  electron  content) 

4 

UZIM 

degreei 

F  10,  b 

3  1-40 

Azimuth  angle  of  observation 

5 

UT 

hour* 

F10.  7 

41-50 

Univer»al  time  of  ob*ervation 

6 

OBS 

MHs  or 

e  /ma 

D  1  5 .  8 

5  1-65 

Observation  to  be  u*ed  for  update 

7 

IT  Y  PE 

I  5 

66-70 

Obiervation  flag,  *  1  for  f0F2,  -  2  for 

vertical  electron  content,  «  3  angular  elactton 
ccntant 


Card  Type  1 1 

Satellite  and  time  information  for  lacon-  ob»ervation-u*ed  for  range  differencing 


Format  (2F10.  6,  F10.  0,  F10.  7) 


Word 

No, 

Prog  ram 
Variable 

Unit* 

Format 

Column 

De«c  riptio.n 

i 

ELEV 

degree* 

F  1 0,  6 

1.10 

Elevation  angle  to  eatellite 

2 

AZ 

degree* 

F  1 0,  6 

11-20 

Azimuth  ar.ule 

3 

HS 

F  1  C.  C 

2  1-30 

Height  of  nateihte  above  »uriace  of  earth 

4 

TIME 

hour* 

F  1  0.  7 

3  1.40 

Univer»al  time 

Bt' 


3,  3,  1,  4  lonoapheric  Data  File  with  f„F2-h,  Tablet 

The  file  with  £0F2-h,  tables  i*  generated  in  PROGRAM  TABGEN  for 
u*e  in  the  alternate  ionospheric  veriion  PROGRAM  ION1,  It  conaiiti  of 
fixed  length  record*,  a*  many  a*  wore  generated  in  PROGRAM  TABGEN, 
terminated  by  a  lingle  aid. of. file,  Each  record  contain*,  in  354  word*, 
the  date,  the  (tation  position,  the  daily  aolar  flux  and  value*  of  critical 
frequency  f0F2  and  corresponding  height  h,  .  The  value*  for  f0F2  and 
h,  are  tabulated  for  the  given  date  for  14  different. time*  at  each  location  of 
a  25  point  pattern  around  the  »tatior,  which  cover*  the  ionoephere  viiible 
from  that  station, 


Word 

Mode 

Fortran  Nam* 

Description 

1 

Integer 

IYMD 

Datei  year*  10000  t  month  *  1 00  +■  day 

2 

Real 

FLAT 

Latitude  of  »t*tion  in  radian* 

3 

Real 

FLON 

Longitude  of  station  in  radian* 

4 

Real 

FLUX 

Value  of  daily  *olar  flux  (if  the  daily 
flux  it  greater  than  13  0,  the  limit 
value  of  130  i»  substitutad ) 

5-354 

Integer 

IFH 

Array  of  dimon»ion  14<25  containing 
packed  integer  tabulated  value*  for 
f0F2  and  h,  for  14  local  time  hour* 
at  each  location  of  the  25  point  pattern 
around  the  station,  Each  integer  ha* 

8  digit*,  the  first  4  digits  define  h,  in 
units  of  km,  the  last  4  digits  give 

f0F2  in  units  of  ^  MHz, 

3,  3,  1,  5  Output  Lifting  from  TABGEN 

The  only  line  printer  output  from  TABGEN  i*  the  printout  of  the  input 
data  condition*,  In  addition,  the  following  error  messages  may  occurt 

Printed  in  PROGRAM  TABGEN,  Error  in  *olar  input  data  for  year  =  , 
and  month  ■  .  ,  where  upon  the  computer  run  i*  terminated, 
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Printed  in  FR0GRAM  TABGEN,  "Coefficient*  not: found  on  tape  for 
year,  month,  day  ■  .  ,  where  upon  the  comparer  run-ia  terminat#dr 

3,3.  l,  6  Input  Data  Deck  to  TABGEN 

The  input  to  TABGEN  conaiat*  of  card  type  12,  shown  on  the  next  page 
containing  date  and  itation  Information  and  of  card  type*  6,  7,  8  ai  described 
under  3,  3.1. 3  c)  apecifying  the  aolar  data, 


Card  type  12  i  IYR,  MON,  IDAY,  FLAT,  FLON,  ysar,  month,  day, 
latitude  and  longitude. 

**  If  the  yea;'  and  month  of  this  condition  are  the  lame  as  the  year  and 
month  of  the  previous  condition,  iktp  cards  6,  7,8. 

Card  type  6  !  IYR,  MON,  FLX(  1  )-FLX(  16  ),  date  and  daily  value*  of 

obierved  lolar  flux  for  the  firit  16  day*  of  tho  month, 

If  ‘uture  prediction!  are  to  be  evaluated,  leave  array  FLX 
blank, 


Card  type  7  ,  IYR,  MON,  FLX(  1 7  ). FLX(  3 1 ),  date  and  daily  value,  of 
obierved  aolar  flux  for  the  latter  part  of  the  month,  If 
there  are  leai  than  31  daya  to  the  month,  the  additional 
opacei  are  normally  left  blank, 


Card  typo  8 


IYR,  MON,  SIS,  SIF,  date  and  12-month  running  average 
of  lunipot  number  and  lolar  flux, 


Preparation  of  the  lolar  data  ia  diicueied  under  3.  3,  1,  3  b), 


**  Repeat  card*  12,  fe,  7,  8  for  any  number  of  condition*  deiired, 

**  T B rmlnatn  with  card  12  containing  a  aero  or  negative  value  for  the 
year  IYR, 


3,  3,  1,  7  Output  Llating  from  10N1 

The  typical  output  format  of  thi  reaulta  from  ION1  ia  ihown  for  lome 
teat  caaea  urder  Saction  4,  1,  In  addition,  the  following  error  meaiagei 
may  occuri 
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t; 


I 

i 

j 


Date  and  itation  evaluation  condition 


I  I  II  I  I  I  I  l 


10 


JO  50  45 

I  i  .  I  I  i  ,  ■  .  i  ■ 


at 

>> 


Z  > 

5  q  flat 


FLON 


Format  (  315,  2F10.  5) 


Word 

No, 

Program 

Variable 

Unit* 

F  ormat 

Column 

Oeacription 

1 

1YR 

-  - 

I  5 

1-5 

Year  (lut  2  digit! ) 

2 

MON 

I  5 

6-10 

Month  (  =  1  through  12) 

3 

IDAY 

1  5 

11-15 

Day  (  ■  1  through  3 1 ) 

4 

FLAT 

degree! 

F  1 0.  5 

16-25 

Station  latitude 

5 

FLON 

degree! 

F  1  0,  5 

26-35 

Station  longitude  (poiitive  eaat,  0-360  degree* 
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Printed  In  SUBROUTINE  RJEFRC  1,  "f0FZ.h,  tahto*  for  thi»  station 
and  date  not  found  in  file,"  where  upon  control  i»  traowferred  to  PROGRAM 
ION  1  to  proceed  with  the  next  data  caae, 

Printed  in  SUBROUTINE  BETA,  "Ray  i»  reflected  at  ionosphere  or  near 
reflection  condition,  elevation  angle  correction  i*  not  computed,"  where 
upon  control  la  traneferred  to  SUBROUTINE  REFRC  k  to  proceed  normal  with 
the  remaining  computation*. 


3.  3,  1,  8  Input  Data  Deck  to  ION1 

The  input  data  to  ION1  involve*  only  card  type*  3,4,  and  5  as  they  are 
described  under  3.  3.  1.  3  c )  to  epecify  the  evaluation  condition, 


Card  type  3  ;  FS,  FLAT,  FLON,  atation  informations  wave  frequency, 

latitude  and  longitude.  Set  FS  =  Q  or  positive,  if  refraction 
correction*  are  not  requested. 


Card  type  4  •  ELEV,  AZ,  HS,  EDOT,  HDOT,  satellite  information: 

elevation  angle,  azimuth,  height,  elevation  rate,  altitude 
rate,  If  the  instantaneous  range  rate  correction  is  not 
desired,  EDOT  and  HDOT  are  not  used  and  can  be  left 
blank  or  set  to  any  value, 


Card  type  5 


IYR,  MON,  IDAY,  TIME,  ume  information;  year,  month, 
day,  time. 


♦  *  Repeat  cards  3  through  5  for  any  number  of  condition*  desired, 


*#  Terminate  with  card  3  containing  a  negative  value  for  the  wave  frequency 
FS, 
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3,4  Computer  Program  Functional  Flow  Diagram 


The  functional  flow  diagram  of  the  Bent  Ionoepheric  Program  ION  i»  pre 
aenteu  at  well  at  the  diagram*  for  the  alternate  version  TABGEN-ION 1 , 

The  label*  to  the  right  top  of  each  block  *pecify  the  prcgram/*ubroutine* 
that  perform  the  function  described  in  the  block,  Lower  level  flowchart* 
disclosing  more  detail*  are  listed  under  the  individual  computer  program 
component  description*  in  Section  3,2,  1,2, 


_ V  ION 

Retd  ind  prepare 
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L  solar  data  as  l 
needed  from  card* 
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PROFL1 


Compute  critical 
frequency  and 
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haijjht 


Compute  prediction* 
matching  update 
obeervatioai 
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I 

Compute  final 
update  ratio  and 
update  critical 
frequency 


PRCFL2 


|  Compute  remain-  | 
ing  lonoaphenc 
profile  parameterij 


llevatiorTs^^  \o 
^-kngle  correction 
desired 
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elevation  angle 
refraction 
correction 


Update 
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Rnad  aolar  data 
from  carda 
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lit  sunriae 
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r—  ——————————— 
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f0F2-hm  table,  date 
atation  and  flux  to 
file 


More  N 
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4,  Q  Quality  Assurance  _ 

All  sspects  jrf  the  ionospheric  model  were  tested  thoroughly  during  end 
after  the  development  phaee  and  some  of  the  result* are  shown  in  Section  6,2, 
The  shape  of  the  electron  density  versus  height  prof®*  was  compared  with 
actual  composite  profiles  compiled  at  NASA/GSFC  ud  they  were  always  in 
close  agreement,  The  integrated  electron  content  was  compared  extensively 
with  the  vert  .cal  electron  content  derived  from  Faraday  rotation  measure  - 
ments.  The  results  of  thl«  work  performed  for  SAhfSO,  are  described 
in  Reference  2.  The  predictions  alone  accounted  for  70  to  SO %  of  the  actual 
electron  content  and  after  updating  with  ionospheric  obae rvationa ,  up  to 
90%  of  the  ionosphere  was  estimated.  The  ionospheric  refraction  corrections 
were  tested  in  orbit  determination  work  performed  at  NASA/GSFC,  The 
iterative  leant  square  reduction  programs  were  run  with  and  without  iono¬ 
spheric  corrections  and  the  final  RMS  values  of  the  measurement  residuals 
were  greatly  reduced  by  30  to  75%  upon  use  of  ionospheric  corrections, 

After  modifying  the  ionospheric  program  to  its  current  form,  a  number 
of  test  cases  were  run  and  the  results  including  all  possible  outputs  were 
compared  with  results  from  previous  runs  before  modifications ,  The  same 
test  cases  lilted  under  4,  1  should  be  checked  out  whenever  the  program  is 
duplicated  and  transferred  to  another  computer  system  to  insure  thar  all 
parts  of  the  program  are  in  working  order. 

4,  1  Test  Plan/Procedure 

The  following  page*  show  s  list  of  the  input  card  deck  and  the  corres¬ 
ponding  printed  output:  results  for  test  cases  L  through  5  and  a  cross  reference 
list  in  Table  4  of  the  various  conditions  tssted.  The  five  test  cases  evaluate 
the  function!  of  the  ionospheric  program  for  various  possibilities  in  latitude, 
longitude,  local  time,  season,  and  solar  activity  effecting  the  ionospheric 
profile  and  therefore  also  electron  content  and  refraction  corrections.  Each 
of  the  five  tuwt  cases  computes  all  possible  output  results!  critical  frequency 
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And  corresponding  height,  the  values  of  half  thickness  and  the  decay  constants 
for  the  shape  of  the  profile,  the  profile  plot  and  list,  vertical  and  angular 
electron  content  and  refraction  corrections  to  elevation  angle,  to  range, 
to  instantaneous  range  rate  and  to  range  differencing. 

For  the  standard  ionospheric  PROGRAM  ION  the  input  is  listed  in 
Table  5  for  all  five  test  cases,  and  the  ouput  in  Table  6,  For  the  alternate 
version  of  the  ionospheric  program,  the  input  and  output  of  the  preprocessor 
PROGRAM  TABGEN  are  shown  in  Table  7  and  the  input  and  output  of  the 
reduction  PROGRAM  ION1  are  giver,  in  Tables  8  and  9  respectively.  Only 
test  cases  1,  2,  and  5  are  presented  for  the  alternate  program  since  the 
update  capability  tested  in  cases  3  and  4  is  not  included  in  this  version, 

4,2  Other  Quality  Assurance  Provisions 

Whenever  the  program  is  reproduced  for  use  on  another  system,  the 
program  card  decks  should  be  duplicated  and  verified,  If  the  program  is 
transferred  to  a  system  with  compatible  binary  coding,  the  binary  magnetic 
tape  containing  the  ionospheric  coefficients  should  be  copied  and  verified, 

If  the  program  is  to  be  used  on  a  computer  with  different  binary  word  or 
record  structure,  tne  binary  i ape . should  be  copied  to  a  BCD  tape  and  at  the 
new  location,  transferred  back  onto  a  binary  tape,  Care  should  be  taken  that 
during  the  binary  to  BCD  tape  copy  process  no  loss  of  significance  will  occur, 
which  means  the  format  (E  17.  1 1 }  is  requi red  for  the  general  f0F2  coefficients 
and  the  format  (E14,  8)  is  required  for  the  general  M(3000)F2  coefficients. 

The  binary  tape  format  is  described  under  3.  3,  1.  When  tape  and  card  deck* 
arc  available  on  the  new  system,  the  test  runs  described  in  4,  1  should  be 
performed  and  the  results  compared  with  the  output  results  in  the  tables  for 
ag  reement, 
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I-NVT  ••  SECOND  SATELLITE  POSITION  USED  EON  RAHtE  DIFFERENCING 

FLEVATION-  80*000000  DEG*  AZTrtUTM.3*  1.000000  KCi  HEIGHT-  P0O7.0  RH*  ll*TI  HE»15 .002*770  HRS 
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range  RATE  CORRECTION  FOR  RANGE  DIFFERENCING  OVER  10*0001  SECONDS  »  **2D30AAE»01  H/SEC 
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5,  0  Preparation  for  Delivery 

The  completed  CPC-1  (Computer  Product  Configuration  Item)  for  the 
Bent  Ionoepherlc  Model  coniiiti  of  three -part*  which, ere  pecked  and  ehippad 
•  eparatelyi  a  magnetic  tape,  card  deck*,  acd  k  documentation  manual.  The 
tape  ii  mailed  first  clan  or  airmail  and  ie  marked  with  ''Special  Handling- 
Electro  Magnetic  Item,*'  The  card  deck*  and  manuals  can  be  shipped  third 
clan,  For  etorage  of  the  magnetic  tape  and  thw  card  decki,  a  cool  and  dry 
place  ehov.ld  be  lelected  to  insure  that  the  good  condition  of  the  item*  ie 
preserved.  The  following  lilt  described  the  delivered  items  in  detail! 

s)  Magnetic  tape  containing  the  general  ionoipheric  coefficients  in 
either  BCD  or  Binary  code  depending  un  the  compatibility  of  the 
computeri  between  which  the  tremfer  occur*. 

b)  Card  deckn 

1 )  Fortran  card  deck  to  copy  the  BCD  tape  with  ionoepheric 
coefficient!  to  a  Binary  tape  of  the  proper  form,  Thii  deck 
ia  not  needed  if  the  required  Binary  tape  ii  supplied  in  place 
of  the  BCD  tape, 

2)  Fortran  card  deck  for  PROGRAM  ION,  itandard  veraion  of 
the  ionoipheric  program, 

3)  Data  cards  for  teitrun  of  PROGRAM  ION. 

4)  Fortran  ca^d  deck  for  PROGRAM  TABGEN,  preprocenor 
for  the  alternate  veriion  of  the  ionoipheric  program, 

5)  Data  cards  for  teitrun  of  PROGRAM  TABGEN, 

6)  Fortran  card  deck  for  PROGRAM  ION1,  reduction  program 
for  the  alternate  veriion  of  the  ionoipheric  program, 

7)  Data  cardi  for  testrun  of  PROGRAM  ION1, 

8)  Additional  data  cardi  of  iolar  input  data  from  1962  to  1973, 

c)  Manuali  "Documentation  and  Description  of  the  Bent  Ionoipheric 
Model,  1  For  the  letup  and  checkout  of  the  programi,  Section  3,  3,  1 
File  Description  and  Section  4,  1  T*»t  Plan  ihould  be  consulted. 
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6,  0  Notes 

Section  6,  1  describes ..th*  development  of  the  ionospheric  model, 
the  data  base  on  which  the  was  found# d  *nd  th»  justification# 

for  th#  derivation  of  each  etep  in  tha  davelopment,  In  Section  6.  2  the 
accuracy  and  the  limitations  of  th#  model  arc  outlined;,  justifications  of 
approximation!!  ussd  in  the  model  arc  given  along  with  estimates  of  the 
resulting  errors,  . 
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6.  1  lonoapherlc-ModeTDevelopment 

For  several  year*  acientiata  have  investigated  many  different 
approach**  to  modeling  the  ionoapheric. profile  on  a  theoretical  H*ia, ,  T.he 
name*  and  type*  of  the**  .method*  are  well  iqiowa -end  will  not  be  dUcuMed  here, 
but  it  ia  obvioua  after  all  the  year*  that  a  good  theoretical  ionoapheric  profiia 
still  doe*  not  ex-. at, 

The  object  of  our  paat  inveatigationa  was  to  come  up  with  an  ionoapheric 
profile  that  could  give  much  improved  raaulta  for  refraction  correction*  in 
aatellite  communication*  to  ground  or  to  another  aatellite  than  had  bean  obtained 
with  the  Chapman  and  many  other  theoretical  profile*.  It  would  have  been 
pointleaa  for  ua  to  ait  down  and  investigate  another  theoretical  approach  when 
10  many  more  competent  acientiata  are  working  on  thia  problem.  For  thia 
reiaon  we  decided  that  in  thia  preaent  time  of  computer*, an  empirical  model 
taken  from  a  vaat  data  baae  may  provide  ua  with  the  profile  we  were  looking  for. 

It  waa  our  intention  to  acquire  ionoapheric  date  of  any  kind  that  helped 
u*  build  up  a  data  baae  covering  minimum  to  maximum  of  a  aolar  cycle  and 
providing  information  up  to  1000km,  The  lower  layer*  of  the  ionoaphere  were 
neglected  in  termi  of  their  ir regularise*  although  their  electron  content  waa 
added  into  the  larger  F  layer;  thia  waa  done  to  aimplify  the  approach  and  a* 
the  prime  objective  waa  to  obtain  refraction  correction*  through  the  ionoaphere, 
or  at  least  to  a  point  above  150  km,  euch  an  elimination  would  not  be  very 
detrimental, 

Data  from  bottomaide  ionospheric  aoundera  was  obtained  over  the 
year  1962  through  1  969  covering  14  station*  approximately  along  the  American 
longitude*  having  geographic  latitude*  76  degrees  to  -12  degree!  or  magnetic 
latitude*  85  degrees  to  0  degrees.  Thia  data  waa  in  the  form  of  hourly  profile* 
of  the  ionoaphere  up  to  the  f0F2  peak.  Topalde  aounding*  were  acquired  for 
the  year*  1962  to  1966  covering  the  magnetic  latitude  range  86  degrees  to 
-75  degree*  and  providing  electron  density  profile*  from  about  1,  000  km  down 
to  a  height  just  above  maximum  electron  density,  A»  the  topeide  data  waa 
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not  available  near  the  solar  maximum, electros  density  p»cb#  data  wet  _ 
obtained  from  the  Ariel  3  satSlilt#  over  the. period .May  1.96 7  to  April  1.96 8 
from  70  tie'greet  north  to  70  d*gr»ts  eouth  geographic  latitude  and  linked 
in  real  time  to  foF2  value*  obtained  from  13  itatiom  on  the  ground, 

6,  1,  1  Ionoepherlc  Profile 

In  order  to  analyte  the  vett  amount  of  data  that  wat  obtained  a  number 
of  aieumptioni  had  to  be  made,  In  the  firtt  ■  the  topatde  aounding  date 
did  not  geograpnically  cover  the  entire  gl./uU  and  the  bottoms  Ida  date  wai 
only  eveilebla  for  land  maaset  and  not  over  the  ocsansi  however,  a*  a  local 
time  effect  ii  far  more  significant  than  a  longitude  effect,  the  data  waa 
analysed  a*  a  function  of  latitude  and  local  time,  Geographic  longitude  was, 
however,  taken  into  account  for  the  determination  of  maximum  electron  density 
by  uiing  the  ITls  coefficients  for  f0F2  which  are  a  function  of  latitude,  longitude, 
time  and  lolar  activity,  Secondly  a  theoretical  profile  wai  determined  to  which 
the  data  would  'it,  This  profile  which  is  u»ed  in  the  evaluation  diicussed  later, 
ia  *hown  in  Figure  2  and  is  the  reeult  of  earlier  work  by  Kasantsev  (Reference  7 
and  unpublished  work  of  Bent  (1967)  while  at  the  Radio  and  Space  Re »e arch 
Station  in  Englr.nd  and  requires  the  knowledge  of  the  parameters  k1  ,  k,  ,  k3  , 
y  \  >  y  i  >  foF2,  and  h,  ,  The  equation  of  the  upper  topeide  ie  exponential,  namely 

N  =  Nce**'  , 

the  lower  ionoiphereia  a  bi -parabola, 

N  ■ 

and  the  top  and  bottomside  are  fit  together  with  a  parabola, 

.  N.  (  1  .  i  )  , 

y» 
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where, 

N  it  the  electron  density 

Nt  is  the  maximum  value  of  electron  density 

N’c  is  the  maximum  electron  denaity  (or  each  exponential 

Layer  ■.  •  "-■■■'  . ... .  _ .  • 

a  and  b  are  vertical  distance* 

y  is  the  half  thickne  ss ,  of  the  lower  Layer 

yt  is  the  half  thickneas  of  the  upper  parabolic  layer 

k  is  the  decay  constant  for  an  exponential  profile. 

The  upper  parabola  extends  from  the  height  of  the  maximum  electron 
density  up  to  the  point  where  th*  slope  of  the  parabola  matches  the  slope  of 
the  exponential  layer.  The  data  investigated  included  over  9>Q,  000  topside 
soundings,  6,000  satellite  electron  density  and  related  f0F2  measurements, 
and  over  400,  000  bottomside  soundings. 

6,  1,  2  Topside  Ionosphere 

The  initial  approach  was  to  take  the  topside  soundings  and  break  them 
down  into  tones  5  degrees  of  latitude  by  40  minutes  of  local  time  eliminating 
data  in  the  same  zor-s  that  have  similar  times  and  profiles,  and  therefor* 
are  duplicated,  1  h i *>  resulted  in  over  1,  200  different  areas  in  the  northern 
and  southern  hemisphere  with  a  reasonably  constant  density  of  data  in  each 
area.  By  these  means  it  was  possible  to  investigate  the  decay  constant  k 
in  the  exponential  topside  profile  as  a  function  of  local  time,  latitude,  solar 
flux,  sunspot  number  and  season,  Ons  of  the  major  concerns  was  whether 
the  decay  constant  k  would  bs  uniform  for  each  sounding  over  the  rang* 

L,  000  km  to  thy  minimum  height,  and  investigations  showed  that  such  an 
exponential  profile  does  not  exist,  The  layer  was,  thereto. -a,  divided  into 
three  equal  height  sections  from  1,  000  km  to  the  minimum  recorded  height 
and  the  exponent  k  computed  for  the  center  point  in  each  section,  Figure  2 
shows  such  a  division  where  the  values  under  investigation  are  the  decay  con¬ 
stants  k,  ,  kg,  kg,  In  most  cases  the  topside  soundings  do  not  rsach  the  height 


of  maximum  alactron  deasity  and  therefor*  the  gradient  at, this  lower  point  was 
mathematically  equated  to  the  point  where  the  gradient  o£th*  'noae'  parabola 
wad  the  ume,  Exta&iiva  analysis  of  the  acquired  data  ^owid  theae  gradient* 
to  be  similar,  on  »veJ?ag|c  %t  a  hsighf  yt  /4  above  the  maximum  electron  denaity. 
i\ t  thia  point  the  value  of  fkF2,  which  define*  the  low.est  jjgint  of  the  topiide 
founding,  if  0,  03  f0F2,  (Np  in  Figure  2  i*  the  equivalent  electron  denaity  to 
the  frequency  fkF2), 

For  an  initial  tast  the  decay  coaatante  k  for  each  of  fha  three  layer*,  upper, 
middle,  and  lower  topside  were  plotted  aa  a  function  of  magnetic  latitude  and 
fkF2.  Values  from  th*  northern  and  southern  hemi»ph||*s  were  treated  indepen¬ 
dently  at  first,  but  the  analyaia  showed  that  thara  was  excellent  correlation 
between  the  two,  Figure  3  shows  th*  relationship  between  th*  three  decay 
constant*  k  and  magnetic  latitude  for  all  local  times,  sol^r  activity,  and  season. 
The  equatorial  inomaly  and  a  40  degree  trough  show  in  the  lower  topside  layer, 
The  65  degree  trough  is  not  as  evident  as  it  is  whsn  the  same  analysis  is  done 
for  various  lcciil  times  which  suggests  th*  physical  varlancaa  of  these  anomalias 
should  be  investigated  in  mars  detail. 

It  was  foend  that  correlation*  in  k  for  specific  fkF2  did  not  bear  any 
further  local  time  correlation,  hut  bore  a  significant  variation  with  solar 
activity  and  magnetic  latitude.  However,  the  correlation  with  solar  flux  wsi 
considerably  bitter  tfean  that  with  sunspot  number,  even  allowing  for  the  delay 
in  the  affect  re  tching  the  ionosphere,  so  all  further  correlations  were  with 
the  Ottowa  10,  "  cm  solar  flux,  All  thas#  correlations  wirs  then  plotted  m 
graphical  form  to  enable  final  interpolation, 

Ujafortum-tgiy  th_  Aluuett*  data  did  not  cover  the  period  at  the  peak  of 
the  solar  cycle,  but  the  Director  of  the  U.  K,  Radio  U  Space  Research  Station 
made  available  electron  density  data  from  the  Ariel  3  satellite  to  cover  this 
period.  The  data  had  already  been  reduced  thoroughly  and  the  satellite  electron 
denaity  at  abottc  5  3  0  km  wa*  provided  with  th*  sub- satellite  f0FZ  value  obtained 
from  13  station  around  the  world,  If  th*  satellite  was  not  directly  over  an 
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ionoeonde  at  the  time  of  observation,  the  f0F2  values  from  two  or  three 
transmitter*  in  the  general  area  had  bean  interpolated  in  time  and  position  to 
giv#  the  sub. satellite  value.  These  interpolations  had  been  carriod  out  taking 
care  to  modify  the  valuta  for  uneven  ionospheric  gradients.  Data  that  wn  in 
doubt  was  eliminated.  While  these  values  did  not  giv*  the  three  exponential 
decay  conetnnts  at  each  point,  it  was  found  that  for  similar  conditions  of  solar 
flux  and  position,  ths  Ariel  3  data  fit  very  closely  to  the  profiles  deduced 
from  Alov  atte  1.  The  profile  equations  developed  for  the  lower  solar  activity 
period  related  to  the  topside  soundere  could,  therefore,  be  extended  to  the 
larger  solar  flux  values  and  still  be  in  good  agreement  with  the  Ariel  3  data, 
Typical  results  from  this  analysis  are  shown  in  the  graphs  of  Figure  4,  The 
original  data  UKck  were  less  regular,  and  since  the  variations  ware  mainly 
caused  by  the  relatively  lew  data  density  in  each  group  after  division  of  the 
large  data  base,  the  data  was  smoothed  by  the  fitting  of  straight  linas.  In 
order  to  interpret  iheie  graphs  and  obtain  a  profile,  we  need  the  value  of  f0F2, 
and  the  magnotlc  latitude  position,  These  values  will  indicate  which  graph 
relates  the  10,  7  cm  flux  to  the  decay  constants  k  for  the  upper,  middle,  and 
iower  portions  of  the  topside  ionosphere,  f  igure  4,  therefore,  shows  the  basis 
of  obtaining  the  3  independent  slopes  of  the  topBide  ionosphere  as  a  function  of 
f0F2,  latitude,  and  soler  flux, 

A  furthsr  correlation  to  investigate  the  seasonal  effects  on  k  was  carried 
out  with  some  15,  COO  totally  different  Alouette  soundings  and  fluctuations  in  the 
k  values  of  t  15%  were  noted  from  the  average  spring  nnd  autumn  values,  The 
seasonal  variation  is  monitored  by  observing  the  chang »  in  the  daily  maximum 
solar  zonith  angle  from  the  equinoctial  mid-day  valua,  Figure  5  shows  the 
seasonal  fluctuation  in  k  for  each  of  the  three  layers  in  ti  e  topside  profile, 
There  is  coniideiable  evidence  that  this  seasonal  relationship  has  an  added 
local  time  factor  and  this  point  will  ehortly  be  under  investigation, 

Examination  of  the  upper  part  of  the 'nose'  of  the  N-h  profile  is  difficult 
because  topside  sounding  information  rarely  gives  any  values  in  this  region, 
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Evidence  from  many  leading  ecientieta  aleo  imphiee  tha.t  topaide  profilet 
have  about  a  +4%  error  in  the  effective  distance  from  the  eounding  eatallite 

c 

indicating  the  obtained  topeide  profilet  are  too  low  near  &e  peak.  Thit 
evidence  it  bated  on  comparitont  with  two-frequency  data,  backicatter 
reiulti,  Faraday  rotation  and  overlap  teat#,  ate,  Preliminary  reaulti  in 
thit  empirical  model  thowed  that  a  parabola  in  thie  regioa  gave  the  better 
compariton  with  integrated  total  electron  content  when  com  pared  with  two- 
frequency  and  Faraday  rotation  data,  A  aimple  parabola  having  a  half  thick- 
naei  y,  wa»  fitted  between  the  bi-parabola  and  the  exponential  layer.  Upon 
initial  teat  y,  wee  »et  equal  to  the  half  thickneie  of  the  bi-parabola  y,  for  foF2 
valuee  below  10  5  MHi,  and  y,  increaeee  with  f0F2  valuea  rieing  above 
10,  5  MHn,  Furher  inveetigatione  of  thie  problem  are  planned  in  future  work. 

The  final  itep  in  predicting  the  ehape  of  the  ionoephere  1*  arranging  for 
the  gradient  in  the  upper  parabolic  layer  to  be  the  eame  ae  the  gradient  in 
the  lo  we  at  part  if  the  topeide  exponential  layer,  Thie  ie  the  case  at  a  diit&nee 
d  ■  1/k  i  (i-f-y?  k;:)l  -  1]  above  the  height  of  the  maximum  electron  deniity, 

6.  1,  3  Bottorniic.e  Ionosphere 

Modeling  'he  bottomeide  ionoepheric  profile  wn  a  eomewhet  eaeier 
t&ek  became  foi  each  profile  the  value  uf  f0F2  wai  known  and  the  electron 
deniity  veriui  height  profile  from  hgl  e  to  h|M  wai  alio  known,  Once  more  the 
geographic  effect  of  longitude  wai  eliminated  and  replaced  with  the  more  eimple 
local  time  correlation.  From  Figure  2  we  lee  that  the  equation  of  the  lower 
layer  ii  a  parabola  iquared  or  a  bi-parabola,  Thii  woe  found  in  general  to 
fit  the  real  prof  le  eomewhat  better  than  a  eimple  parabola,  The  unknown  in 
thii  equation  n  the  half  thickneei  of  the  layer  yg  and  in  the  reduction  of  the 
data  the  y,  value  wai  treated  in  a  limilar  way  to  a  topude  k  value. 

The  irreg uiaritiei  in  the  ionoionde  data  due  to  the  lower  iayeri  of  the 
lonoiphere  wen  imoothed  out  becauee  the  prime  objective  of  the  work  wae  to 
umplify  the  mo  lei,  but  keep  the  total  contmt  as  accurate  ai  poaalbla.  The 


sounding  data  was  therefore  integrated  up  to  the  peak  electron  density  (N,)  and 
forced  to  fit  the  bi-parabolic  equation  along  with  the  value  of  Nt  obtained  from 
the  sounding.  In  each  instance  the  value  of  y,  was  computed  ready  for  further 
correlation, 

A  number  of  real  profiles  from  various  stations  at  different  local  times 
ware  compared  with  the  computed  profile  and  excellent  agreement  found. 

A  further  12,  000  soundings  from  p  11  14  stations  were  analyzed  and  the  computed 
value  of  y,  compared  to  the  actual  measured  value.  These  results  are  shown 
in  Figure  6  along  with  the  RMS  errors.  The  two  tests  indicate  that  the  bi¬ 
parabolic  profile  is,  on  average,  ir.  close  agreement  to  the  real  profile, 
Investigations,  similar  to  those  carried  out  for  the  topside  decay  constants, 
correlated  y,  with  solar  flux  foF2,  local  time  and  season,  Surprisingly, 
no  direct  correlation  was  found  between  y,  and  solar  flux,  but  a  definite 
correlation  existed  in  local  time  and  also  in  the  solar  zenith  angle  at  local  noon 
which  represents  the  season, 

Figure  7  indicates  how  y,  can  be  determined  from  local  time  and  f0F2, 
and  Figure  8  shows  the  seasonal  update  as  a  function  of  local  time  for  the 
sunrise,  sunset,  night  and  daytime  period,  In  the  cases  where  f0F2  was 
larger  than  10  MHz  the  local  time  curve  fluctuated  very  little  from  the  10  MHz 
curve.  All  of  the  curves  displayed  have  not  been  hand  smoothed;  uue  to  the 
large  data  base  tha  average  of  all  values  taken  every  hour  fit  precisely  on 
the  lines  shown, 

The  remaining  unknowns  which  ars  needsd  to  computs  the  profile  am 
f-)F2  and  the  height  of  that  value;  by  far  the  most  important  of  these  being 
the  value  of  f0F2, 
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6,  1,4  Predicting  fnF2 

Severe  horizontal  gradient*  in  fQF2  exist  within  the  ionosphere  as  can 
be  seen  by  examining  Figure  9.  In  fact  even  if  the  value  cA  fciF2  is  known 
directly  above  a  station,  it  can  change  considerab.y  cvcrdfoe  whole  visible 
ionoaphere  from  that  site,  Figure  9  is  a  predicted  status  of  f0F2  over  the 
world  at  6.  0  am  luring  August  1968  and  two  types  of  sever-e  gradients  are 
immediately  noticeable,  one  due  to  sunrise  causes  rapid  changes  in  f0F2  in  an 
east  to  west  direction  and  the  other  situated  around  the  equatorial  anomaly 
occurs  primarily  during  the  afternoon  and  early  evening  causes  severe 
gradients  in  the  north  to  south  direction,  Two  hypothetical  stations,  A  and  B, 
are  marked  on  F.gure  9  along  with  the  ionosphere  'visible'  from  those  sites, 

In  case  A  the  value  of  f0F2  change*  from  11.  5  MHi  directly  overhead  to  5  MHz 
on  the  southern  horizon,  This  change  must  be  squared  when  converting  to 
electron  content  nonce  a  difference  of  a  factor  of  ever  5  in  the  vertical  content 
arise*  before  co  meeting  for  elevation  angle  effects,  Similar  gradients  exist 
over  half  the  earth's  surface  at  some  time  of  the  day  and  it  is  therefore 
imperative  to  me  del  these  gradient*  in  any  ionospheric  model, 

For  many  Tear*  NO  A  A  (formerly  CRPL  and  ITSA)  have  been  engaged 
in  the  development  of  numerical  methods  and  computer  programs  for  mapping 
and  predicting  characteristics  of  the  ionosphere  used  in  telecommunications, 

The  most  advanced  method  for  producing  an  f0F2  model  undoubtedly  comes 
from  their  work,  Jones,  Graham  U  Leftin  (Reference  5  )  duscribe  their 

techniques  on  hov  a  monthly  median  of  the  F2  layer  critical  frequency  (fcF2) 
was  developed  f:om  an  extremely  large  worldwide  data  base,  In  fact  the  gradient 
map  shown  in  Figure  9  is  a  result  of  this  work.  We  have  already  shown  that  it 
is  important  to  include  the  horizontal  gradients  of  {q^  2  in  any  analysis  and  the 
work  by  Jones  ol  al  is  undoubtedly  the  only  satisfactory  approach  to  th.s  problem, 

The  docurr  ent  by  Jones  at  al  describing  this  work  include*  a  Fortran 
program  wnich,  with  monthly  coefficient*  obtainable  from  NOAA,  enables  the 
monthly  median  value  of  f0F2  to  be  computed  above  any  point  in  the  world  at 
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any  time.  This  program  w?a  primarily  written  to  accept  monthly  coefficient, 
using  an  average  aunapot  number,  but  more  recent  work  by  Jonea  k  Obitts 
(Reference  6  ;  has  deacribed  a  more  generalized  sat  of  coefficient! which 
provide!  annua'  con  inuity  and  uaea  more  extenaive  analyaia,  .Tha.se  _g$ntr %}ired 
coefficient!  can  be  obtained  from  the  lonoapheric  Prediction  Service!,  NOAA, 
Boulder,  for  a  sunspot  numbe;  or  a  aolar  flux  approach,  The  value  of  a  monthly 
median  fQF2  can  be  computed  on  a  worldwide  baaia  centralized  around  the  apecif 
day  in  question  rather  than  the  15th  of  the  month;  it  can  alao  be  baaed  on  a 
12 -month  running  average  of  solar  flux  or  aunapot  number.  Private  communi¬ 
cation  with  Mri,  Leftin  at  NOAA  indicate!  that  the  aolar  flux  approach  is  likely 
to  provide  more  accurate  ^aluea  cf  f0F2  than  the  uae  of  the  aunapot  number, 

For  the  ionoapheric  profile  under  diccuaaion,  it  waa  decided  to  uae  the 
generalized  f 0F 2  coefficient!  from  NOAA  incorporating  aolar  flux  thereby 
eliminating  any  need  to  purchaae  monthly  data  from  them.  The  program  waa 
made  self-contained  and  enabled  a  monthly  median  fQF2  to  be  produced  above 
any  surface  position  for  any  time  of  day  or  reason  and  any  twelve  month 
running  average  uf  solar  flux. 

The  question  now  arisec  as  to  how  good  these  monthly  median  values 
are  and  how  m  :ch  error  is  introduced  by  day  tu  day  tluctuations ,  Many  daily 
soundings  were  analyzed  and  the  monthly  median  value  computed;  these  were 
compared  with  “'ho  monthly  median  predicted  valueB  and  the  actual  day  to  day 
fluctuations,  Some  typical  results  are  shown  in  Figure  10,  It  is  sear,  that  the 
monthly  median  predicted  values  are  indeed  very  close  to  the  actual  rr.tnsu’-ed 
value,  but  the  day  to  day  fluctuations  can  be  as  large  as  *75%,  A  technique 
therefore  had  to  be  derived  to  bring  the  computed  monthly  median  value  closer 
to  the  actual  value. 

It  would  be  pointless  to  U8e  the  daily  value  of  solar  flux  in  the  generalized 
coefficient  set  which  had  been  built  up  using  a  twelve  month  running  average, 
but  it  was  thought  possible  that  thure  may  be  a  relation  between  the  difference 
in  f  F2  from  monthly  median  to  daiiy  value  and  the  difference  in  the  12-month 
running  average  of  solar  flux  to  the  daily  value, 
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Approximately  6,  000  real  values  of  f0F2  from  13  stations  widely  spread 
in  latitude,  longi^u’e,  and  solar  cycla.  war*  compared,  with  the  predicted  value* 
using  the  NOAA  solar  flux  method,  A  very  surprising  r*»ult  emerged  and  can 
ba  explained  by  referring  to  Figure  11.  Eliminating  the  data  from  stations 
close  to  the  magnetic  polea  which  did  not  quite  follow  the  trend  of  the  other 
stations  a  comparison  between  the  difference  in  daily  end  12-month  Tux  value 
and  the  percentage  difference  of  computed  and  measured  f0F2  shoved  all  stations 
having  a  very  s  nilar  bias,  Figure  11  shows  this  comparison  where  the  stations 
having  similar  latitude  were  averaged  quoting  their  rm»*n  magnetic  latitude,  The 
fact  that  the  mas  did  not  pass  through  the  zero  points  in  the  graph  undoubtedly 
indicates  an  srrouetmt  bi’9  in  the  NOAA  predictions,  but  results  help  one  to 
update  subat;  ntially  the  monthly  median  fQF2  value  on  a  daily  basis,  Further 
comparisons  were  carried  out  with  two  years  of  hourly  f0F2  value3  obtained  near 
solar  maximum  from  H  wail  end  the  results  fit  perfectly  m  the  latitude  position 
expected  in  I  igure  11,  By  these  means  it  is  possible  to  come  somewhat  nearer 
the  actual  daily  value  of  f0F2,  Further  accuracy  can  be  derived  by  update 
from  stations  within  the  general  area  if  this  is  available  and  the  investigation 
oi  this  appro.. ch  will  now  be  explained, 

In  order  to  investigate  the  size  of  an  area  from  which  ionospheric  values 
'veil’d  show  s  milar  cjeviations  from  normal,  many  comparisons  of  throe  or 
more  stations  were  investigated  for  random  dates,  It  is  well  known  that 
magnetic  disturbances  can  effect  the  ionosphere  above  one  station  in  one 
direction  and  a  nearby  station  in  an  opposite  direction,  For  this  reason 
investigation!  of  disturbances  were  not  carried  out  near  to  the  magnetic  poles, 
Over  10U  proipa  of  stations  frerr  various  continents  and  having  similar 
longitudes  wera  compared  in  similar  ways.  Figure  12  is  a  typical  result  of 
such  a  test  aid  shows  f0F2  disturbances  being  recorded  simultaneously  at 
s.tes  1  ,  COOkn l  apart,  The  perce nt age  error  in  the  predicted  f&F2  value  when 
compared  to  he  real  value  was  noted  to  be  similar  in  90%  of  the  cases  where 
stations  were  within  2,  000km  of  one  another  in  a  longitudinal  direction  od 
inve Jtigationi  over  the  quiet1  North  American  continent  shew  improvement 
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in  9  out  of  iO  casta  when  foF2  was  updated  with  information  from  across  the 
continent;  or  3,  000  to  4,  000km.  However,  in  general,  the  update  procedure 
is  restricted  to  information  from  within  2,  000km  of  the  evaluating  station. 

6,  I,  3  -  Predicting  the  Height  cf  tha  Maximum  Layer 

In  order  to  predict  the  real  height  of  f0F2  the  wi(JOOO)F2  predictions  from 
NOAA  were  used,  To  explain  the  terminology; 

M(3000)F2  ■  M  FACTOR  ■>  MUF(3000)F2  /fQF2, 

where  MUF(3Q00)F2  it  the  maximum  usable  frequency  to  propagate  by 
reflection  from  the  F2  layer  a  distance  of  3,  000km.  The  M(300Q)F2  predictions 
can  be  calculated  on  a  monthly  basis  from  a  generalised  aet  issued  by  NOAA 
and  provide  the  monthly  median  value  as  a  function  of  sunspot  number. 

Knowledge  of  this  factor  along  with  the  f0F2  value  enables  the  height  of 
the  lnyar  to  be  calculated  using  the  equations  of  Appleton  U  Beynon  (Reference  1 
If  M  is  the  M(3000)F2  factor  and  one  assumes  that  y,  divided  by  the  height  of 
the  Dattom  edge  of  the  lower  layer  is  greater  than  0,  4,  then  it  is  possible  to 
derive  the  following  polynomial, 

h,  =  1  346,  92  -  526,  40M  +  59,  825MS  , 


whero  h,  is  the  required  height. 
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Fig.  3  The  mean  fluctuation  of  the  decay  conitxnt  k  with  magnetic 
latitude  for  the  upper  (U),  middle  (M)  and  lower  (L) 
portion*  of  the  top*ide  lono«phera. 
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Fla,  9  The  predicted  global  itatm  of  a  monthly  median  fxF2  at  6.0  a.m.  UT 

Auguit  1968  ihowlng  area*  of  vlilblllty  for  two  hypothetical  ground  itaflorn. 
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Fig.  <  I  An  error  in  the  NOAA  fQF2  predictions  as  a  function  of  magnetic 
latitude  and  daily  solar  flux  minus  the  12  month  running  avevqge. 


u .  L  Model  Accuracy  and  T  mutations 


As  a  means  of  testing  the  accuracy  of  the  model,  an  intense  comparison 
with  Faraday  rotation  data  has  been  performed  as  well  as  tests  with  two 
fiequency  data,  actual  ionospheric  profiles,  and  use  in  orbit  determination 
programs . 

Remarkable  improvements  have  been  noticed  in  precise  orbit  determination 
rystema  and  the  model  has  reduced  the  numbe-  of  iterations  needed  for  the 
program  to  converge  as  well  as  the  size  of  the  residuals  by  upto  a  factor  of  four, 
Excellent  results  have  beer,  noted  with  orbit  programs  using  elevation  angle, 
range,  and  range  rate  systems. 

The  most  extensive  tests  were  carried  out  by  comparing  Faraday  rotation 
data  for  seven  stations  from  Hawaii  to  Puerto  Rico  to  Alaska  looking  at  thr. 

ATS1,  ATS3,  and  SYNCOM3  satellites,  In  all,  ever  100  station  months  of 
continuous  data  were  used  during  the  years  1  96  5  and  1  967-  1  969  with  data  taken 
ev«ry  hour.  The  integrated  model  data  was  compared  with  these  actual  results; 
vour.ie  situations  were  also  investigated.  The  results  are  shown  in  Figure  13 
where  the  percentage  of  the  ionosphere  removed  with  the  model  is  shown. 

In  general,  between  7  5  and  90*0  of  the  ionospheric  effects  are  removed  and 
the  je  circumstances  are  for  solar  maximum  conditions, 

o  ■  0 .  1  Basic  Misconceptions  in  Ionospheric  Modeling 

During  the  course  of  developing  this  ionospheric  model  thorough  investi¬ 
gations  were  carried  out  on  a  number  of  other  ionospheric  models  as  a  means 
to  nndi  ig  their  basic  inaccuracies,  The  limitations  and  inaccuracies  were  then 
considered  ■  the  tinal  development  of  the  Bent  Ionospheric  Model, 

Among  the  oasic  simplifications  in  the  models  leading  to  inaccuracies, 
v-'c--  formulae  related  to  a  flat  earth  and  ionosphere  as  well  as  little  consideration 
for  the  height  of  the  ionobphere.  Each  of  these  approaches  causes  f0F2  to  be 
evaluates  at  an  incorrect  position,  consequently  produces  an  error  in  f0F2 
which  propagates  into  electron  content  and  the  refraction  corrections,  and  ,n 
addition  large  errors  in  elevation  angle  correction  can  result  from  the  in¬ 
correct  geometric  conditions. 
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The  bud  effect  of  a  flat  ionosphere  on  low  elevation  angle  satellite*  is 
obvious,  and  serious  problems  also  exist  for  satellites  at  large  distances, 
Elevation  angle  correction*  cannot  be  obtained  for  satellites  at  infinity,  and 
errors  of  a  factor  of  2  in  elevation  angle  still  occur  with  satellites  at  5000  km 
altitude,  The  height  hMF2  quite  commonly  change*  by  over  175  km  during  the 
course  of  a  day  at  low  latitudes,  Ignoring  the  importance  of  the  h,F2  computa¬ 
tion  can  give  rise  to  an  error  of  a  factor  of  2  in  elevation  angle  correction, 
and  at  low  elevations  also  to  &  difference  of  3  degrees  in  earth  central  angle 
between  the  observer  ^nd  the  ionospheric  point,  which  in  turn  can  produce  a 
change  in  (f0Ir2  f  of  20%, 

6,  2,  2  Errors  in  Range  Rate  Computations 

A  problem  can  occur  in  computing  range  rate  correction*  through  the 
ionosphere  to  a  sate.lite,  Many  Doppler  satellite  tracking  syjtems  integrate 
cycle  counts  over  a  few  seconds  of  time,  The  ionospheric  correction*  for  such 
a  technique  are  beat  obtained  by  range  differencing  the  ionospheric  correction* 
anc  dividing  by  the  integration  period;  hence  time,  elevation  and  azimuth 
changes  are  incorporated,  Atypical  ionospheric  range  rate  correction  can 
be  significantly  changed  by  the  sixth  digit  in  the  ionospheric  range  correction; 
precautions  have  therefore  to  be  taken  to  ensure  that  no  irregularities  occur 
m  computing  the  two  adjacent  range  corrections,  Furthermore,  the  ionospheric 
height  at  the  ray  intersection  point  must  be  computed  to  1km  convergence  in 
order  to  obtain  a  precise  ionospheric  latitude  ana  longitude  for  f0F2 
computations.  An  error  of  over  1km  in  h,  F2  will  cause  the  fQF2  value  to  be 
very  slightly  different  and  from  this  a  change  in  the  5th  or  6th  digit  .n  range 
can  easily  arise  leading  to  very  large  errors  in  range  rate.  It  is  not  claimed 
that  h,F2  hr  a  to  be  accurate  to  1km  as  true  is  an  impo s  sible  predicti.cn,  but 
the  values  of  h,  F2  should  be  cons : stent  in  their  calculation  to  1  km  convergence, 

The  theoretical  approach  to  range  rate  correction  either  by  differentiating 
range  or  using  the  deviation  angle  of  arrival  at  the  satellite  is  in  no  way  accurate. 
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The  differentiating  technique  yields  a  correction  to  an  instantaneous  measurement 
which  can  vary  greatly  from  the  correction  to  Doppler  range  rate  measured  over 

a  finite  time  interval,  from  a  fraction  of  a  second  up  to  over  a  minute's  time. 

In  addition,  the  range  rate  correction  i«  not  only  influenced  by  the  change  in 
the  aatellite  position,  but  also  by  the  changing  ionoephere  below  the  moving 
satellite,  which  has  mostly  been  neglected  in  either  approach.  To  explain 
this  fact,  consider  the  range  correction  ^R  as  given  by 

ad  knt  .  „  40.  3 

sinE  f“ 


Nt  is  the  integrated  vertical  content  and  E  is  the  local  elevation  angle  in  the 
ionosphere,  D.  ffe rentiating  AR  while  considering  the  case  where  the  satellite 
passes  directly  overhead  where  no  azimuth  change  is  observed; 
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In  this  equation  the  first  term  is  in  many  cases  the  only  one  ueod,  but  it 
applies  only  to  the  instantaneous  change  in  the  satellite  posit'on,  The  other 
two  terms  are,  however,  often  dominant.  The  second  term  is  due  to  the 
positional  change  in  the  ionosphere  and  the  last  term  represent"  the  time 
variation  of  the  ionos,  here,  For  instance,  with  a  high  satellite  moving 
east-west  across  the  north-south  ionospheric  gradients  a:  sunrise,  the 
time  variation  is  dominant  as  these  gradients  move  towards  the  west  with 
time,  For  a  satellite  moving  north- south  across  the  east-  west  ionospheric 
gradients  near  the  equator,  the  time  variation  in  the  ionosphere  is  very  small 
because  the  gradients  change  little  in  position  while  the  ionosphere  rotates 
with  time.  Thu  second  term  which  indicates  positional  change  in  the  ionoephere 
is  dominant  foi  lower  oatelLites  where  the  ay  path  to  the  observer  moves 
faster  through  the  ionosphere,  In  cases  where  the  satellite  does  not  pass 
overhead  the  aumuth  change  must  also  bs  considered, 


The  Bent  3c  lospheric  Model  was  developed  tor  general  use  even  at 
frequencies  close  to  critical  frequency  and  therefore  all  these  basic  mis¬ 
conceptions  were  eliminated  as  much  as  possible,  The  limitations  still 
present  in  the  system  are  now  discussed  in  more  detail. 
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6.2.3  Electron  Density  above  lQQQKxn 

The  topside  sounding  data  used  to  derive  the  data  base  for  thi»  mode) 
wai  taken  from  aatellitea  at  altitude*  of  about  1  000  km  and  analyei*  *howed 
that  the  ionosphere  above  haF2  ii  not  truly  exponential;  in  fact  at  time*,  large 
deviation*  from  a  perfect  exponential  layer  exist,  In  the  u*e  of  thi»  model 
it  ii  recommended  that  the  decay  constant  from  the  uppermost  exponential 
layer  i*  the  value  that  *hould  be  taken  for  all  analyei*  between  1000  and 
2000  km,  At  time*,  however,  thi*  value  wili  be  too  large  thereby  giving  a 
lower  electron  density  than  actually  exist*, 

Some  »cienti*t*  have  reported  that  10  to  20%  of  the  ionosphere  lie* 
above  1  000  I'm,  hut  there  i*  not  conclusive  evidence  to  support  thi*.  Further 
studies  aro  now  underway  using  satellite  topside  sounding  data  at  3  000  to 
4000  km  altitudes  and  the  model  will  oe  improved  accordingly. 

6.2.4  The  Uncertainty  in  the  Profile  ju»t-above  h,F2 

An  uncertainty  existed  in  defining  a  profile  for  the  area  ju»t  above  h,  F2, 
Topside  sounding  data  provided  a  profile  to  a  short  distance  above  h,F2 
and  bottomside  data  provided  accurate  values  to  the  height  of  h,F2.  In 
order  to  investigate  this  unknown  region  both  parabolic  and  bi- parabolic 
profiles  were  incorporated  into  that  part  of  the  model  and  extensive  testa 
carried  out  with  total  electron  content  data  provided  from  Faraday  rotation 
experiments,  The  model  was  used  to  predict  total  content  to  2000  km  where 
Faraday  rotation  probably  ceases,  The  mean  value  of  the  residuals  between 
Faraday  computed  electron  content  and  model  integrated  electron  content 
indicated  the  accuracy  of  the  profile  just  above  the  peak.  This  region  was 
found  to  have  diurnal  and  seasonal  dependency.but  these  characteristics  have 
not  yet  been  well  enough  defined  to  incorporate  into  the  model,  It  was  found, 
however,  that  a  parabolic  layer  with  half  thicknesa  a  function  of  f 0 F" 2  gave 
significantly  improved  result*,  but  further  work  will  be  needed  to  define 
thi*  region  more  accurately  a*  a  function  of  time  and  *ea*on. 
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6,  2.  5  Profile  Inaccuracle*  in  the  Lower  Layer* 


The  model  wa»  developed  primarily  for  u»e  near  to  or  above  the  height 
of  the  F2  layer  of  the  ionoaphere.  For  thi«  reaion,  it  wai  not  neceiaary  to 
model  the  E  and  FI  layer*  into  the  profile,  but  their  deneity  value*  were 
included  in  the  total  electron  content  below  h,F2,  Thi*  total  content  wa* 
then  used  in  the  deviation  of  tho  lower  layer  bi-parabola,  Care  mult  be 
taken,  therefore,  in  uiing  the  model  if  a  profile  of  theae  lower  layer* 
i*  required,  but  if  the  requirements  only  involve  total  electron  content  or 
refraction  correction*  for  value*  cloae  to  or  above  f0F2,  tho  model  is 
quite  accurate, 

6.  2,  6  Maximum  Limit  on  Solar  Flux  ir,  the  Dnrlvatlonof  the  Topside  Profile 

It  can  be  *een  from  Figure  4  that  the  top»ide  exponential  decay  constant* 
are  a  function  of  the  1C,  7  cm  *olar  flux.  The  graph*  shown  ir.  this  figure 
indicate  values  only  when  the  flux  i*  below  130,  Thi*  i*  primarily  due  to  a 
lack  of  largij  amounts  of  data  in  the  original  data  base  for  condition*  of 
higher  aolai  activity,  It  i*  not  recommended  to  extrapolate  the  exponential 
decay  con*tunt«  beyond  thi*  value  of  flux  a*  it  is  po»*ible  they  may  become 
negative  giving  an  erroneou*  increase  of  electron  den»ity  with  height.  It 
is  suggeitec  that  the  value  of  flux  be  kept  at  13C  even  when  measured  values 
are  larger, 

6.  2,  7  Limitations  in  the  Computation  of  h,  F2 

The  calculation  of  the  height  of  the  F2  layer  i*  achieved  by  knowledge 
of  M(30C0)F2  and  the  use  of  the  Appleton- 3eynon  equations  (Reference  1  ), 
Basically  a  parabolic  modei  is  fit  to  the  none  of  the  F2  layer  and  knowing 
the  half  thicknes*  y,  ,  the  lower  limit  h^  of  the  bottom  layer  and  the  value 
ofM(3QCQ)F2,  the  i implified  equation*  permit  the  calculation  oi  h,  ,  The 
equation*  oi  Appleton  and  Beynon  permit  the  con»truction  of  a  family  of  curve* 
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•  howing  the  variation  of  the  M  factor  MUF( 3 000 )F2 HqYZ  with  distance 

for  &  range  of  value*  to  the  height  of  maximum  electron  density  +y,  ) 

and  for  different  values  of  the  ratio  y, /h^  .  Such  a  family  of  curve*  i*  ihown 
in  Figure  14.  The  equation  u»ed  in  this  model  for  computing  h,  F2  (lee 
Section  6,  1,  5  )  is  derived  from  the»e  curve*  whsro,  for  a  particular  h,  ,  the 
M  factor  i*  constant  over  a  wide  range  of  y,  /\  .  Thie  condition  hold*  for 

y„  /h^  >  0.  4.  Examination  of  the  curve  *  ihown  in  Figure  14  indicate*  that 
in  general  for  accurate  value*  of  M(300C)F2,  hp  will  be  accurate  to  10  km, 
If  thi*  M  factor  i*  in  error  by  j-  5%,  w#  can  have  error*  in  h,  a*  large  a* 

*  20  km.  The*e  error*  will  increase  in  the  uncommon  *ituation*  where 
y,  i*  imalier  than  0.  4, 

6,  2.  8  Limitation*  in  the  Application  of  the  Dally  Solar  Flux  Update 

Figure  11  ahow*  the  results  of  analysing  thouiand*  of  actual  value*  of 
critical  frequency  againat  the  predicted  value »,  taking  Into  con  lidc  ration  the 
daily  and  monthly  lolar  flux,  Theae  value*  are  typical  for  the  fallowing 
thirteen  observing  atationi  from  wh’ch  the  data  wa»  reduced;  Godhavn. 
Churchill,  3oulder,  White  Sand*,  Hawaii,  El  Cerillo,  Kenora,  Paramaribo, 
Coco*  Island,  Bueno*  Aire*,  Hobart,  Port  Stanley,  and  Argentine  Island. 

The  only  *tation«  Lilted  that  are  not  on  the  North  and  South  American  chain 
are  Hawaii,  Coco*  Island,  and  Hobart,  but  the  results  from  these  stations 
closely  resemble  the  pattern  set  up  by  the  American  stations. 

In  using  this  update  procedure  outside  the  Amt  rican  chain  one  must, 
therefore,  bear  in  mind  that  the  pattern  displayed  in  Figure  1  1  is  not 
necessarily  a  worldwide  pattern,  However,  the  reiult*  from  Hawaii,  Coco* 
Island,  and  Hobart  indicate  that  thi*  update  procedure  can  be  used  elsewhere 
with  caution. 
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Figure  14,  Family  of  Curves  from  the  A ppleton - Beyr.on  Equation, 
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6,  2,  9  The  Errors  due  to  Neglecting  Angular  Refraction  in  the  Computation 
of  AR  and  Aft 

The  computation  of  total  electron  content  for  determining  the  ianoapherlc 
range  and  range  rate  correction  assumes  the  ray  paaaea  through  the  ionoaphere 
undeviated,  Thia  aaaumption  waa  made  becauae  the  majority  of  the  work  for 
which  the  model  waa  being  developed  waa  for  /HF  and  S  band  frequencies, 
Should  the  actual  path  length  of  the  ray  be  much  different  from  the  undeviated 
ray,  aa  will  be  the  caae  at  lower  frequenciea,  Maliphant  (Reference  8)  gives 
the  following  equation  for  true  path  length  d  and  apparent  path  length  d1, 

d  =  d'  *  a  R,  cob  E 

where  3  ia  the  angular  aeparation  of  the  true  ray  path  above  and  below  the 
ionoaphere,  R,  ia  the  radiua  of  the  earth  and  E  ia  the  obaerved  elevation  angle 
at  the  earth's  aurface,  Maliphant  (Reference  8)  alao  givea  a  formula  for 
computing  d  in  wavelengtha  , 

6,2,  10  Limitationa  in  the  Computation  of  Angular  Refraction 

In  the  computation  of  ionoapheric  elevation  angle  correction,  we  have 
used  the  technique  of  Maliphant  (Reference  S).  Anyone  wiahing  to  uae  this 
technique  at  frequenciea  cloae  to  critical  frequency  should  read  the  above 
reference,  in  particular  where  the  deviation  factor  ^  »ee  cj*  J  ia  larger 

than  0.  9,  s,  ia  the  angle  of  incidence  of  the  apparent  direction  of  propagation 
measured  from  the  vertical  at  the  height  of  maximum  electron  density. 

In  the  Maliphant  formula  the  exact  equation  for  ray  deviation  haa  been 
simplified  by  separating  the  functions  that  are  sensitive  to  distribution  changes, 
and  then  approximating  these  functions  for  a  typical  electron  distribution. 

The  resulting  functions  vary  by  only  small  amounta  with  changes  in  electron 
distribution  of  the  earth'a  ionoaphere  so  that  the  equation  mi"  be  used  for 
moat  of  the  values  of  the  deviation  factor,  However,  when  the  deviation 
factor  is  larger  than  0,  9,  the  deviation  angle  thus  obtained  should  be  used  with 
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caution  as  the  error  may  be  quite  large, 


Ray  trace  compariion  at  VHF  with  the  model  described  in  thi»  report 
have  »hown  po**ible  error*  in  elevation  angle  correction  of  only  a  few  percent, 
and  theae  or  cur  only  clo*e  to  the  horizon, 

6,2,  11  Additional  Limitation*  to  the  Alternate  Veriion  ot  the  Ionospheric 
Program  due  to  Interpolation  of  tho  Preproce*aed  f0F2  -  Table  a 

Section  3  2,  1,  1  CPC  No.  12  de«crlba*  how  the  table*  with  value*  of 
foF2  and  h,  art  computed  and  »tored  for  *pectfic  time*  at  one  hour  interval* 
around  aunriie  and  two  hour  interval*  otherwise,  and  for  the  location*  around 
the  station  defined  by  the  25  point  grid  pattern  ihowr,  in  Figure  1  of  that  section, 
f0F2  and  h,  for  any  *pecific  condition  are  later  extracted  from  the  table*  by 
interpolating  in  time  and  »pace,  Interpolation  over  »table  ior.oapheric  zones 
such  a*  North  America  provide*  quite  accurate  result*,  but  problem*  can 
arne  at  aunri*  s  and  at  place*  with  lower  magnetic  latitudes, 

A.  nurr.be :  of  *imulation*  were  performed  for  situation*  where  the 
ionospheric  grt  dient*  were  changing  rapidly  in  time  around  *unri»e  and  in 
?o*ition  around  the  equatorial  anomaly,  The  following  error*  were  detected 
when  comparin;  the  reiult*  from  the  time  and  »pace  interpolation  with  the 
actual  model  v« lue  a ,  In  general  h,  was  interpolated  to  only  two  percent  error 
or  butter  than  0  km,  The  interpolation  in  f0F2,  however,  provided  larger 
errors.  At  uur  rise,  the  grid  was  computed  at  or.e  hour  intervals  and  the  largest 
possible  time  interpolation  over  half  an  hour  orovided  or,  the  average  ar.  RMS 
error  in  (l'oF2J"  of  8%  with  a  maximum  excursion  to  16%  for  all  values  of  f0F2 
large  r  than  6  \  Hz,  But  for  critical  frequencies  smaller  than  &  MHz,  the 
percentage  values  can  be  quite  a  bit  larger.  Around  the  equatorial  anomaly 
where  the  ionot  phera  changes  faster  with  position  than  with  time,  the  grid  was 
computed  ever”  two  hours,  allowing  for  the  largest  time  interpolation  over 
one  hour ;  again  an  RMS  error  of  8%  ho  (f0F2  f  wa»  noted  and  the  maximum 
excurcion  was  1  9%, 
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1  0,  0  Appendix  1 


The  instruction  lilting*  in  thil  lection  ipecify  the  exact  configuration 
of  the  Bent  lonoipheric  Program  ION  and  the  alternate  veriion  TABGEN-ION1, 
The  main  program*  and  lubroutinei  are  lilted  in  order  of  their  CPC  Number*, 

Requirement!  for  veriion  IONi 

PROGRAM  ION,  CPC  No.  1 

SUBROUTINE  REFRAC,  CPC  No.  2 
SUBROUTINE  FLOTNH,  CPC  No.  3 
SUBROUTINE  PROFL1,  CPC  No,  4 
SUBROUTINE  PROFL2,  CPC  No,  5 
SUBROUTINE  BETA,  CPC  No,  o 
SUBROUTINE  SICOJT,  CPC  No,  7 
SUBROUTINE  DKSICO,  CPC  No.  8 
SUBROUTINE  MAGFIN,  CPC  No.  9 
SUBROUTINE  GK,  CPC  No.  10 

SUBROUTINE  DKCK,  CPC  No,  11 

Requirement*  for  veriion  TABGEN-IONlt 

PROGRAM  TABGEN,  CPC  No.  12  (subroutine!  of  the  above  lilt 

required  are  SIC OJT ,  DKSIC O, 
MAGFIN,  GK,  DKGK) 

PROGRAM  IONI,  CPC  No,  13 

SUBROUTINE  REFRC1,  CPC  No,  14  (lubroutinei  of  the  above  lilt 

required  are  PROFL2,  BETA) 
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*lOX<3hH‘-«,Fa,3J3G'H  <*,  CPfTlCAu  FRECuESCY  rOF  2  ■  #  F  7  .  3 ,  a  -  v^j) 
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C  Y*  •  HALF  T  H  I C  <  N  E  5  S  8F  TW£  BBTTBMslDE  eiPARABS,.A  M) 

C  V T  ■  HALF  T*!C<N£SS  BF  TW£  TSPSIDE  P ARA0BL A  ( M  ) 

C  XK  •  AKRaV  CONTAINING  DECAY  CBN'STaNTS  EBR  THE  LBWER*  MIDDLE  AND 

C  UHPeR  SECTJBN  BP  THE  TBPSIDE  EXBBN£NT I AL  LATER  (I/M) 
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Cf.  CRSiLuF  ^  1 
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R  A  T  ■  !  K  /  ( >  *  » •  *  !  )  •  *  ? 

CE^E*fll*PA^#  CF**p 
DEk ■  S  ’  J  R  T  1  ^  L s  2  ! 

T ; T  s»  x  \  -  ■  M*.;s,-«rc;P?**E 
T8T\*«^l  \/QFS 
I  P  (  JAP’  .  L  ' «  3  1  0  TQ  1  “C 

,)Wa\TiP^a’'.R\3«  *Nt'  vVLN\ 

C 

C  CC^^vTP  »(  a  V5  E  »ATE  CRPRFCTI01  DbAtE 

i  I  bPT ,  L  r  .  -OR  •  PRD'E'J.  1  )  r,U  Ttt  1“C 
DR  ATE  uCRAMtEror*pAT*5p#CF:/r'EA? 

DR  a  T  t  »  UR  a  T  r  -  f  R  a  ’  •  RS  3*  hrP  T  *  R  ov/^v 
l-C  CBsTlNLL 
RLT  JR“- 
E\n 
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PL?T  A K  '  P:  IS-  r  <=  I  " v  VfP'SJS  -^FJ'ihT  PPHEILE 

I  S  c=u  T  f  PL1-'?  IS  yH?»  ^Jv-(V'!  ]K  ^  r  T  E  s  J  X<  IS  1/5ETE0 

SUal?eLTIM  PuSTSH  <rcrg,  M^v^vT^o 
C I *ESS  1  5\  A<(3)>J(73J<h(?)i'^(?)>xS(?liHT[^)/E’)(5l 
DATA  LOiLl i IBLAS^i vAR</c<  ,1.  /1H  , 1h#/ 

Data  D3,L:C/"?»7it.i?4E,-if:j?i:4s:D??e,':2C?5E#R25E/3.  >1C. 

*:.?‘*F:C^iL:2.E1<17?t'*E3i?CEP*E.P»PP*EP/ 

*s*it  =  ({,, :  j 


CS'-Pw’l  PPpFiLE  CSSSTAsTf, 

0  •  -:Dl-  F :pt ( Cl* ! x< I  1 ) *Y- 1 *♦?  )  ! /X< (  :  ) 
*T ( 5  )  •  M'-'-v-- 

l“T  (  4  )  * 

“T  (  3  )  i 

DE  l  K  ■  (  '>  1L  1  c F  •  h t  (  3  )  )  / L  3 

iT  (  2  )  «,-T(j)OF_* 

hT( • > .hi ( ? )OElm 

ED  (  5  )  «  Ll?«E  *P"CF?**? 

ED  t  a  )  *t  D  (  T  ! 

ED(3!*E  :(-!•(  L'l-C/YTU.n 
E  D  (  2  )  ■  E  3  (  3  )  *  ExP(X<(1)*i^-,p(3)-hT(?)  j  ) 
ED(i)^3(el»  ))) 

i s : t ; a l.. : 7 e  „or=  f  ^  p  l  : > 
m  f 1 5 ■ c i:?hc 
h  t  ? :  •  l  r :  ?  b  r 
I  h  (  1  )  .  1  D  E  S 
I  h  (  ?  )  .  2  D  E  h 

3  ?  ;  3  c  1 . : , »  r 

3  5  pc  «. . :  /  ? 

”1  «■)»"!<;  -  '.‘cPF 

’  H  (  K  I  «'..(<)  .  ’r'b 

C^P.^r  ..;CTP‘-k  DE^SI-y  AT  -  E  !  3  h  ’  ^ 

L5  I C  u ■ 1 / t 

:  p  (  H  !  K  '  .  r-e  ■  ht  (  _  ;  ~,p  T.  ?c 

iC  CSV’IS-E 
Z  S  ■  C  3 
L*  *■»  *> 3 

EC  Ch.Y-f.-)"’-1!.) 

3  5  TT:  (3Li<*C/SL/fcC/7C)Jt 
3C  Z  S  *  t  (  -  X  (  3  )  •>  ; 

3  ?  ’  J'.  M  H 

<•  L  Z  S  * 

35  Tm  ko 

5C  ZS«  E X (  -Xk  i  1  )  .3”  1 
35  TP  e;j 

6  C  Z  S  •  G  1  *  i  3  ►-  /  v  T  )  •  •  ? 

35  Ta  *3 

7C  Z\.(L1"(M-i'h/vy!**?)**E 
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gC  XK ( <) »fc D ( L  ) 

9C  CSMIMjE. 

C 

C  PL?T  AM  P P  I  S' T 

XNL^O 

I  F  (  XM  1  !  ,LF.0C>  u*  T3  100 

xm  ■  l y 0 1 c  t  x s  <  1 1 ) 
ice  cevTisufc. 

db  no  » » n 7 3 

UC  J(L  >  •  IMaM 

NB  « ( XNL -U 1 C ) *3?7+0l 
ir  (N?*LT.l.’iP.‘v?-r!T»735  1*  T9  130 

J  (  '  B  )  ■  M  A  fi!  K 

12  0  rnR  1  T  2  {  6  >  2  !  n  (  n  •  Ji  x  ^  (  1  1 1  1*  (  2  )  fXK  (  2  ) 

130  CBM  I  SUF 
vs9  T  ^  F  :fc>3  ! 

1  F0P*AT  (X,l]hif!(,Hl  (  K.  -. )  ,B1Xa 57wVFRSj5  fi.tCTRBN  [U  M  I  T  Y  (f/w**3 
•  )  h  L  ;  i  - '  .3,  Lu.BSASITy  ! 

£  CJPmaT  (  IXi  1  hi  *<73A1  ,  4r*..  in’ll«*»i5X/  lAjfH  *  0  1  1  .  A  ) 

3  ABPwATl7A;?(lH^n7(lH-l,lHY#6(lH-)  1>1  Hf  ,  1  7  (  1  K-  )  i  1  H  +  / 

*bX>5"l*t3Ci  2  2X  /  •Elli2?XiB-il‘E’l?/ 

OCX,  47-M  C-.SOA.F  -  FLtrT^B'  OEMITV  (EM. *3)) 

*ET^, 
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C?Mw-TE  FP?ST  RaRT  rf  Pk^PlLPi  CRITICAL  FREQUENCY  FOP’S  AKO 

C35«FS3»,',j  I  No  HEIGHT 

SiiORliLMM  RRBFtl  (PL^T*FuSS,h.LFV.  A Za T  !  *Ea OPLUX » U# 0ma 

*  ■>  LAT1BL5\ir0f'S>HM^i.*T) 

DIMiNSlBN  <  (  10!  #'-(  I3i76) »  <N  (10Wj“!Ri‘*R!iC9T<fe)<S:T(fc)/P(3)/CBMi3) 

•  >C(3)^(7f)^OF(7fc)jG,<t4S)>0MM*R) 

31MENSI5  >  0 j ( 1  )  , CFNT !3) 

0ATA  K/U,  }a»S3»63#ft7,h9,?l,  73/75,6/*  K*,1CM/ 

0  A  T  A  KS/li7,l3>?b<37M8/B5it0^^>7  2/iNFP^\Hr/7^,(.9/ 

OAtA  Cl*:.:  0C3*  03T5/1.  ilCCO-  a3qOCCQ»  / 

0ATA  OlSvi  s V3.  1*1592653*>  ,1.C?97**426  ,ir|&69219  A-»l7Sq5R65  / 

OAT  A  RaSPi,‘1'T,C  =  ,.AT,P1.Q\/A,37i«2F  3  a»97°924(SaM  99  3684/5*078?  DA/ 

OATA  r<  1  * '  .j‘  i  r>  3/  1  l<*fe  •  92  ,l'?fr.4  ;K9»8?5  / 

OAya  Pt-,,:;- S’/ ,001  33  »  1  •  03c  »*9Bi7  #»a  / 

^  C  3  )  «Ga  r  5 

Slat*  Si',  r u * T  1 
CLAT*  CSs  RJ'i 
">EL«  S  !  S  ( £  v  ) 

CE-«  C 7 s  ,  '.'..tv  ! 

3A 2*  Sis;  k’ ) 

c a z -  cog;  k? ) 

c  t^wvT.  r  r  ^  F  0F3F.  S0f\'r  £  ,sCT!9sS  PS?  Fer?  *s0  ^3000 

T.*rF-ji  J ; 

C  A_G  3  I  O1*  O  (  (••  r°T,  c  lr  t  T  ; 

ca_l  o‘-s:  :■■•('-- p » -a ( i :i ^ r p, ob t/ > 
cap.  d< s : : a- t nw* ^ ^  1  c 

^•iTwOr,  _,T MjT.Oe  IMr  I  S' SS^-E® i C  pSIS^  PuAT,8L8S 

23  t :svf 

S  F  .  R »  c  h  .  /  !■  ♦  J  (  3  )  ) 

C  r «  s  0  "  ’ ;  i  •  S  p  *  8 c  1 

SA  «  Cb  .  *  0r  •  S-EC  * 

CA  ■  S3.  »  Cr  -  C  F  „  *  S» 

3S.AT«S^A1*CAfC.AT.SA*rA? 

Csp’»  SO?*  <  01-Ss^a  T‘S'  T  ! 

9  „  A  T  •  '''p.(ijRp/C‘PT! 

SO>ps«bA<;*SA/C'  pT 
C  0  l  6  s «  S  r.  ^  ?  ( '!  1  •  y  0  L 61  ^  *  S  G  l  f  S  ! 

9w»\iFLS'-*  A  T  A'  {  SOL  PS/:;,!.  P-  5 

Ct'^.TE  (J  •*  G  :  ■’■  I  ^  'JF,Pf‘Ok‘T  P.  S't  T  !  BSF  rv-'  PQrS  *'0  *3C0C 
pi : ) «blat 
?!  ?  !  »Bl6‘ 

call  Vi:,F  i ' 

T-PbCP''  (r).CMM?)+CI°M(3:  •CI-"  I  3  ! 

C?c ) *p < 2  ) 

C  1  3  !  » p  I  1  ) 

Cl!)*  ATaM  aTas  t  ■ct1*  I  1  )/  !'-4',;™Pl!/  S';»T  (  r^LAT  )  ) 
call  u * ( ► a  c / s  > 
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I 

t 


<K  ■  C 

oe  is  i  i.i,  m,p 

I  I  ) 

:  2-KS,  <  I  1  +  1  ) 

D9  15  j »  1 1  >  ;  ? 

<K  *  ♦  1 

15 

Ct,vp,  Tf  n  3  [,  c  0  AkC  HplRWT  SF  UAX  ,  F'LE’CTRP1'.  D  FN  S  !  TV  m1*1 
CALL  L*Cf  (‘  vp/'!M<Dw>H3oCC) 
hT  »(F,l^g4-3CC0«'H?*H3yDC*F  3000} *c IOC 0 
I  r  (  A55(E.;3i-  w«  )  ,lT»U10CC  1  G5  T6  2*. 

P( 3) ■  h- 
35  T9  e 3 

cet,l\Tfc.  err?  asc  adjust  rOF?  rpR  CAJlv  variation  *ith  fi  ^ 
2u  cr v ’ i \» • 

call  SKOnc.FFirwJ'F^CFai 

S*i  ■  5  *■  _  AT  *  qPLAT  ♦  C^L  AT  «  CPLaT  *  CR$(  ?L  *  N  - p  l  9  S  ) 

cml«  smft  ( i -sml*s*l ) 

hLat.  atm'  (SmL/Cm1-) 
uAT  1  ■  1 

L  A  ^  S  ■  1 

IF(hLAtI0EO3(_A12)  )  Gr  TO  9\ 

LATP  .  ? 

;r  (  ><l*t  k;'i(  _  a'2  } )  fi?  t*  ?\ 

L  AT  1  •  c 

IF'hLAT.-L.OKuA-R)  l  Gr  T9  pi 

LAT?  I  S 

I f < - lat.^",:'3{-ats)>  r,e  75  ?i 

L  A  T  1  •  3 

21  CS”  «  C  F ‘,  ’  !  ^  A ^  1  ) 

;f  (  ,A*1  -J.I.A-?)  30  TP  ?? 

c\T  «  ’  *  !Cf\T(l*t«!) -rr>  ■'•(l* T:  i  )  *  (og(lati  ;  -h. at i 

*  /  . r'  j  (  La ’l  )  -LG  ( i_A’«? )  ; 

2?  F  Or  2  ■  FGFT  *  (PLR*CFluX  ♦  CNT  ) 

-<FT„RS 

E\r' 
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CP'-Pj’L  SrC°S3  "ART  Of  »r?Hr  T  l.  E  {  PaRa-TT-RS  yw/yT,XK<  SATJ9  eF 
EL.C'.'nTsT  r  M  F!L  •  D  C  N.  O I  yy  X,\|T»"S  7  A\GF  RATF.  “ULTJPLIFR 
SJPR3JTI\t  -’Pgr  L,2(5L*T<t!i.ek<6S/,fIME/!0Ay<M?'j,rF_wX/F  0F2<HI4<HLAT< 

<  VN1/  yT t  X< <  7  R  H  <  X  \  T  S' M  ) 

D  I  wENa  1  yyTAi>(l?.9),,CE;p,'f‘*i3,3)iSL*'p!<*<313^w(4)<Dh!3!>X<,3) 

5jyE\SI0\  ~E3  (  3  ! 

Dlk‘FSSl9‘.  paT<(4<4<3) 

D!kTk.Sl;is  vraT(7 <6! 

DATA  S0<C.:,  J2<  :3>'5A/J5('ih<:4A>^P4,r>37/JlOO0,:>p.O5i3P1133/3P9^A32p5i 


#310R5<  l*8Sir;/  C 
■w  3  7  •  <  1 3  C  3  *  <  •  0  ^ 


<  1  •  < ? *  <3i 

<•133333  <*95 


<  4  •  <?,  /  6  *  <8*  <  2  4  i 

<  ?  •  5  >10*5  <  *933333333 


:a‘ 


'*<  ^  lC<0i6/03C,'j  13501  80/ p!H<Pr2<0E3/ 


■  3F7PA6..A2-  <■  1  308996917*  <•  n96?fc  3*  /•  1*7*532925  < 

•  279?52(>b  <.??!59*775  <2.15619449  <  3 . 1  4  1  *5926536  < 

10  7079*53  268  ,  6*?K'jl893C7F  <1*3389969375  /.7853981625  < 

.  2  A  1  7  9  9 1 A  7  5  / 

3aTa  SM<.S^2<pk,4<Hl3l2/*409*749593  <*0172142063  <.  9375  < 

1 0 1 2  0  0  C ,  / 

OA'a  Cr>  ',/iA.2r-fc<H,73r-6<  ;  ■* .  *♦  5  F  •  6  <  <c.45F«A<8ifi8E-6<  10‘3iiE'"fc<9.E*,6 
<13,9**L*'><2.3F,6<9.56t«!j<il*'7E“A<i?.A4F»6<7.62fc',6<4.67r,,6<f5»86F»6< 
5.  *  6 1 «  6  <  1  >  ‘  ;‘*fe<  3  •  1 4  r  *  6  <  .,i4P-"6<9.08fc>6/  *  8  F  •  A  /  3»  1 6  E  ••  6  <  4»4lE"6/4»54E 

■  6/  3.92  *_ -  6*  J*5wF“6<4.Chi.*6<  »  •  0  A  r  •»  6  #  •  6  5  E  ■  6  /  1  *  5  0  E  ■  6  /  [  »  4  4  “  ■  A  <  1  .  9  5  L  •  6  < 
'  5  F  •  6  /  2  •  7  *» :  •  A  <  1  «  2  2  r  •  6  <  1  .  *5  *  P  «  6  / 

3  A  a  b„3F/*7»b“  «  6 / *3 • 6F  •8/«H.E',a/“9.E*fe/"3.lfc."8/ 

■  3  *  A t  *8 < ■ 1 iE“8<"l.ar-H< 

4.5r:-a<"l.AE-6<»3.5F*tt<-?.^f-8*-3.*»E-8A».4K-A<-1.3E-fc/-1.2f’“6A 
J.t-R<l.FE-*5<.hF-8<1.4r-K.2,5L.«<i.3r->5<1.3F-8<l*7r-a<-.9E-8, 
•.fc'-a<-.7r«8<«.7r-8<:.hF-feil.:F»8<1.7F»?<i.E»8<l.F-8<.Sk'»8<i.6E«8 
/ 1 • ?t •*/ 

DAta  RA»p,/,<?/,9b<1.37,i.i4..8A/.g9/,9/*,2fc;J,a8<.975<1.05il.,i25< 

.  9  4  <  1  .  : ,  !•>  '.  5<l  .llr)<'9  =  <.P7,.9*-<1.34<.94<.o«5<.975<i.cc5<lM25< 

1  •  1  i<  '.  ' :  j  a  A  >  .  *3A5<  0  9.  1  .  1  7  5 ,  *  *3«<  *9>*<  *  985  .  • n  <  *8A<  •  9  9  5  <  •  9  2  5  <  1  .  ;  b  c  < 

•  9  7  »  >  9** .  l  >  ^  2  b  <  1  •  C  9  <  1  •  3 1 :  *5  <  •  0  7  <  ."•*''<  *895/  *83  /  .8**/ 

^Ata  /‘'TA.)/87.7<5?.C<q7.a<i:2.r<lC?'7<99,,.9b.l<9l.3<«R.3<y6*8< 

8  6  •  3  /  8  b  »  c  >  9  fc  .  R  <  9  8  .  *  <  *  j  *3 .  <9  *  109»5<  112.5/  1 1 2  •  5  <  1C7»5<  101  *2.96.2/ 

95.h<97,l'<o».3<  107*A<  1  *  7  *  7  »  1  h  T  «  1  <  1 5  0  » *•  <  1  b  3  •  3  /  l54*rJ<  1 5  0  *  0  a 

1 4  *  •  2  <  1  ?  7  >  1  »  1 1  5  .  t  <  l  -  9  .  ? ,  ^  o  fa  .  5  <  *  1  h  .  4  <  1  t  9 . 5  <  1  4  4  .  ?  <  1  6  2  •  7  <  *  7  b  .  6  < 

18  0*6/  1  1  **  •  M  <  1  b  7  •  b  <  1  J4.7.  1  lt'.J<  1  10.1.  1  13. 7<  113*3/ 120*7/ *34*9/ 

1  5  8 . 2  <  U  1 . 4  <  1  9  "  .  b<  1  f  7  .  ,  1 1?  .  b<  1  5>  3 . 3 <  1  1  9 . 4 ,  M  1  •  R <  1 1  0  •  3  <  113.5/ 

125  ►**<lj9.^<lb-.b/l  99.  b  <188* 3/ 1*3. 3/166. A<  136.9/119. 9<lll.9il08.0< 
1 1 4 . 0  *  ll8*2/l?5.6< 157. 7/211. 4/232*9<21  1*2/  188.3  /  142.5/  126.8/ 

1 1  6  •  8  /  1 1  r  .  H  /  122*7, 132. ".l‘*3.3.158«3/187.1i214M4/l  96 .8/185*5/ 

152*5  /  130*7/  120*  7<  l  1  7  ,  ^ ,  iO-  8  /  1 47.  5,  1*5 » 0/167. 8/2C0*J/195.6/ 

187.3/164.  j< 144 , 4, 1 38. 7/117. 7/1 77. b/ 

0  A  7  y*'HT/«,25<:  •  1 2  /  1  •  J**  <  OK.  *  9£ .  •  9? .  »9?< 

11. 1/  1 . 06  /  *  99  /  .*8  /  .78  /  *73/  .7<  1 . 3  < 1 • 2 1 <  1  .02/  »88< . fc 1 <  .  7 1 < .78/ 

21 . OS < 1  *  0  4/ 1  *  01 / *  98/ *98 , . 99, * , <  .95/. ^6<« 97/1*/  1.04/1*0° /l* 13/ 

3 1 . 2  *•  /  1*24/  1  *  *4  1  1  .  2  4  <  l»jj<  10  3/  1.64/ 
hi,at«  a y s  ( *  v * T  ) 

Tl " C ■  T  I  -1 1  ♦‘V6N.J12 
TL*C*  Amd;:  ,  TOC  <  p  I  2  ) 
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C 

C  COMPUTE  HALF  TKICKNE83  YM 

T12«TLflC/03C 
L  T 1  ■  T 1 2 
Ti.LTl 
LT2.LT 1+1 

I^(LT1*EQ*12)  L T 2 • L T 1 
I F ( L T 1 1 3 E  o i )  a*  T9  55 
L  T I  ■  1 2 

55  Tl.T12*-l 
IFi.F0F2.QP35 
IF2.F0F2.QO05 
IF(  IP1  •  UT • i  )  I F  1 »  1 
I F  (  !  F 1  •  QT  •  9  )  I P 1 
IF  (ira.LT.lJ  IF2«1 
irtira.3T.9J  IF2M 

y*|»  t  VMTAfi  (lTI  J  IFl  )♦  (  VMTAB(.ra<  1F1 )  •VMTA5ILT1/  1F1 )  )  *T1  )  *Q  1 000 
I F  (  IF1.EQ.IF2)  38  tjj  60 

YP2«(  V AB  (  l.T  1  ^  I F2  )  ♦  <  VMTAB  {  lT2#  IF 2 )  -YHT AB  (  LT 1  i  1 F2  )  )  *T  1  >  *Q1 000 
F  1  ■  I  F  1 

YM  YM-»(YY2«YM)#(F0F2t1-Q1  ) 

60  C6NTINUE 
C 

C  C8MPJTE  DIFFERENCE  BETwEEN  AVER.  AND  DAILY  TB'.AR  ZENITH  ANBLE  DBZA 
DAY. ( m8n« 1 ) *  30+ 1 D AY  «8C 
DSZA-SBl*  SIN(S82«5AV) 

I  F (  ABS ( 9L a  T )  .LT.981 )  39  T9  61 
I  F  (  BL  AT  i  u'r .  30  )  D9ZA«*DSZA 
39  T?  62 

61  3AN3«BLAT/S81 

CANQ«  SQR T ( 3  1  •  ABS ( SAN3*SAN3i ) 

DAN3«  ATaN( SAN3/CAS3) 

ASZ A*S9l* (CAN3*3AN3»DAN3> /PTN 
D6ZA»ASZa-  ABS i  9LAT-DSZA) 

C 

c  APPLY  SeaSBML  EFFECT  9F  DS7A  T9  HALF  THICKNESS  YM 

62  S12«QA*DSZA/D8 
I  F 1  ■  S  1 2 

S' ■  IF1 
S1.S12-S1 
IF2.3FU1 
RAT. 30 

IF(HLAT,lE.Q5)  39  "9  63 

T12.tTL9C*D?B5)/PlH 

LT  l  . T  1 2 

Ti.LTl 

LT2  «LT 1  + 1 

1F[LTZ»3T.h)  LT8«1 
irtLTl  •  LT  •  1  )  LT1M 

RATl.YRAK  IF1.LT1  )  ♦  <  Y«  AT  (  I F P,  U 1  I  . YR A T  1  IFIaLTU  >•■! 
RATf.YHAT(IFl#i.Ti)-*(YRAT(lrPiLT2)-YRAT(IFl<LTI>)*il 
RAT.RATU(RATt.RATl  )•( TlB.TJ | 

I F  f  HLAT  *  QE , QE3 { 3  )  )  09  T9  6* 
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63  Ti8«UU&C*Dl3B)/DlB0*0l 
IF(T1B*GT*Q1 )  T12«Q8*T12 
If  ;T12*LY.Q0)  Tl2»*75.H 

fUTl«YRAT< IF  1/ 5 ) * ( YRAT [ IFgj$) »YRAT ( IFt/5) ) *S  1 
RAT8*YRAT{ IF  1*  6 )  +  ( VRAT ( IFg,*) *YRATtlFi/ 6) >•$! 
R*TMbRAT1Y(RAT2.RAT1>*T18 
RAT«RATm*(<?aT»RATM)*(HlAT-DT>/510 
IF  (HL.AT  .LEiDB)  RATiRATw 

64  Y5«YK#RAT 
C 

C  CCMPJTE  K-PARAMETERg  X< 

FCF2  •  RN4  *  F[;F2 
11-2 
1 2  ■  2 

I  F  (-HL  AT  *2Eq  (  2  )  )28/3Qif9 

28  M  •  3 

I F  (  HL  AT  , lE * OEQ  (  3  )')  !2«3 
3B  TC  3c 

29  S  2«  1 

IF(HLAT,gE»DEQ( 1 ) ) 11*1 

30  J  ■  (FUP2  *  31)/Q3 
XF.QC 

IFfj.QE.l)  39  ^9  3? 

J»  1 

QB  TB  **5 

35  I  F  (  j .  L  T  »  4  )  39  *8  4 0 
J*  4 

30  TB  45 
40  Fl-J 

xr  «  (fqf2  4  Qi)/03  •  n 

45  DO  51  Mil, 3 

SLR*  (  SL9PI J4  1, J 1,M ) -SLOP! J, 1 1,M ) !#XF48L9b ( J,  1 1,M) 

CPT«(CEPT(J4l  /  I  lMi-CEP*!  J,  T  !/*■<)  >#XF*CEPT(J,  [1,M) 

I F  (  I1.E5-.I2)  35  Te  50 

DfUi(  HLAT-SE3( 1 1 ! )/(0E3( 12) -0E3(  11  )  ) 

SC.PiS.P4(  (Bl9P  (  J*i/  I  2/ M ) »SLBP ( J/  12/ *>  )  *XF4SL8p(J/  12/M)  •SLP>*DEL 
CpT.car+(  (CEPI(  j*1/  I2/mj.ccpT(  J,  12/M)  )  •  X*‘*CE°T  (  J/ I  2/ M  ) -CPT  )  «dEL 

50  XK(4)  •'  SUP  •  PLY*  ♦  C°T 

51  CBM  I  SUE 

c 

C  APPLY  SEASONAL  EFFECT  8f  OS?A  TB  DECAY  CONSTANTS  XX 

T12.T.BC/DE3(3)-Q8 
IF  (  T1  2 •  LT •  30  )  T12«T124-Q8* 

T12iT12/36*jl 

LT1.T12 

T1.LT1 

LT8-LTW1 

IF(LTJ«3T»4)  L T2*  1 

31BiQJP5-0Si*/D16 

IF  1 »S  12 

SI . JF  l 

S1*S1  2 - S 1 

IF?«  3  •  Ul 
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D8  52  M#l,3 

RAT1.RAT<<  !F1*LT1,M)*(S<ATK(  !  F2  <  LT  1  *  H  )  •*  ATK  <  J  F 1  *  LT  1  *  M  )  )  *31 
9AT2«RAT<( !F1>lT2/^ )*(RAT<( ?P8*LT2#M)#RATK(lPliLT2#M) ) *81 
RAT  •  RATI  *  !RAT2-RAT1)  •  (T12-T1) 

52  XK ( M 1 ■ X< ( M ) #RAT 

C9MPJTE  HALF  THICKNESS  8F  TBPSinE  RARAB0LA  YT 
CBNV»Q 1 

I F ( TOF 2 . lE  t J10®5>  39  Tft  71 
CBN V# CP  1333# (FCr2-C10P5)+Ql 
71  CBNTINLE 
YT«C9NV# YH 

C9HPJTE  HOST  MULTIPLIER  P8R  RAN3E  RATE  CBMPUTATJBN  RRM 
CBHPwTE  T9TA.  ELECTRON  C9NTENT  /  rLECTRBN  DENSITY  XNTNy 
XNTNY#UC 

RSJ  m  a  Q  Q 

D«- (31-SSRT(Q1#(XK(1)#YT)**2) )/XK(l) 

Hfl ) «HM#2 

irinSiLE.H!  1  )  )  39  tb  8o 
RRm.31 

DElH  ■  (—1012  -  m ( 1 ) )/33 
H r 2 )  •  h ( 1 >  *  DEl h 
Hf3)«H<2)#2ELH 

Hf k  ) «HS 
M»  3 

65  If(HS*oT.H(Y)  )  38  T 8  7 0 
1|H),H(H»1) 

M«Mr  1 

;  F  (  M  .  3 T  .  1  )  (39  TB  65 
70  DM ( M  S  #H ( M#  1  ) #H ( M  ) 

RK#Q1/XK ( H  ) 

EX'CO 

ARG«XK ( ( H ) 

I F ( ARQ  »lt ' 037 )  £A#  EXP ( -ARG  ) 

RRM.RRf-*E  X 

XNTNM»Rx#EX* ( XNTNM.RX ) 

M#**  1 

I F ( M . a T . 0 )  39  Tfl  70 
TEfc'P*389l5.YY#H.D#*'j/'!03#YT#YT) 

TE^Pl »U1 - ( D/*T ) #*2 
RBm«RRH#TEwP  1 
XNTNM.TE^P1#XNTNM+TEHP 

QB  TB  110 

80  I F ( hS  t  LE  *  !  H M - y M  )  J  03  TB  110 
DJRT#  Hm.hr 

I F ( hS • lT  .  hH)  3?  T9  30 
XNTNHi089l5*YH*DlST+Dl8T*#3/( Q3*YT#YT ) 

R F M ■  0  1  -  (  |H-.H8)/YT>##P 

39  TB  110 
90  CB^TINuE 

XNTNYiJ8Hl5#YM«Dl8T  +  r32.DIBT#*3/tO3#VH**2).0!8T«#5/(DB#YM*#9) 
RR'4#(Cl-<(>-ri«HS)<'YM)##2)##p 
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110  CBNTINUE 
RETim 


r>  r>n 
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C 

SUBROUTINE  BETA(FPAT«XNTNH,H8iMH,vM#*E#CF.#MLCV) 

BETA  COMPUTES  IONOSPHERIC.  ELEVATION  ANGLE  CORRECTION  TO  BE 
SUBTRACTED  from  measured  ELEVATION  ANQLC 
DIMENSION  XAX(B>iTAX<5) 

DATA  R#  Qo<  QI/Q5333«Q2/6371i0E3#0»  »  1  *EO* «B*33£Qi B.EO/ 

DATA  XAX/D.gO,  .ZEO,  #*EO,  .6E0,  »B1E0/ 

DATA  vAx/1 #EO# .92*EC# •88AEC<*7E0# .S53E0/ 

RE«R*R 
RS-HS4R 
C 

C  COMPUTE  SQUARED  DEVIATION  FACTOR  XCOM 

RflMiR+  HM 
SP1M*R*CF/R0M 
CP  I M »  SORT lQl*SFlM*tS) 

XCOM-FRAT/CF IM*#2 
C 

C  INTERPOLATE  TABULATED  VALUES  VAX  TO  3eT  vcom 

DO  30  laws 

I P ( XC6M#XAX ( I ) 1 20, 10/ 3D 
10  YC6M»V AX ( I ) 

38  TO  40 

20  VCOM- Y  AX  (  I  )  +  <YAXU*l  )-YAX(  I  >  )  •  I  XCOM.XAX  (  I  >  )  /  (  XAX  (1*1)  *  X  A  X  (  I  )  ) 

30  TO  40 
30  CONTINUE 
38  TO  50 
40  YC8M»(Jl/VCf>M 
C 

C  COMPUTE  DEVIATION  ANGLE  ALPHA 

R08»R8M*c5333*ym 
SF I8"R4CE/R00 
CP  1 8 ■  SQRT ( 3  1-SF I  8**2 ) 

ALPnA»FraT«yC8m#XNTNM*sfIB/ I 32*R88*CF 10**3 ) 

c 

C  COMPUTE  ELEVATION  ANGlE  CORRECTION 

C A*  COS(ALPHA) 

SA*  SIN(alPhA) 

X3»R«CE*SA/(  QWCA) 

X2*R«Sfc.*X3 

Xl*  SQRT(RS**2*R2*CE**2WX3 

CTE*(X1*CA*X21/  SQRTt  Xl«*2*v2**2.3?*Xl«X2*CA  > 

STE-  SURTI  A88<G1-CTE**2)  ] 

DELEV*  ATAN(STF/CTE) 

RETURN 

50  wRITE(6,l  ) 

1  FORMAT (  112H  **•  PAY  18  REFLECTED  AT  IONOSPHERE  OR  NEAR  REFlECTIO 

*N  CONDITION,  ELEVATION  ANGLE  CORRECTION  JB  NOT  COMPUTED  *•*} 
DlLEViQO 
RETjRn 
En^ 
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SBBR8UTINE  SICfUTlUC/tiT) 

C8MPUTE  SlMJT>*C8f ( JT ) < J» 1, . ♦ F88  ANQLE  A 

31HEN8I8S  S(1)*CU) 
cee t t > 
eu>«  simt> 

D8  10  1>?<L 

C<1  >»C(1)*C{  !-l)-SM)*S<  I-l) 

10  8?  niCU)*9(l«li'HI  tll»C!!*t) 

HfTUSN 

EKiC 
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SCJBMI'TIKiE  o^SJCO  (HX#LM#0*aitTliE/CeT'lMEj3io 
CB«PUTt  0  SUB  X/CBE^JCICNTS  W  Vf  !XW  TIME 
0!*EN$!*n  DC  U  #CBTIME(  i  8ITIHP (1  )*DK(  i.) 

UHAX*LHtg*i 
LK*l*lMAX 
DB  5  K»1,MX 
LK-LK+LMAX 
DK(X).D(LK) 

DB  5  L«l#Lw 
NK»LK+L*2 

5  .DRCXJ  liJlTl^E  ;u^^(NK  5  *C0T  I  ME  <  L  > 

RETURN 
END 
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SUF3S0T ! \t  M  AGr I N { P9S/ UNF ! 

C8VPLTE  N' A •; A  MAGNETIC  p  JfLO  COMPONENTS 


r 

c 


5 

L 

C 


DIMENSION  P{  7/7)/DP{7/7)/CF><71/A0Rm/£P<7>/Pfi8C3>  »'jNE<3>/CT(?/7>/ 

•  GC7/7)/M7/7) 

DATA  CT/2*Ci# *33333333* .£6* 66667/ .2571 4gP6/, 25356825, *?525?525, 

•  3*C  •  / *20000000,  »?2f,A71-*2»  .23609523. . 242 4?4£4/ 

•  4*C.,  .  14285714,  .19U7619,  »?1?12121/ 

•  5*Ci/*1UU1.UmI'6UUU* 

•  6*0 » /  .09090909/  .  '  .  ;  -  ' 

•  7*C»/7*0»/ 

DATA  13/  0*/  •  3Q4112/  •  0HAC35>  *.  03151  8/  -  *0M 79<*>  • C 16256/ -• 0 19523/ 

»  C  *  /  .021474/-, C51293  /  .C6213C/-.045298/-.C344C7/ -*004853/ 

»  2*0. / -*013381/ -. 024638/ - *021 795/- .Q19447/  *  003212/ 

.  3*0./-* 006496/  . CC7C08/-. 000608/  ■  .021413/ 

.  4*0./ -.002044/  ,002775/  *001051/ 

»  5*0, /  .000697/  , CCQ227/ 

.  6*0,,  .001115/ 

DATA  h/7*G  »  j 

.  0,/-» 057989/  .0335  94/  » Cl  4B7Q/ - . 011825/ • *000 706/ -• 005758/ 

•  2*0  ./  -  *001579/  -.004075/  *010006/ -i002OC0z 008735/ 

•  3*0  ./  .000210/  .0004.30/  .  004597/ -  .003406/ 

•  4,0./  .001385/  »CC24?1/*. 000118/ 

•  5*0. i -.001218/ -.0011 16/ 

•  6*0. , -.000325/ 

DATA  |>(  1/ 1  )  /  TP  (  1/  1  ) ■/  SP  f  1 )  #cn<  1)  f\  .  ,0./  0./1  ./ 

DATA  i<F /RO  / ^895/637  l2C0»t  0/0*  EC  /  1  •5A90&3S98EC/ 

P?*P5:i(2) 

Pi  ■P0:3  ( 1 ) 

I  F  (  AIj Ej  (  PI  1  •  L E  •  W890  :  Gt‘  T«  .. 

PI  *  s  !'jN(»P9P/f5i ) 

P2  «  DO 

4  CjMluOg 

AR.PF.  '(PE  +  300. E3) 

C«  S I M  n  1  ) 

S»  S5l*T(f.P(  1)-:»C) 

A  8 1’  (  1  |  »  a  P  *  A  R 

CB  ‘-'PU  r  L  SlN,Cfcc  Fflr?  ML'lTIPLP  L'^GITGOE  AN3LE 
CALL  3ICeJT( A,:a(2) /SP(21/8 ?) 

D5  5  <»2,7 

A5P (M ) ■ (  m.1  ) 

CLFAP  L'UTEP  Sj"S  A"  D  Jjt  T  Jt  Lef'P 

fW.OO 

9Ni20 

ept- 1 »;(' 

05  6  >*2,7 
F\  »  N 


C 
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C  CLEAR  INNER  SUMS  AND  SET  UP  L8BP 
S6HM90 
St)MT*8C 
8cjmP*Q0 

oe  7  m»wn 

CfMPUTE  FUNCTIONS  AND  Df*lVATlVt«  8f  MULT. A8tiLC3fNDR|  FUNCTlSN 
II  THIS  LAST  C8NTRJ8UT18N  T0  INNER  IUM 
If ( Hi NE iN )  SB  T8  6 
PtN>N)»S*P(N*IiN*i ) 

DP(N*H)«8*0P(Npl#N*l)+C#P!K-i #  M- 1 5 
08  TB  JO 

I  PlNiM)PC*PtN*liH)*CT(NjMJ*PtN»?*H>  ' 

DP  t  N# H  >*C*DP(N*1#H )-■•!»(  N.i#H>  «CT(N*M)*Of 
10  FR«M*1 

T8«3(N/H)»CPlM>-»H<NiH)«8P:M! 

C 

C  SDM  INTO  INNER  SUMS  F8R  ZjX,Y 

SWHMSUHR*PtN,M)*TS 
IOMTbSUHTmDP ! N<  M ) »T8 

7  IOHP«»UHP*FM*P  )  *8P  { H )  *H  ( N>  M )  #CP  ( H  )  ) 

C 

c  SClM  I  NTS  BuTER  SUH8  T8R  Z#X,  Y 

BVaBV+ABR(N)»FN»8UMR 
BN«BN*ABR(M*8UHT 
6  B*H I«BPH1«a8R(N)*SUMP 

C 

c  set  magnetic  field  components  ?. vertical  uPiX»nbrtm,y.east 

UNE ( 1 )  «»BV 

UNE ( 2 ) »BN 

UNE ( 3 ) •*BPU I / S 

RETURN 

END 
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C 

8B&R0IJT  1 NE  3<  ( K/CjQ ) 

CBMPUTE  CBMDINATE  FUNCmNR** (I ) # I* it • . **K  +  1 
Ctl>*K8D!F ! CD  LAT!TU0CiC(*J*C(l)«ai8atU8N#!TU0t/UATlTU0C 
a  I*  THE  ARRAY  FAR  3EB3RAPH  l.C  PUNCTIBNS 

DJmEN8I8N  < ( 1 ) iC ( 1 ) i 3( 1 ) 

DATA  fit/li  /#  N/B/ 

X*C£1) 

Y*C(2) 

ZlC(3) 

<0»K ( 1 > 

sx«  sin' no 
c 

C  SET  TERMS  DJE  T9  MAIN  LATITUDINAL  VAR  I  AT  1 8N 
3*8) *BX 
3( 1 ) »G1 
DB  10  1 ■ 2 / <  0 
10  3*1  +  1) ■SX«3(  I  ) 

KD IF •<  <  2 ) »<0 
U«1 

CX1«  COS  £  Z  ) 

cx«cxi 

TpY 

18  KC«X(J)+4 

C 

C  CBMPUTE  FIRST  2  T£RMg  8F  J.TH  PRDfR  LATITUDINAL  VARIATION 
Q( VC-2) «CX#  COS ( T) 

G»K01)*CX*  SIN(T) 

C 

c  are  bnlv  2  terms  t*  be  cbmpijte-  frr  this  brder  lbnqjtude 

IF<K0IF.EQ.2i  39  TP  28 
KN«< ( J* 1 1 
C 

C  CBMPUTE  REGAINING  terms  PP  J.TM  PRDEp  LONGITUDE 
DB  22  I«aCj<N*? 

3 ( I ) »SX*G< 1-2) 

22  3<  I  +  1)'«SX*3I  I -1 ) 

C 

C  ARE  TERMS  P9R  MAXIMUM  RRCE&  L0NQITUDF  CSMPjtED 
28  IF ( U  »EQ  »N  >  39  T9  8C 

C 

C  PREPARE  r9P  NEXT  9RDER  lSNSTTlDE  CBmRUTATI9NS 
kb  IF«K ( J*2  )  •< ( J*1 ) 

IP ( KD IF . CO . Q )  39  TP  8C 

cx«cx*cxi 

J»U+1 

FU*J 

TiPj*y 
3b  TB  18 
80  RETURN 
END 
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8BBN8UTINE  D<a<(MX,a#D<STAR#GHCaA) 

C8MPUTI  B«taA4  8UHMIN3  THE  3E93RAPHIC 

OIHSNBieN  fl(  DiDMTARtl! 

BBE3A*Q( 1 )#0<STAR{ l > 

DB  5  K«8iMx 
5  GrtEGA»BMEQA-OKSTAR(K}*Q(K ) 

RETURN 

ENT  ■  . 
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•RB3RAM  T  A83EN ( IMPUT#  OUTPUT*  T A*E5" 1  NPUT*TAPE6*BUTPjT,TAPEl*TAPE2) 
PREPROCESSOR  WNE*6TlMt»  P0F1*MM  f AtLCt  0N  TAPE  TB  BE  JSED 
with  PRO-AM  ttfcl 

ITP«UNJT  A$#4#MMCA*T  9*  TNRUT  TAP*  WJTW  IONOSPHERIC  C»£Fr  J  C ICMTS 
JTP-UNIT  A§*!SM*f*T  8T  iUTniT  TAPE  W!TN  IBS*  F0F2-HM  TABLES 
DIMENSION  JAZ(4>*IFH( 14*88 >*FLX< 31) 

DIMENSION  K ( 10 ) *KN< 10 )*C9T (4 ) *5! ? f 6 ) /P ( 3  >  *CBH( 3 ) 

**Ct3>*3(76),Dr(  765  *3WU9)*DM-<  4-*> 

0IMEN8I9N  DB< 3)  * CENTtt) 

01  MENS  1 0N  WCSEPdj  13*74)  #U(lfc*7*)#UM(S#P»)AUiH*AA») 

BATA  JAZ/J *4>8*12/a  MPaJW1i|/ 

DATA  M»NDV,M»SD#LVIW*/O>10e00#0/ 

SATA  k/1 1*35* 53* 63* 47* 69* 71*77*75* 6/*  <M10/4/ 

I A TA  KM/1, 7*  13/28*37* *|i 55, 60*65* ?8/* Nrr,NMF/?6* 49/ 

DATA  01 *310* 0100*0 130#S3TB> 341 #3*8/  1*  *10.  *100*  *130* 

*  * 3O0C00 •  * • 1  i * S  / 

SATA  DR*PI2*D7,DHRi*D*R2/.0l7*53»925  *6. 2*3115308  *.  1221730476 

***2617993878  **5235987756  / 

BATA  0180*33/3.1*15316536  *1*01*74426  *.4886*219  *-.57995865  / 

BATA  PERCENT/. 00133  *1  .035  j.*S7  *.*  / 

I  AT  A  SPLAT*:3WAT,PL9V  .979*146.  .1993614*5.078908/ 

I A T A  M1.H2.H3/1346.91  *5l6»4  .59*8*5  / 

P(3) *L3T5 

LBtta  0v£9  CONDITIONS 
1 00  CSNTJNJE 

PEA 3  Date  anO  STATION  POSITION  FRO*  CARD 
PEADI6.1)  I  VR.MBN, IBAy.FLAT.FL9S 

1  FORMAT ( 319,2^10*5) 

jr { IYR.lE.O)  39  T9  900 

KR ! TE ( 6 . ? )  I MR* MON* IDAy,rL*T,pL9N 

2  PBRMAT(//75H  jENEPATE  RECORD  9N  9JTPUT  TA»E  C9STA I S I N3  IONOSPHERIC 

•  F0T2-H-  TA3LES  P9R  /6M  TEAR*, Ig, 8H,  MONTH* >  1 2, 6H,  0AV*,J2,11H,  L 
*ATITU0E*,P  10.8,27m  0E3,  .SN3ITJIE  Or  ST A T I  ON* , F 1 0 . 5* 4H  3E0) 

plat*f.at.or 

* LPN"P-8N-3R 

1ELAQ*0 

ISKIP.O 

IVRM8*;yR#100TM0N 

lM0Oy*MBN*1004lOAV 

RkAi  CBErri:lENT  TaPE 

10  irnweOT.LE.MBSDV.AND.IMBDV.ar.MBNO)  39  T9  29 
20  RE  AD ( I  TP !  L9SD,L?NDT*-C»tr,UM,U91 
ir(liF.lTP)  23,22 
22  I SK I P* 1 
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*!&N0«LJSiD 
*!?NDY«L»MDY 
SB  T8  10 
23  REWIND  1  TP 
Iflag*  irua+i 
JF( IFUAQ.LE* 1 )  38  T8  80 
W»ITE(6,25)  I YR/M9N/ IDAY 

25  f?RMAT(54H  « **C6ErF I C ! £NT8  NBT  rgtJND  BN  TAPE  F8R  YEAR#m8NTH,DAY«* 
*313) 

88  T8  400 

29  IH  IS*IP.EJ.O«AND«  lYRMBtES.LYfHB)  38  T5  10 
C 

C  RE. A 3  SSlAR  DATA 

ir  <  tY*M5.E2.LVi?Y8)  38  T8  55 

RFA0(Si7)  !YM1#  { e“UX  ( Z  J  ,  I  ■  1 ,  1 6  ) ,  1  YM2,  (  FLX  {  !  i ,  I  •  1 7,  in  1 1  YM3, S 1 8#  8  \f 
7  rpR^AT { IAj4X# 16P4.1/l4#15F4»l/I*i?r5»l) 

IP(  !Yf<i.E3*!YMj«AN0' JvM8.raiIYMS#AND*IV*13»CQilYRM8)  08  T9  50 

WRITE  <6>ft) JYR,Y6N 

a  p«RYAT(//39H  *t*ERR8R  in  59*A»  IN^UT  DATA  FB»  YEARp,1Z/11H  AND  N8N 
*TH«ij£) 

8?  T8  400 
50  Ly*Ye»lvR*j 

0 

C  RRE.RARP  SPECIFIC  CBERHCIE^T  8ETB 

55  D8  62  J « 1 >  4 9 
16  6?  I « 1 j  9 

UM( U J) »JV( \t j)*( JH 1 ( 1,J) .JM( I* J) ) *8 1  S/D  loo 
62  C6NT I Nj£ 

88  70  J - 1 i 7 6 
•6  70  1-1/13 

7  0  0  (  I#  J  )  ■wCBEr  (  1*  \t  J  )♦  (  *C6EP  IS/  Ti  J)  *WCBEF  (  3* 1 1  J I  *8  ID  *81? 

C 

C  PREPARE  S9_AR  OATA 

<jO  PLXD'FlX  (  I  0 A Y  ) 

WPITE^lflir.XO/SIFfSlS 

15  F9RMAT(12^  Daj.Y  r L JX . , f 6 , l , 4  In,  12.M8NTH  RUNNING  AVERAGE  6F  8BLAR 
•  PwUX«,F6, 1<20M>  8F  8 1  j  N  S  P  ft  T  N  J'iBER"  »P6»  1  ) 
r.UX^LXD 

Ir(PL4X.LTiDPl  )  PLUX«SIf 
8FLUX'PLUX*SJP 

!P(fLUA.3T, 31301  'LUX'QISO 

GENERATE  25  POINT  PATTERN  ARBuND  8TATI6N 
LB6P  6VER  EARTh  CENTRAL  ANQ^Jg 

EC  A ■ *07 

9*0 

§5  300  I C  A  *  1 »  4 
ECA«ECA*D7 
Sa.  BIN(ECA) 

C  *  ■  CtJB(EOA) 

NAZ-JAZI ICA) 

9 A  Z  *N A Z 
0AZ-PI2/DAZ 
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L99P  AVER  AZIMUTH 
D9  300  1AZ-1/NAZ 

MAM*  1 

A* »AZ*CAZ 
S*Z«  s I N ( AZ ) 

CAZ»  C88(AZ) 

C 

C  C9MPUTE  LATITJOEa  L8N3ITJDE  9P  I9\98»WERlC  P9INT  9LAT<919n 

SMLAT.  S  ?  N (F.AT ) *CA*  C88 ( FL AT ) *8A»CAZ 
CMLAT"  93RT!31*8NLAT*8NJLAT' 

9LAT«  ATANI  JSNLAT/CNLATI 
S9L8N«Saz»8A/Cnl*T 
C9L9N*  SORT(31*S0L9N«8OU9Sl 
9l®N*FL9n*  ATAN(8DL9N/CDL9N> 

C 

C  CBMMJTE  P9SITI9N  BEP£N5FNT  r U*MC T I  9NS  rBR  F0F2  AND  M3000 

P ( 1 1 BBL AT 
P ( 2 ) «9L9N 

CALL  MA3FIN(°aC9M) 

TlDtC«M(2)*c9M(2)AC8M(3)*CB'*(  3> 

C(?)«P!2) 

C? 3 )  "P ( 1 ) 

C ( 1 ) •  ATAN (  ATAN(-CBM( 1)/  SQRT(TAP))/  S3*T(CNLAT>> 

CALL  QK ( <  a  C  #  3 ) 

<K  •  0 

09  85  I I»1j 10a2 
I l ■ <n : ii i 
I  * ■ <  N  !  I 1*1  ) 

09  85  J ■ 1 1 i I  2 
«  ■  KK  t  i 
85  3T(<<)«3(J) 

C 

C  CAMPJTE  ma3\ET!C  .AtiTJDp  9^  I BN8SPHEP 1 C  P9  I  NT 

SMl  «  SNI.AT  •  SaLAT  ♦  CNLAT  ♦  CPLAT  ♦  CASt  9L9N-PL9N) 

CML*  SORT  !  01  ■S,*l*9ML  > 

HLAT«  ATAN ( SML/CML ) 

C 

C  wBOP  9VEP  1“  lBCAL  WSURS 

TL9C»«0h«2 
D9  200  I H ■ 1 #  i  A 
OMR ■ DHR2 

If(  IH«3E.a,aN0. I h • LE • 7  j  0AR.3HR1 

T.9C«TL9C*9HR 

Tl-E«TL9C*fiL9N*!»l2 

TJME«AM9-5(TIME,B1?) 

C 

e  C»MPjTE  TIME  0EbENOenT  FjNcTJ9N8  P9R  rorf  AND  H30CO 

TAT  IME»3l80 

CALL  IIC9JT(6<C9Ti8IT#T) 

CALL  D*8?Cfl I MFr,*( 1 0 > # U# B I T , C9T, DF ) 

CALL  D*0  t  Cfl (  NMF * <Mi 0  #JM<8TT/C9T/DM) 
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C 

C  C8MPJTE  H3000  AND  HEIGHT  BF  HAX»  ELECTRON  OENBITV  HH 

CALL  0<3<(NiHP/3H/0H#H3O00) 

NT  ■  H1*h2*h3000*H3#H3000*h3000 
C 

C  C9MPUTC  F0F2  AND  ADJUST  F0F2  FJR  DAILY  VARI AT18N  WITH  FLUX 

CALL  DK3<(\Ff*3>DF/FQF8) 

L ATI  •  1 
LAT8  ■  i 

ir (HLATiG£,03(uAT2)  )  3fi  T5  91 
LAT2  •  2 

IF(HLAT,QT»33(uAT2) )  38  TS  §1 
LAT1  ■  2 

IT (HlAT.EQ, q3(laT2)  )  Gfl  T!  91 
L ATP  ■  3 

IF(HLA?i.JT,03(LAT2)  )  3ft  T9  91 
L ATI  •  3 

91  CRT  •  CENT  <  l  AT  1 > 

IF(LAT1.EQ.LAT2)  Q9  78  92 

C*dT  «  CNT  *  (CENT(LAT2)-CENfCLATl>  >  •  03  (  L  AT  1 )  •H\m  AT  ) 

#  /  ( H3 ( LAT1 ) »03I LAT2 )  ! 

92  F0"2  ■  PCF2  «  (PER*DFLuX  ♦  CNT) 

I FU( I M/ M ) ■  hH*310*QP5 
If2.r0Fp»3l00*3F5 
IFw(IH#H!*!F-( 1H/H) *10000* I F2 

200  C9 NT  I NUE 
300  CONTI NUE 
C 

C  WRITE  9JT-JT  RECBRD  9F  T ONOSPME^ I C  F0F2-UM  TABLES 

iv*o»  i  vrmogoo*t8n*ioowOav 

WRITE  !  JTP)  I VHOjFlaT^FlBNiF f.UXi I*H 
38  TB  100 
WOO  C 0 N T I NUE 

END  PILE  UTP 
REWIND  JTP 
STBP 
EHD 
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! 


PRBQRAM  1 8MI  <  INPjT/»uTPjTiTAPE5«INt>UTiTAPC6«BuTPuT>TAPE2) 
C9MPUTE3  l9N98»HER!C  PR9F  I LE  PARAMETERS  AND  REFRACTJ9N  CBRRECTI9NB 

UTILIZING  ®re:bmputed  fqf2.hm  tables 

*####  T9  BE  USED  9NLV  F9R  8TRJN3EN?  C9RE  SPACE  AND/9R  RUN  TIME 

REQUIREMENTS^  SINCE  INTERP9LATI9NS  BE  THE  PREC9MPUTED  F0F2.HM 
•  ♦#**  TABLES  CREATE  LESS  ACCURATE  RESULTS  THAN  TH9SE  BBTA I  NED 
*»«**  FR9M  PRB3RAT  J  ©N 

C8NTENT  ftF  C9MM9N  BLBCKS  EXPLAINED  JN  S JBR9UT I NE  REFRCl 
C9MM9N  /LVAL1/  fS/FlAT/FLBN/ElEU/ kit H8/ ED9T/ HDCT/ T IME / 

*ITR>M9N>  I -D A T i  J T P 

C9MM9N  /C8RR1/  DRAN3/ DRATE/ DElEV/FOFP/ HM/YM/ YT/ XK/ T8TN/ T9TNA 
DJHENSI9N  X<  (  3  S 

DATA  QO/2lOOO/Q3bOOiDR/HR  /0»  /10C0.  *3600*  / .017*532925  / 

*.2617993675  / 

JTP«2 
NCM»0 

WRITE(6/26) 

26  F9RMAT ( 1h1  ) 

10  C9NTINUE 

C  Rf  :  AND  PRINT  EVALUATI9N  C9NDITIBN 
RE'-D(5/3)FS<*’LAT<Fl9N 
3  FBRMATtFiO. A#2F10i5) 

IFfFS.LT. 30  )  38  T8  100 

READ«5AP)EL6V<AZ#HS#iD9T*HB9T 
A  FBRHAT(2Fl0.6>FlO.O*2El5t8) 

READ ( 5/ 5 1 IYR/M8N/ ISAY/ T IME 

5  FBRMAT ( 3I3/F10.7) 

WRITE (6,6  >rs/FL AT/F19N/ ELEV/AZ/ MS/ EDBT,  I YR, m9N,  I  DA y, T I  me , HDBT 

6  F9BHATI  jgM  «•  INPUT 

*  1 1  H  FRr3UEN:v»/FlO.<*/  15H  MHZ,  L  A  T  I  T  UD  E  ■  /  F  1 0 . 5 , 
#27W  DEO/  LBN3ITJ0E  8F  STAY  I  AN.  /  F 1 0 .5/ ah  DE3/1H  ELEVAT  I  BN./ r  1 0  •  6/ 

•  IBM  DEO,  *Z IMJTM./F10.6/ 27H  DE3/  HE  I 3HT  9F  8ATELL !TE«, *1 1 i 1, 

* 2 1 H  KM/  E.EVATI8N  RATE., E1B.8/8H  RA0/SEC/6H  YEAR*, 12, 8H,  M9nTH., 

*  I  2 / 6H/  OaV«, 12, iOH/  U. TIME. /F10. 7/5^  *RS,,39X,15h  ALTITUDE  RATE«, 

•  EI5.8/6H  M/SEC  1 
C 

C  CBNVERT  UNITS 

flat  «flat  *dr 
Fl8N*FLBN*0R 

elev«elev*or 

AZ*AZ*DR 

HS"H8*Q1000 

TI“E,T]ME*MR 

c 

C  C9MPUTE  AND  PRINT  1 9N8BPHER I C  DATA 

CALL  REFRC1 
IFdTR'LT.O*  38  T9  10 
XHM«HM/01000 
WRITE16/21!  *NM/F0F2 
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21  FORMAT!  /13H  «♦  OUTPUT  ♦♦//16H  HEIGHT  AT  MAXIMUM  ELECT*#*  DENSITY# 
•l0X/3HHM»#F*t3*30H  KM#_  CRITICAL  FREQUENCY  FO F*«#fT*3#AW  HH|» 
XYM«YM/91000 
XTT«YT/9l^00 

WHITE  <*#«>  TOTNiTOTNA#XVM,XYT#XK 

IE  FiRMATUSH  TOTAL  INTEGRATED  ELECTRON  CONTENT#  VERTICAL  NTffCili*/ 
•WH  C/(M.M),  A  MULSH  NTA**Ett»4#lVH  COLUMN)/ 

#4|H  HALF  THICKNESS  OF  SOYTOmSIDE  IIPARAROLA  YMi#Fl«3# 

•J0W  KM#  8F  TOPSIDE  PARABOLA  >Tl|FS.3#*H  KM/ 

•58H  DECAY  CONSTANTS  FOR  TOPSIDE  EXPONENTIAL  kAYERS#  LOWER  *!•# 
•EUiSiIEh#  MIDDLE  K2P#Cl?*B»iiH,  UPPER  K3a*ElS«5#4N  1/M} 

TEL EVaDELCV* 33400  /OR 

WHI?E(4<EI)TELEV 

23  FORMAT  ( 5#H  IONOSPHERIC  REFRACTION  CORRECTION  TO.ELEYAUON  ANGLE*# 
«tl3* OillH  SEC  OF  ARC) 

WRITE (4#0A)  DRANG 

24  FORMAT  <  43H  IONOSPHERIC  REFRACTION  CORRECTION  TO  RANGE# iOX# 1M»> 
•EI3i6#Sh  M) 

WRITE(*#2B)  5RATE 

25  F0RMAT(54H  IONOSPHERIC  REFRACTION  CORRECTION  TO  RANGE  RATE  *# 
AE13«fe#4H  M/SEC ) 

N0M«NUM*1 

1 F ( NUM • LT • 3 )  33  TO  27 

WHITC(6#tS) 

NQM«C 
GO  TO  10 

27  WR I TE ( 4/ ?• > 

20  FORMAT!//) 

30  TO  10 
100  CONTINUE 
STOP 
END 
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C 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


IBNB8PHERIC  REFRACTJ9N  hBDEL  UTILIZING  PRECOMPUTED  fOFf.W* 
TABLES  F*R  INTERR0LATI9K 

##**•  TB  BE  used  9NLY  reR  STRINGENT  CSRE  SPACE  AND/BR  BUM. TIME 
BcouiHMfMTf  . 

SBBWU riM  tt'Mc'l 

INPDTl  CBMMBN  /EVAL1/ 

9UTPUTI  C9MM9N  /CBRRi/ 

CBMMBN  /EVA U/  rB*FLAT,FL8N>ELEV*AZ*M**ED»T*HD9T*TIMEx 
•IfRlMSN* IDAY* JTP 

FS  ■  TRANSM1S$IBN  FREQUENCY  IN  MmZ 
FLAT  •  BTATIBN  LATITUDE  IN  RADIANS  BF  ARC 

FL9N  •  STATIST  L6NQ IT JOE  JN  RADIANS  BF  ARC  (POSITIVE  EA$T*0  TO  360  0) 

ELEV  «  ELEVATION  SF  SATELLITE  rM  RAlTANS  BP  ARC 

A Z  r  AZIMUTH  BF  SATELLITE  IN  RADIANS  BF  ARC 

HS  *  HEIGHT  BP  SATELLITE  I*  hETERS 

EDBT  »  ELEVATION  RATE  IN  RADIANS  BF  ARC/BECBND 

HDBT  •  RATE  Bp  CHanOF  IN  h£I3HT  BF  SATELLITE  IN  METERS/SECBND 

TIME  •  UNIVERSAL  TIME  IN  RADIANS  BF  ARC 

IYR  •  YEAR  (LAST  2  DIGIT*  BNLY) 

MBN  ■  MBNTh  <«1  THR9U3M  1 2 ) 

I  DAY  ■  DAY  t«l  THRBUflH  31) 

JTP  •  UNIT  ASSIGNMENT  BF  IONOSPHERIC  TAPE  NtTH  FOFB.Hm  TABLES 

COMMON  /CBRRl/  DRANa*DR4TE*OLLEV*F0F!*HM,YM, YT*XK*T9TN*T9TNA 
DRANB*  RANGE  CORRECTION  In  M£TFR» 

DRATI-  RANGE  RATE  CflRRECT IBN  In  METER8/SECBND 
DELEV*  ELEVATION  ANGLE  CBRRECTIBN  IN  RADIANS  BF  ARC 

RAN3E*  RAN3E  RATE*  AND  ELEVATIBN  ANGLE  C6RRECT I BNS  ARE  TB  BE 
SUBTRACTED  FR9M  THEIR  RESPECTIVE  BBSERV AT  1 8N8 
FQF2  ■  CRITICAL  FRE-UENCY  (MmZ) 

HM  •  HEIGHT  at  maximum  ElECTRBN  DEN8ITY  (M> 

ym  ■  half  thickness  bf  the  bbttbmside  bi»arabbla  (mj 

yt  •  half  thickness  bf  the  topside  parabola  <m> 

X<  •  ARRAY  CBNTAINING  DECAY  CONSTANTS  F9R  THE  LBwER,  MIDDLE  AND 
UP*ER  SECT1BN  BF  The  TOPSIDE  EXPBNENTIAL  LAYER  (1/M) 

TBTN  •  VERTICAL  E-ECTRBN  CONTENT  Ie/M**2> 

TBTNA.  ANGULAR  ElFCTRBN  CBNTENT  ( r/M**2 ) 

01  HENS  I  BN  XK (  3) *LYM0C4)* ALAT( 4  >#  Al8N| A) *FLXD<  4)# I Fh { 14# 25# A) 
DIMENI19N  LT(2)/MP(Zl*FT(?)*HT(Z)*FI(2)*Ml(2)*FAtf)*HA(E>*JAZ(4)* 
•KIZ ( A  ) 

DATA  JAZ/1 *4*1*12/*  KA2/1,2,6*  14/*  LYMD/0*0# 0*0/*  NB*NR/4#0/ 

DATA  R* SPLAT* C»L AT* PLBN/6171.EF3* .9799246*. 1993614* 5 .07*90*/ 

ORTA  RM*T0L/667tl00*  *>001726*463  / 

DATA  6D*Cl*G*#3>iDlC0<'lJP5i  Q4P5*  QnM*  RN3*  P  IE*  DR*  HR/O •  *  1  •  *2.  * 

*7*  *100.  *3*8  *4.8  *1.24C10* *49972  * i .f*3lS53U7t  * 

.0174*32915  *  .2617993*7*  / 

I8JI VALENCE  (LT( D.LTll* (LT(2) #LTf)#{HP(i )*MPi)* (MP(2),MP8I 
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C  INITIAL  I  ?  E  "6 N STANDS 

DFlEVOC 

SRAN3«Q0 
.*RAT£«B0 
T0TNR  «0 

t#t*a«co 

jruA&»0 

MYMD»IYR«1000D*MBN*1004.13AY 

c  ;  r- 

C  *EO  FoFg-HY  IMTCRP©t  AT  TON  TABLES  FILE/  SELECT  MBFeR  UT 

1  5B  ?  1 » i / V B 

I F ( SYMQ , Nf . LY YD ( ! ))  8ft  TB  P 
IFIABSi ALAT( I (-FLAT) «3T*TBL)  38  TBS 
1 F { ABS ( AL8N  <  I  )  •F'LONi  >  •3T»TBLi  38  TB  2 

9*  T9  6  .  a-,  L  .. 

2  SBMJA  jp. 

NR»NR+i  ■  ..  T 

I P ( NH«OT*MO )  NR»1 

3  BEAD  ( JTP  )  lYMD(NR)#ALATtN*nAL&N<N*)|PlOtNR)<:<  M*H(LaU.iNR-> 

*#L"1a  lA)  <L.L*1^25> 

IP(EBFjjTF)  4,1 
k  REWIND  jTp 
!FLAG»1PLaU*1 
I P  <  IFLAC.«LC  •  1 J  CjB  T9  3 
WRITE (6,5) 

9  r9RHAT (  63H  **#  roFi-HY  TABLES  FBR  THIS  *TATIBN  AND  DATE  NBT  FBUND 

*  IS  FILE.) 

ITR  •  .1 
RF.TJMN 

6  r.uX.FLXD{ I ) 

FBM  A  Z  I T  J  T  H  AZ,  EARTH  TENTRAL  ANGLE  STATION  TB  SAT.  ECA, 

lesasP^ERi:  ;.at./Lwn.  plat #blbn#  masnettc  lat,  bf  ion*  rbint 

HLAT,  a n"3  _BCAL  the.  TLflC 
I P  (  A  7  •  L  T  •  '1 0  >  *Z«AZoP:8 

S.ATi  SIN(F.AT) 

C„AT.  C“S(F_AT) 

5fu«  S  !  N  t  f LEV ! 

C  r  L  ■  C  3  S  1 V  L  F  V  ) 

S*Z»  sin:  ,\l  ) 
caz»  rest  u ) 

SF«R*CEl/RY 
CF«  SURTC31-SB,»SF) 

B*  •  CEl  *  :r  •  SEL  *  SF 
CA  •  bEl  •  :f  4.  CEL  «  SF 

ECA«  atan ( sa/;a)/DR 
SNuAT*SL A^#t A*CwAT*SA# CAZ 
C  Nta  AT  ■  S3'ma>SNLAT*S\LAT) 

BuAT«  ATAMf SN.AT/CNLAT) 

SDL8N»SAZ«8A/c:NwAT 
CDLBN«  S3RT(0HSDL6N#sDLBN) 

B.8N«pLHN*  ATAN(SDLBN/cnwhN) 

SML  •  8NL  A T  *  flPLAT  ♦  C  N  L  a  T  •  C^L  AT  *  CBB(  BL8N*PLBN> 

CYL*  SURT  ;  U-SYL*BML> 

WLAT«  ATAN(SYl/CYL) 

TlBC»THE*(J.8sj*R!8 
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TC8C*A*8D(T<.9CiM>  .  :  :-v  , 

TC8CAL»tL8C/Ht 

InterpBlatc  rort«*4W  TAfc.ES  — 

CSM^UTE  INDlCf®  V7t*LT£,  INCRCMENT.  Bv?  if*  L$CAL.  JVTf*F|U 
X;T«?L8CAw/3^3i  '  v  w 

■■■'■ .  uft.*xk.t  ,  •  ? .  ‘ 

&LT-L.T1 
DwT*XUT.OLT  • 

I*(XLT.3E»33*5)  LfltUTm  ! 

I*  (  XLT  *  3E04BB 1  LTi»LTtM 
imfi-LT-a-M *u tuatvi y  a®  ts  id 

OtT«8LT#ot  -  -  .  ■  -  I 

1*<DLT.»3E01)  5LTOlT«Oi 
1C  19 i LT1.3T.14)  ,T1«1 
L?*»LTWi 

I «“  <  uT2 •  3T •  i **  >  .T2*l 

CBMp JTt  EARTH  CENTRAL  ANGLE  INDEX  1*19,  INCREMENT  DAL** 

At F -EC A /Q7 + j 1 
lAL'-ALr 

IK  IALF..QT.A)  !  AuF- A 
0ALF»ALF»  FLBAt i I ALF ) 

<!  ■  1 

CBMPjTE  AZImjTH  INDICES  MPl/NRS#  INCREMENT  DELAZ  ! 

20  NAZ  ■  JAZ ( I AlF !  ; 

M»1.<AZ(IAlF) 
xr  c  MP1 . 3T« 1 >  3?  T8  30 
.DCLAZ-00 

MFB-1 
39  T5  60 

30  DAZ IM.P !?/  F.BaT (NAZI 
AZIM.QC 

DB  AO  L50p ■ 1 / N  a  Z 

HF2-MP1 

NA1-MP1-M 

IK  L8»F,EQ,NAZ  !  mp  1  «»AP  1  -NAZ 

AZ I M» AZ I H*DAZ I m  ! 

I F ( AZ  I M» OE • AZ  )  38  TB  50  ! 

AC  C8nT INuE 

5C  DtLAZ«( AZI^-AZ  )/DAZ  IN 

*C  C9NTINUE  I 

I 

DB  SO^IPT^l ^2^  ^  TI1E  BBINTS  NPi#MPt  TB  3£T  rj#Ht 

MPT-MP< XPT I 
08  70  1.-1#  P 
LTM-LT1L) 

IH1-1FH(lTm#mPt#  n/10©o0 

MT(LJ-  FlBaT( I -1 ) -510C 

IF  5  -!FM(lTm#mPt#  I!  •  HUIOCCO 

70  FT  f  L  )  •  FlBa?(  Ifl  >/01W  i 


n  n r?  r» n  or»  rin 
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fiurti  *rru  trYtitwT  1 i >  >*dct 

10  H? ( I*T ) «HT( 1 >4 < MT {t >fHT< J) )*9LT 
C  . 

C  INTIMfll.ATI  I-N  ?§  «fT  FA^WA 

r»uvi»rim*(r2tf)*ri{x>  u^eux 
>♦<*!*»>•«! n>  l*KU*I 
2F(KltCQ»ll  90  TO  100 
Kt*» 

IIL?*Ialfai 

iraALF.ar.4)  39  to  »o 

30  T8  80 
90  PK*>»FA<1) 

H»(t)*HA(l> 


INTERPOLATE  IS  EARTH  CENTRAL  AN3LE  T9  SET  F0F»#HM 
100  F0F8.FA{l)f  {•TAjai.FAfl  J  J«DALP 
H1«WA(l)*(WA(f  )»HA(U  )  *DALF 

COMPJTE  iECOND  PART  9F  PROFILE 
CALL  M8FLI(8LAT#»L8N^l,Tmt*  !DAV*M8NjFLUX#FQF8aHMa MLAT# 

«  VM,yT*XK#RRM>XNTNM! 

IF(XST^M,LE»30)  38  T8  1A0 

COMPUTE  ELEVATION  AN&LE  C6RPECTI0N  OELEV 
FRAT»!F0F2/FB)*»E 

CALL  84 TA<r*AT,XNTNM#HS,KM,VM,8IL,CEL* OELEV) 

C0MP*U  VERTICAL  AND  AN3JLAR  ELECTRON  CONTENT  T8TN/T8TNA 
COMPUTE  9AN3E  CORRECTION  DRANQ 
RAT«{R/{R*  HM) > «*B 
DEN0»B1-RAT«CEL*CEl 
OEN«  BORT(DENI) 

T0TN»kK'TNM#3NM*F0F8*»E 
TOTNA.TOTN/OEN 
DR*N3»FRaT*RV3*XNTNM/0EN 

C  COMPUTE  PAN3E  RATE  CORRECTION  ORATE 

DRATE»3RAN3«f.00T*RAT*lEL*CFL/DEN2 
ORATE «DRATJ  .'■RAT»RN3*HD0T*RRM/DEN 
UO  CONTINUE 
RETURN 
END 
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